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FOREWORD 


I 


Icing  of  unprotected  aircraft  components  is  one  of  the  major  problems  engineers  are 
faced  with  during  the  development  phase  of  an  aircraft.  Ice  accretions,  as  well  as  their 
shape,  have  to  be  predicted  in  order  to  investigate  their  effect  on  aerodynamic  flight  safety 
and  performance,  and  to  decide  on  the  need  for  provision  of  anti-or  de-icing  systems. 
Typical  components  where  hazardous  accretions  may  occur  are  wings,  tails,  engine  inlets, 
hubs,  antennas,  propellers  and  rotor  blades. 

The  Fluid  Dynamics  Panel  of  AGARD  held  a round-table  discussion  of  the  fluid 
dynamic  aspects  of  icing.  This  took  place  following  the  FDP  Symposium  on  Unsteady 
Aerodynamics  in  Ottawa,  Canada,  on  the  30th  September,  1977. 

Seven  contributions  two  from  Germany,  two  from  Canada,  and  one  each  from 
France,  UK  and  USA  - were  presented,  covering  abroad  spectrum  of  topics  such  as  icing 
parameters,  icing  of  wings,  tails,  inlets,  helicopter  blades  and  icing  test  facilities  within  the 
NATO  countries  of  Europe,  the  USA  and  Canada.  These  papers  have  been  collected  in  this 
AGARD  Advisory  Report.  It  gives  an  excellent  state-of-the-art  survey  of  this  subject.  The 
subject  of  engine  and  inlet  icing  has  been  treated  separately  in  the  Spring  1978  Meeting  of 
the  AGARD  Propulsion  and  Energetics  Panel,  and  is  reported  in  AGARD  Conference 
Proceedings  CP-236. 

It  is  hoped  that  by  this  report  further  investigations  into  the  establishment  of  atmos- 
pheric icing  parameters  will  be  initiated,  that  the  full  survey  of  icing  test  facilities  will  allow 
a proper  choice,  and  that  the  described  procedures  and  data  on  wing  icing  will  be  helpful  to 
predict  ice  accretion  on  future  aircraft. 


B.LASCHKA 

Chairman,  Programme  Committee 
(AGARD  Fluid  Dynamics  Panel) 
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SOME  RECENT  RESULTS  ON  ICING  PARAMETERS 


W.  Kleuters  and  G.  WOlfer 
MESSERSCHMITT-BOLKOW-BLOHM 
2103  Hamburg  95 
West  Germany 


SUMMARY 


Meteorological  icing  data  on  forty-seven  encounters  in  cumulus  cloud  over  Western  and 
Central  Europe  are  here  analysed  statistically,  and  compared  with  similar  data  quoted 
in  FAR  25  Appendix  C. 

The  data  for  FAR  25  Appendix  were  obtained  by  the  "Rotating  Cylinder"  method,  a 
measuring  device  with  a capacity  limited  to  relatively  low  liquid  water  content  and/or 
low  temperatures.  To  overcome  these  limitations,  the  liquid  water  content  was  measured 
in  these  tests  by  employing  a liquid-nitrogen-cooled  Test  Rod,  enabling  measurements 
of  liquid  water  content  of  up  to  3 g/m  . 

It  is  thus  evident  that,  in  the  temperature  range  above  -15  °C,  the  MBB  analysis  yields 
values  of  liquid  water  content  considerably  higher  than  quoted  in  FAR  25. 


1 . INTRODUCTION 

The  design  of  ice-protection  equipment  requires  some  defining  of  the  meteorological 
conditions  in  which  the  system  is  expected  to  operate.  The  conditions  are  governed  by 
liquid  water  content,  mean  effective  droplet  diameter,  and  ambient  temperature. 

The  investigation  of  icing  parameters  began  with  the  development  of  heating  wings, 
leading  by  1948  to  an  accumulation  of  some  1000  individual  measurements  of  water  content 
and  effective  droplet  size,  obtained  in  two  hundred  forty  separate  encounters  with 
icing  conditions.  After  review  and  discussion  by  NASA  scientists,  the  Technical  Note 
1855[1]was  published,  providing  industry  with  tentative  standards  to  fill  an  immediate 
need.  These  recommended  design  conditions  were  adopted  by  the  Civil  Aeronautics 
Administration  and  are  still  in  effect  today,  embodied  in  FAR  25  Appendix  C. 

The  design  envelopes  for  1 ) stratiform  or  layer  clouds  (the  continuous  maximum)  and 
2)  cumulus  clouds  (the  intermittent  maximum)  are  shown  in  Figures  1 and  2 respectively, 
in  the  form  of  curves  drawn  from  combinations  of  liquid  water  content  and  effective 
droplet  diameter,  for  selected  values  of  temperature. 

Examination  of  the  statistical  distributions  of  the  parameters  indicated  that  the 
probability  of  all  three  variables  - as  represented  by  a single^point  in  Figure  1 
and  2 - being  exceeded  at  any  one  time,  was  of  the  order  of  10-J[21.  Considering  the 
data  were  obtained  from  one  hundred  sixty-seven  encounters  in  layer  clouds,  and  seventy- 
three  in  cumulus  clouds,  the  statistical  extrapolation  to  a probability  of  10-J  was 
indeed  questionable. 


2.  THE  ROTATING  MULTI-CYLINDER  METHOD 

The  majority  of  flight  measurements  of  liquid  water  content  and  effective  droplet 
diameter,  upon  which  the  curves  in  Figures  1 and  2 depend,  were  made  using  the 
"Rotating  Cylinder"  method.  A set  of  five  cylinders  (Figure  3)  of  different  diameters 
(0.125"j  0.5"i  1.25"i  3.0"  and  4.5")  wbb  exposed  normal  to  the  airstream  while  slowly 
rotating.  Theoretical  relations  involving  droplet  size,  cylinder  diameter,  collection 
efficiency,  and  airspeed,  were  used  to  derive  average  liquid  water  content  and  effective 
droplet  diameter (3) . 

The  measurement  of  liquid  water  content  and  droplet  diameter  by  this  "Rotating  Cylinder” 
method  is  subject  to  the  following  sources  of  error:  1)  blow-off,  2)  run-off 
3)  evaporation,  and  4)  variation  in  water  concentration  in  the  vicinity  of  the  cylinders, 
due  to  disturbing  Influences  from  the  aircraft. 

In  a statistical  analysis,  however,  reasonably  reliable  results  may  be  obtained  provided 
that  the  errors  are  randomly  distributed  and  small  compared  with  the  quantity  being 
measured.  But  the  problem  of  blow-off  merits  some  further  consideration  since  it  Imposes 
a systematic  error,  as  with  all  instruments  measuring  rate  of  icing,  the  "Rotating 
Cylinders"  are  absolutely  limited  by  latent  heat  effects.  The  following  equations  show 
the  heat  balance  at  the  surface  of  a cylinder  assuming  that  the  surface  temperature 
-i£  is  0 C and  the  total  temperature  iQ  is  below  zero. 
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The  various  contributing  factors  affecting  the  heat  balance  at  the  stagnation  point  are: 

a)  The  density  of  heat  flux  due  to  convection:  t 

I 

q1  s <&  ^ ) with  ^d'1  ♦ r (1 ) 

where:  •JL  = Static  air  temperature 

r * Recovery  factor  ~ 1 
cp  = Specific  heat  of  the  air 
*-  = Heat  transfer  coefficient 
V®  = Air  speed 

b)  The  density  of  heat  flux  due  to  the  temperature  rise  of  the  water.  Since  the 
temperature  at  the  surface  ist  0 °C,  a fraction  f of  the  supercooled  water  changes 
into  ice,  and  if  L'  is  the  heat  of  fusion  of  ice,  then: 


q2  --mw(cw^«fL')  with  mw-E-W-V®  (2) 

where:  cw  = Specific  heat  of  water 

f = Freezing  fraction 
h = Collection  efficiency 
W = Liquid  water  content 

c)  The  density  of  heat  flux  due  to  evaporation: 


q3 


0.62^  • ^ 

c p B 


(3) 


where:  P - Vapour  saturation  pressure  of  water  or  ice  at  the 

surface 


P « Partial  water  vapor  pressure  at  the  edge  of  the 
boundary  layer 

L"  « Heat  of  vaporization  of  both  water  and  ice 
B “ Barometric  pressure 


d)  The  density  of  heat  flux  due  to  kinetic  energy  of  the  droplets: 

V2 

q4 ^ mw 

If  the  surface  of  the  cylinder  is  neither  heated  nor  cooled: 


(4) 


q 1 ♦ <l2+  q3  '^4"° 

or : , , ( 5 ) 

For  any  combination  (W,Va,,-^)  the  freezing  fraction  f can  be  calculated  by  the  equation 
(5).  The  limit  for  given  values  of  (W. , ,-J^ | corresponding  to  total  freezing,  i.e. 
f “ 1,  is  called  the  Ludlam-limit,  which  is  of  major  importance  when  employing  devices 
that  measure  liquid  water  content  of  a cloud  by  means  of  the  icing  speed. 

Thus  the  "Rotating  Cylinders"  cannot  operate  precisely  when  the  Ludlam  limit  is 
exceeded,  i.e.  when  f < 1,  because  part  of  the  impinging  water  is  blown  off,  leading 
to  smaller  measured  values  of  liquid  water  content.  This  built-in  error  increases  with 
increasing  temperature,  airspeed,  or  liquid  water  content.  Figure  4 shows  the  Ludlam 
limits  for  three  cylinder  diameters  as  used  in  the  "Rotating  Cylinder"  device. 

The  results  indicate  quite  clearly  that  "Rotating  Cylinders"  should  be  used  only  for 
conditions  of  small  liquid  water  content  and/or  low  temperatures. 

Since  the  highest  liquid  water  content  occurs  in  cumulus  clouds  - about  three  times 
greater  than  that  in  layer  clouds  - it  was  expected  that  the  largest  errors  would  show 
up  in  the  intermittent  maximum  curves  of  Figure  2.  Thus  efforts  were  concentrated  on 
measurements  in  cumulus  cloud. 


h 
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3.  MEASURING  EQUIPMENT  USED  BY  MBB 

a)  Liquid  Water  Content  t 

To  overcome  the  limitations  associated  with  the  "Rotating  Cylinders"  at  high 
liquid  water  content,  MBB  developed  a special  cooled  Test  Rod  [41  , [51  . 

This  device  consisted  of  a hollow  steel  cylinder  30  cm  long,  and  only  1 cm  in 
diameter, to  ensure  a rate  of  catch  nearly  independent  of  droplet  diameter 
(Figure  5).  The  cylinder,  with  an  internal  flow  of  liquid  nitrogen,  was  exposed 
normal  to  the  airstream.  Three  thermocouples  were  fitted  axially  at  the  surface 
to  provide  continuous  checking  of  the  required  -160  °C  surface  temperature. 

Thus,  the  latent  heat  of  fusion  was  transferred  and  liquid  water  content  of 
up  to  3 g/mJ  could  readily  be  measured,  except  where  thickness  of  ice  at  the 
stagnation  point  exceeded  5 mm.  The  Test  Rod  was  extensively  calibrated  in  the 
wind-tunnel  at  Artington  (U.K.),  under  the  test  conditions  given  in  Figure  6, 
before  any  actual  flight  test  measurements  were  made. 

To  check  the  reproducibility,  more  than  a hundred  conditions  of  (W,  ) were 

tested.  It  was  found  that  maximum  error  was  in  the  order  of  10  * over  the  test 
range . 

b)  Droplet  Size  Measurement 

To  measure  droplet  size,  the  "Oil  Slide"  technique  was  employed.  A slide  was 
prepared  by  smearing  a strip  of  glass  with  a small  amount  of  a special  oil. 

The  slide  was  next  inserted  into  a holder  and  retracted  into  a special  sampling 
pole.  This  pole  was  then  extended  into  the  airflow,  with  the  slide  surface 
vertical  and  facing  forward.  The  slide  was  now  exposed  for  a short  time  in  cloud, 
removed,  transferred  to  a microscope  and  photographed. 

Figure  7 shows  a typical  photograph  of  a droplet  sample. 

I 

c)  Total  Temperature  Measurement 

The  total  temperature  was  measured  with  a Rosemount  temperature  sensor, Model  102, 
suitable  for  operation  under  icing  conditions.  The  error  due  to  de-icing  was 
approx.  +0,5  ”c,  which  could  be  considered  acceptable. 

d)  Data  Recording  and  Test  Procedure 

Flight  tests  under  natural  icing  conditions  require  a very  high  standard  of 
coordination  in  all  relevant  data  acquisition  instrumentation.  Therefore,  airspeed, 
altitude,  ambient  air  temperature,  exposure  time  of  the  Test  Rod,  and  a time  code 
- which  was  also  registered  on  all  photographs  - were  simultaneously  recorded  on  a 
tape  recorder. 

The  aircraft  used  was  a TRANSALL  C-160,  with  a maximum  take-off  weight  of  some 
100.000  lbs.  The  tests  took  place  over  Western  France  and  Southern  Germany,  at 
altitudes  ranging  from  3.000  ft  to  20.000  ft. 

A factor  which  may  have  had  an  influence  on  the  reliability  of  the  data,  was  the 
systematic  bias  introduced  by  obtaining  measurements  during  flights  in  which  icing 
conditions  were  deliberately  sought,  instead  of  being  accidently  encountered.  The 

net  effect,  however,  is  expected  to  be  negligible,  because  results  have  been  compared  ' 

with  those  of  NASA,  obtained  when  faced  with  the  same  situation. 


4.  AMOUNT  OF  DATA  AVAILABLE 

Thirty-four  flights  under  natural  icing  conditions  have  been  made,  during  which  a total 
of  75  icing  areas  were  encountered  An  icing  encounter,  the  basic  unit  of  data,  consists 
of  all  measurements  made  during  flight  through  a single  area  of  icing  conditions.  An 
icing  area  is  defined  as  one  that  is  separated  by  at  least  40  miles  of  ice-free  conditions 
from  the  next  observed  area  of  icing,  or  by  a least  6000  ft  of  altitude  from  another  such 
area. 

The  icing  clouds  have  been  classified  as  either  cumulus  or  stratiform,  icing  encounters 
having  been  in  forty-seven  of  the  former  clouds,  and  twenty-eight  of  the  latter.  Because 
of  insufficient  data,  results  from  the  stratiform  cloud  encounters  are  omitted  in  this 
report. 


5.  ANALYSIS 

For  direct  comparison  with  FAR  25  Appendix  C,  the  analysis  was  carried  out  strictly  in 
accordance  with  TN  2738.  For  convenience,  the  main  features  of  this  analysis  are  cited 
here: 

Exceedance  Probability 

The  three  fundamental  icing  variables  used  to  define  any  icing  condition  are  the  liquid 
water  content  W,  the  mean  effective  droplet  diameter  Def^,  and  the  ambient  temperature  . 
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Thus,  any  particular  icing  condition  may  be  represented  by  selected  values  of  these 
variables,  i.e.  W,  D ff,  and  x% . The  probability  of  any  random  icing  encounter  of  standard 
extent  providing  simSlCaneous  values  of  each  of  the  icing  variables  equal  to,  or  greater 
than,  each  of  the  selected  values  (i.e.H  ^ W.  Deff1  'l^l—  I will  be  termed 

exceedance  probability  Pe>  related  to  the  icingeconditioAs  W1 , Def  f i » abd  , 

The  exceedance  probability  represents  the  probability  of  a simultaneous  occurrence  of 
three  events  which  are  not  necessarily  indepedent.  In  order  to  evaluate  P , use  is  made 
of  the  theorem  of  compound  probability  [ 6 ] which  states  that  the  probability  of  a 
simultaneous  occurrence  of  three  events  A,  B,  and  C,  is  equal  to  the  probability  of  the 
occurrence  "A"  times  the  conditioned  propabillty  of  the  occurrence  "B"  under  the  condition 
that  “A"  will  occur,  times  the  conditioned  probability  of  the  occurrence  "C"  under  the 
condition  that  both  "A"  and  ”B"  will  occur. 

This  theorem  may  be  expressed  symbolically  as  follows: 

Pe  = PA  . PB(A) .PC(A,B)  (6) 

If  it  is  considered  that  the  probability  of  exceeding  certain  values  of  W,  D andlx&l 
is  analogous  to  the  probability  of  the  occurrence  of  events  A,  B,  and  C,  as  expressed  in 
equation  (6),  then  the  exceedance  probability  P may  be  represented  by  the  product  of  three 
probability  factors  Q,  R,  and  S,  according  to  the  equation 

Pe  = Q.R.S  (7) 

where 

Q llx^l)  is  a term  that  defines  the  probability  Q as  a function  of  x%, , and  expresses  the 
probability  for  a given  value  of  x^£  for  a random  icing  encounter  having  a 
temperature  depression  below  freezing,  equal  to  or  greater  than|t^| . 

RlWjxfcl!  is  a term  that  defines  the  probability  R as  a function  of  W where  is  specified, 
and  expresses  the  probability  for  a given  value  of  W for  a random  icing  encounter 
having  a liquid  water  content  equal  to  or  greater  than  W,  under  the  condition 
that  the  temperature  depression  below  freezing  is  equal  to  or  greater  than 
selected  values  of  It^l. 


SIDeff  ,W  |i£,|)ls  a term  that  defines  the  probability  S as  a function  of  D when  W and 

are  specified,  and  expresses  the  probability  for  a given  valai  of  D for  a 
random  icing  encounter  having  a value  of  droplet  diameter  equal  to  or  greater 
than  D under  the  condition  that  the  liquid  water  content  is  equal  to  or 
greate?rthan  W,  and  the  temperature  depression  below  freezing  is  equal  to  or 
greater  than  1x^,1. 


Figure  8 represents  a graphical  plot  with  the  variables  W,  D 
given  icing  condition  W. , Deffi>lt?Sli»  defined  by  selected  val 
appears  as  a single  point  0.eon1the  plot. 


, and  lx£  I , as  axes.  Any 
ues  of  these  variables. 


Corresponding  to  the  point  O. , there  is  an  exceedance  probability  that  a random  icing 
encounter  will  provide  simultaneous  values  of  each  of  the  icing  variables,  equal  to  or 
greater  than  each  of  the  selected  values  W. , D f f 1 , and  xPj,  . . Referring  to  equation  |7), 
the  value  Pgl  would  be  given  by  the  equation: 


Pei  “ <Wl 

Thus,  when  P regarded  as  a constant,  P . , equation  (7)  defines  a surface  passing  through 
0,,  represented  by  all  the  combinationseof  W,D,f,  andt^ , which  have  the  same  probability 
of  being  simultaneously  equalled  or  exceeded. 


The  preceedlng  theory  was  slightly  modified  by  the  authors  of  TN  2378, for  that  portion 
of  the  equiprobability  surface  that  represents  low  values  of  liquid  water  content  and 
large  values  of  droplet  diameter,  resulting  in  two  portions  of  the  equiprobability  surface. 
No  reason  was  found  to  modify  the  theory  in  such  a manner  in  this  report. 


6.  STATISTICAL  CALCULATIONS  AND  RESULTS 
Q - Q (Ix^l) 

The  “N"  maximum  values  of li^l measured  within  N icing  areas  were  arranged  in  a decreasing 
sequence  of  values  of  1x^,1,  beginning  with  the  largest  values  oflx^land  the  ordinal  number 
m ■ 1.  If  a certain  temperature  occured  p times  in  p icing  encounters,  then  only  one  value 
l«  was  used  with  m ” "V  m (m+p-1 ) , where  m is  the  first  and  m+p-1  the  last  ordinal  number. 

From  these  m values,  Q « Qllxfcltwas  determined  by  the  equation 

3 m - 1 

Q “ 3 N'  V 1 (8) 

where  m is  the  ordinal  number  and  N the  number  of  icing  encounters . 

R - R(W,  IxJil) 

Based  on  the  preceedlng  theory,  twelve  temperatures -Jtn  between  -7  °C  and  -17,5  °C  were 
selected.  For  each  temperature,  the  maximum  liquid  water  contents  were  arranged  with 
decreasing  W and  R(W,I^|) determined  by  equation  (8)  above  as  follows: 
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where  Nn  is  the  number  of  icing  encounters  with  . 

S = S <DefffW,h%|) 

For  selected  values  of  W andhfal,  the  function  S(D  ff,W,l$»l)wa8  determined  in  the  same 
way  as  QlloTilland  R(W,  I >1^,1)  . However , it  could  be  demonstrated  that  S was  independent  of 
R and  Q.  With  these  three  functions  - Q,  R,  S - the  equiprobability  surface  for  selected 
exceedance  probabilities  Pg  could  be  calculated  by  the  formula 

Pe  = Q R S 


Results 


Figure  9 shows  the  fungtion  Q(|\fel)on  normal  probability  paper.  The  curve  intersects  the 
li&laxis  at  approx.  -6  °C,  due  to  aerodynamic  heating  at  velocities  between  100  m/s  and 
130  m/s  of  the  test  aircraft.  For  purposes  of  design  of  ice  protection  systems  in  modern 
aircraft,  however,  temperatures  above  -6  °C  are  generally  of  no  importance,  because  this 
aerodynamic  heating,  even  during  holding  and  descent,  is  greater  than  6 °C. 

The  functions  R|W,I-#dI)  forlt$&>|2  7 °C,  2 10  °C,  and  ^17  °C  are  shown  in  Figures  10,  11,  and 
12  respectively  . The  function  S,  which  as  already  mentioned,  is  independent  of  W and  tSS  , 
is  shown  in  Figure  13. 

For  direct  comparison  with  the  plots  used  in  FAR  25  Appendix  C,  a graphical  representation 
was  chosen,  with  liquid  water  content  and  mean  effective  droplet  diameter  for  different 
temperatures,  for  three  values  of  exceedance  probability  - viz.P  = 10  ,10-  ,and  10~ 

(Ref.:  Figures  14,  15,  and  16).  e 

The  limitations  noted  above  for  "Rotating  Cylinders"  show  up  very  clearly  at  high  liquid 
water  contents  and  high  temperatures.  From  Figure  15  it  appears  evident  that,  for 
temperatures  higher  than  -15  °C,  the  anticipated  exceedance  probability  of  10~  -according 
to  the  data  -1b  °C  from  TN  1855  and  TN  2738  - seems  to  be  in  fact  considerably  greater. 
Figures  14  and  15  indicate  that  an  exceedance  probability  of  between  10~  and  10~  may 
be  assumed.  Thus  the  results  obtained  from  these  tests  tend  to  be  in  general  agreement 
with  the  Ludlam  theory. 

For  low  temperatures  and  small  liquid  water  content,  where  the  "Rotating  Cylinders "were 
expected  to  measure  more  satisfactorily,  test  values  obtained  were  even  lower  than  those 
of  FAR  25  Appendix  C.  There  may  be  several  reasons  for  this: 

a)  Statistical  error,  due  to  extrapolating  from  limited  sampling  of  icing  encounters 

b)  Regional  differences  in  meteorological  conditions  between  US  West  Coast  and  Europe 

c)  Vatiations  in  liquid  water  content  in  the  vicinity  of  the  test  aircraft  fuselage 

The  statistical  errors  for  a confidence  limit  of  68  % have  been  estimated  and  are  shown 
in  Figure  16  for  1 7®  = -7  C (dotted  lines).  However,  the  upper  limits  thereof  should  be 
regarded  with  considerable  caution,  because  statistics  are  no  respecters  of  cloud  physics. 

Variations  in  liquid  water  content  in  the  vicinity  of  the  aircraft  fuselage  have  been 
calculated,  and  account  for  an  increase  in  the  measured  liquid  water  content  of  approx. 20%. 
This  value  has  been  taken  into  consideration  in  Figures  14,  15, and  16. 

Finally,  it  may  be  noted  that  all  the  values  of  liquid  water  content  have  been  corrected 
to  2.6  nautical  miles  to  agree  with  the  distance  factor  in  FAR  25  Appendix  C. 


7.  CONCLUSIONS 

In  conclusion,  when  considering  criteria  from  - say  - the  last  twenty  years,  in  relation 
to  the  design  and  certification  of  any  ice  protection  equipment  for  transport  aircraft, 
it  is  recommended  that  all  available  current  criteria  on  conditions  of  higher  liquid  water 
content,  and  temperatures  above  -15  °C,  be  taken  into  consideration  in  future  projects. 

Further  worldwide  observation  and  measurement  of  these  meteorological  conditions  should 
be  urgently  initiated, .employing  instrumentation  and  technique  which  are  not  limited  to 
measuring  small  liquid  water  content. 

This  would  result  in  more  solid  foundations  for  future  statistical  analysis  of  the 
phenomena . 
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SUMMARY 

Icing  detection  on  a helicopter  presents  difficulties  because  of  the  wide  relative  speed  range  from 
hover  to  maximum  cruise.  High  detection  sensitivity  is  necessary  at  any  flight  speed  because  of  the  high 
susceptibility  of  the  rotor  blades  to  the  effects  of  even  quite  small  ice  accretions.  An  icing  detector 
is  described  that  largely  overcomes  the  problem  of  speed  dependence  and  thus  permits  the  determination  of 
icing  rate  and  cloud  liquid  water  content  with  reasonable  accuracy. 


INTRODUCTION 

The  susceptibility  of  fixed  wing  aircraft  to  cloud  icing  is  well  known,  but  perhaps  it  is  not  gener- 
ally appreciated  to  what  degree  the  flight  safety  of  some  helicopters  is  affected  in  very  short  time 
periods.  This  is  particularly  true  for  light  helicopters  (Ref.  1),  but  Bradley  in  a recent  paper  (Ref.  2) 
has  shown  that  this  can  be  true  for  larger  helicopters  as  well,  and  he  documents  a case  where  Intolerable 
rotor  loing  conditions  were  reached  within  a minute  of  the  onset  of  icing.  If  helioopters  are  to  fly  IFR, 
and  that  implies  flying  in  icing  conditions,  then  instrumentation  capable  not  only  of  quantifying  the 
icing  environment  but  also  of  providing  an  almost  instantaneous  Indication  of  having  encountered  an  icing 
condition  is  essential.  The  quantitative  information  provided  by  such  instrumentation  (i.e.  the  loing 
severity)  is  essential  not  only  to  inform  the  pilot  so  that  he  may  take  the  appropriate  action,  but  also 
for  the  proper  control  of  any  icing  protection  systems  that  may  be  fitted  to  the  helicopter. 

What  sort  of  information  should  such  instrumentation  provide?  The  critical  components  on  the  heli- 
copter are  the  engines  and  the  rotor  blades.  The  icing  qualification  requirements  of  the  engine  should 
ensure  that  no  icing  problems  will  be  encountered  there  — other  than  the  danger  that  oould  result  from 
the  ingestion  of  shed  ice  or  slush.  Hopefully  the  more  robust  engines  used  in  the  new  generation  of 
helicopters  will  largely  obviate  this  danger.  That  leaves  the  rotor  blades,  and  clearly  the  information 
that  is  required  in  regard  to  them  is  the  rate  at  which  the  ice  is  building  up  at  some  representative 
radius.  The  shape  of  the  ice  is  also  important,  but  this  is  a function  largely  of  air  temperature,  which 
is  readily  measurable,  as  well  aa  of  the  liquid  water  content.  The  presence  of  snow  or  ice  crystals  in 
the  supercooled  icing  cloud  also  affeots  the  loing  rate  and  the  shape  of  the  ioe  accretion,  the  implica- 
tions of  which  have  not  been  fully  appreciated  until  recently.  It  seems  probable  that  the  rapid  rate  of 
performance  deterioration  at  high  air  temperatures  as  reported  by  Bradley  (Ref.  2)  is  the  result  of  such 
mixed  icing  conditions.  The  measurement  of  such  mixed  conditions  is  a formidable  taak  for  which  at  pres- 
ent no  means  exists. 

ICE  DETECTION  SYSTEM 

An  instrument  has  been  developed  to  provide  data  on  the  supercooled  droplet  cloud  conditions  which 
represent  by  far  the  most  cession  type  of  icing  conditions.  Because  of  the  wide  relative  operating  speed 
range  of  a helicopter  compared  to  that  of  a fixed  wing  aircraft,  the  detection  of  icing  conditions  has 
presented  special  problems.  The  droplet  catch  rate  and  hence  the  loing  rate  are  a function  of  velocity, 
and  thus  the  response  of  an  ice  deteotor  probe  mounted  on  the  fuselage  of  a hovering  helicopter  is  essen- 
tially nil,  while  at  the  same  time  the  essential  components  for  sustaining  flight  (i.e.  the  engines  and 
the  rotor  blades)  are  experiencing  substantially  the  same  icing  rate  aa  in  cruising  flight.  To  maintain 
a constant  response  to  a given  set  of  icing  conditions  regardless  of  flight  velooity,  some  means  of  estab- 
lishing a constant  relative  velooity  between  the  ioe  sensing  element  end  the  droplet  laden  atmosphere  is 
needed. 

In  the  instrument  which  is  referred  to  aa  the  Dynamlo  Ioe  Deteotor,  the  relative  velooity  is  obtained 
by  moving  the  droplet  laden  air  over  a stationary  ioe  sensor,  rather  than  moving  the  sensing  means  through 
the  ambient  air.  This  is  achieved  by  placing  an  ioe  collecting  probe  in  the  secondary  Induced  flow  of  an 
air  ejector.  The  air  ejeotor  forms  the  body  of  the  instrument  (Fig.  l)  and  comprises  a smooth  entry  for 
the  induced  ambient  air  and  within  this  entry  is  the  primary  air  noxxle  in  the  form  of  aa  annular  wall  Jet. 
A short  cylindrical  mixing  length  behind  the  annular  Jet  acoommo dates  the  loe  collecting  probe  and  the 
sensing  means.  The  primary  and  secondary  air  flows  are  exhausted  to  the  atmosphere  through  a diffuser. 

The  air  required  for  operation  of  the  ejector  noxxle  is  normally  delivered  from  the  engine  compressor 
bleed  at  an  operating  pressure  of  about  50  pslg  (31*5  kPa),  with  an  air  mass  flow  of  about  0.5  lb/ain. 

(0.23  kg/min.).  By  ducting  the  hot  bleed  air  up  the  leading  edge  of  the  strut  and  into  the  plenos  in  the 
nose  of  the  deteotor  body,  anti-lolng  of  the  strut  and  body  is  accomplished. 

The  induced  secondary  velocity  experienced  by  the  loe  collecting  probe  is  a function  of  airspeed  and 
1s  of  the  form  shown  in  Fig.  2.  At  low  airspeed  the  effect  of  the  ejeotor  predominates  the  Induced 
velocity  is  essentially  constant,  while  the  ran  effeot  predominates  at  high  airspeeds , and  the  induced 
velocity  approaches  proportionality  with  the  airspeed.  The  maximum  operating  airspeed  of  a helioopter 
normally  corresponds  to  a point  Just  above  the  knee  of  the  curve. 
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For  sensing  the  ioe  building  up  on  the  probe,  a beam  of  near  Infrared  (IB)  radiation  from  a minia- 
ture emitter  mounted  in  the  duct  wall  is  directed  acroee  the  leading  edge  of  the  probe  to  a seneor  in  the 
oppoeite  wall.  Ae  ice  ie  accreted  on  the  probe  the  IB  beam  ie  gradually  occluded  and  once  a given  thick- 
ness is  accumulated,  logic  circuits  within  the  unit  turn  on  electric  current  to  de-ice  the  probe.  The 
cycle  can  then  repeat.  The  probe  consists  of  a thin- wall  metal  tube  2.8  mm  in  diameter  that  has  suffici- 
ent eleotrical  resistance  to  serve  as  its  own  heater,  thus  no  delicate  internal  filament  is  required. 

An  ice  warning  signal  is  provided  at  the  initiation  of  de-icing  and  its  duration  may  be  preset  to  any 
desired  value. 

Figure  3 shows  the  maimer  in  which  the  voltage  output  of  the  infrared  sensor  varies  during  the 
ice/de-ice  cycle  of  the  detector.  As  ice  begins  to  form  on  the  probe  the  signal  output  level  begins  to 
rise  at  (a)  from  itc  quiescent  level.  When  the  output  voltage  reaohea  a reference  value  at  (b)  as 

a result  of  between  0.25  and  0.1*  mm  of  ice  having  accreted  on  the  probe,  the  probe  heater  and  signal 
lamp  are  activated.  The  probe  is  heated  for  a period,  t,  sufficient  to  ensure  de-icing,  and  after  a short 
cooling  period  ice  again  begins  to  form  on  it. 

MEASURHffiNT  OF  ICING  SEVERITY 

Of  particular  significance  is  the  linearity  of  the  voltage  ramp  between  points  d and  b in  Fig.  3, 
because  this  implies  that  the  slope  in  this  region  may  be  used  as  a measure  of  the  rate  of  icing,  or, 
more  significantly,  since  the  velocity  remains  essentially  constant  within  the  operating  airspeed  range, 
this  slope  oan  be  used  to  determine  the  liquid  water  content  of  the  icing  cloud.  Advantage  was  taken  of 
this  to  develop  two  types  of  icing  severity  indicators.  In  both,  the  time  required  for  the  sensor  volt- 
age to  rise  from  a level  corresponding  to  point  d to  the  level  TMf  corresponding  to  point  b was  measured 

and  used  to  determine  the  average  concentration  between  points  d and  b.  The  value  determined  was  dis- 

played until  updated  at  the  completion  of  the  next  icing  oycle,  or  until  a certain  "time-out"  period  had 
elapsed  after  which  the  display  reverted  to  zero  if  no  further  updates  had  been  received. 

The  first  type  of  severity  Indicator  was  a simple  3-level  device  in  which  one  of  the  three  lights, 
corresponding  to  Light,  Moderate  or  Heavy  icing  severity,  was  lit.  If  the  measured  rise  time  was  less 
than  a given  time,  t^,  then  the  Heavy  light  was  switched  on;  if  the  rise  time  was  greater  than  tj*  but 

less  than  a second  value  tg,  the  Moderate  light  was  lit;  and  if  the  rise  time  exceeded  tg,  then  the  indi- 

cation was  of  Light  icing.  The  reference  times  t^  and  to  could  be  adjusted  to  permit  the  desired  severity 
levels  to  be  set,  for  instance,  the  transition  from  Light  to  Moderate  might  be  set  at  0.25  g m“3  and  from 
Moderate  to  Heavy  at  0.5  g m_3. 

The  second  type  of  indicator  or  Icing  Severity  Meter  takes  the  rise  time,  t,  and  by  means  o*  a cal- 
culator chip  performs  the  operation: 

LWC  - k/t 

where  k is  a constant  that  characterizes  the  slope  of  the  voltage  ramp.  The  resultant  value  of  liquid 
water  content  is  displayed  digitally  by  a LED  display.  The  Icing  Severity  Mater  together  with  the  Ioe 
Detector  unit  are  shown  in  Fig.  i*.  As  with  the  3-Hght  display,  the  digital  LWC  indication  is  displayed 
until  updated  or  until  a time-out  period  is  exceeded. 

It  was  thought  that  the  3-level  system  might  be  more  acceptable  aa  a visual  icing  severity  indica- 
tion to  the  pilot  and  flight  crew  where  regular  military  or  commercial  operations  were  concerned,  while 
the  digital  Icing  Severity  Meter  was  considered  more  applicable  to  test,  qualification  and  research 
applications.  For  control  of  a rotor  blade  icing  protection  system  it  was  originally  thought  that  the 
3-level  system  would  suffioe,  but  the  digital  system  is  so  attractive  and  versatile  that  it  has  recently 
been  adapted  to  provide  this  control  function,  aa  will  be  described  later. 

PROBLEMS  ENCOUNTERED 

It  would  not  be  fair  to  give  the  impression  that  no  problems  were  encountered  in  the  development  of 
this  ice  detector,  for  indeed  there  were  problems. 

The  major  problem  that  had  to  be  overcome  was  that  due  to  contamination  of  the  optical  component s of 
the  sensing  system  by  dirt.  Contamination  was  present  in  the  ambient  air  and  in  the  compressor  bleed  air. 
Often  particles  in  the  air,  i.e.  smoke,  eto. , would  have  no  effect  in  dry  air  conditions,  but  would  be  a 
serious  source  of  contamination  when  water  droplets  were  also  present  in  the  air. 

Many  modifications  were  made  to  the  detector  to  try  to  solve  this  problem,  and  although  it  may  not 
be  possible  to  completely  eliminate  it,  the  rate  of  dirtying  has  been  reduced  to  what  is  considered  an 
acceptable  level  by  ventilating  the  optioal  surfaces  with  filtered  air.  The  venturi  suction  effect  in 
the  detector  duot  is  used  to  draw  this  air  over  the  windows. 

It  was  pointed  out  earlier  that  the  induced  velocity  was  essentially  constant  over  a range  of  air- 
speeds from  the  hover  up.  It  was  found  that  the  rate  of  icing  on  the  probe  is  not  constant  as  might  be 
expected,  but  that  at  low  airspeeds  when  the  entry  conditions  are  less  than  isokinetic,  droplet  concentra- 
tion enhancement  occurs  on  the  centre  line  of  the  duct  where  the  ice  is  sensed,  with  the  result  that  the 
lndloated  liquid  water  oontent  takes  on  the  sort  of  characteristics  shorn  in  Fig.  5.  If  the  induced 
velooity  is  too  great,  the  concentration  enhancement  is  exoeesive,  while  if  the  duot  velocity  is  too  low, 
the  lndloated  liquid  water  oontent  is  too  airspeed-dependent.  The  optimum  induced  velocity  is  that  which 
provides  sn  Indication  of  liquid  water  content  that  lias  within  a reasonable  tolerarce  range  of  the  true 
value  over  the  widest  possible  airspeed  range. 

Variability  in  the  bleed  air  supply  temperature  in  different  helicopter  installations,  and  the  pos- 
sibility of  Inadequate  heat  to  maintain  tha  detector  body  anti-iced  under  all  icing  conditions,  resulted 
in  a small  electrical  in-line  air  heater  being  installed  in  the  detector  strut.  The  power  to  the  heater 
is  modulated  by  a temperature  controller  using  a sensor  in  the  detector  body  to  maintain  the  body  at  about 
♦10  C. 


PRACTICAL  APPLICATIONS 


Prototype  loa  detectors  and  associated  liquid  water  content  systems  have  been  used  in  a number  of 
practical  applications. 

One  installation  was  made  on  a fixed  win#  aircraft  and  not  a helicopter.  Ibis  was  for  the  natural 
icing  flight  stage  of  the  qualification  of  the  DeHaviliand  Dash  7 STOL  aircraft.  In  thin  oasa  the  avail- 
able bleed  air  pressure  was  low,  only  18  psig,  but  because  the  flight  speeds  of  interest  wsre  in  the  rsnge 
150  knots  to  200  knots  the  low  air  pressure  was  not  a serious  drawback.  Before  installation  on  the  air- 
craft,  the  detector  was  calibrated  for  this  air  pressure  and  at  the  airspeeds  of  prime  interest,  i.e. 

150  and  180  knots,  and  correction  factors  were  determined  for  temperatures  near  zero  where  Ludlum  effects 
oocur. 

Pig.  6 shows  a portion  of  the  recorded  data  for  an  icing  flight  on  18  March  1978.  It  shows  the 
detector  output  voltage  and  the  air  temperature.  The  detector  voltage  is  seen  to  climb  as  the  probe 
ices,  to  drop  sharply  when  de-icing  occurs,  and  then  to  remain  low  for  about  a 6-second  probe  cooling 
period  before  starting  to  rise  again.  (The  grid  lines  are  spaced  at  20  seconds).  The  cloud  variability 
may  be  Judged  by  the  varying  slope  of  the  voltage  ramps.  The  mean  Bhope  of  each  ramp  between  the  appro- 
priate voltage  levels  was  (in  this  case)  measured  manually  and  the  liquid  water  content  computed  and 
written  on  the  reoord.  This  value  was  further  corrected  for  the  air  temperature  of  about  -3.5*0  and  the 
corrected  LWC  also  noted.  A vernier  hot  rod  was  also  installed  on  this  aircraft  and  the  thickness  of  ice 
on  it  is  noted  on  the  record  (i.e.  19  mm  at  time  13>10:11  and  20  mm  U minutes  later).  This  icing  rate 
on  the  hot  rod  Implies  a liquid  water  content  of  0.05U  g m"3,  whereas  the  loe  detector  records  give  an 
average  value  of  0.18  g m~3  over  the  same  period. 

A 3-level  icing  severity  indicator  was  installed,  and  its  Indication  is  shown  at  the  top  of  the 
recording.  Light  icing  was  indicated  until  the  un corrected  LWC  value  of  0.31  g m"3,  at  which  time  a 
Moderate  indication  was  given  until  updated  to  a Light  Indication  at  the  end  of  the  ensuing  icing  ramp. 

For  the  continuation  of  these  trials  in  the  following  winter  (1976-77)  «o  Icing  Severity  Meter  was 
installed  in  place  of  the  3-level  system,  so  that  a direct  read-out  of  liquid  water  content  was  obtained. 

A plot  of  the  relevant  data  for  a 1-hour,  26-minute  icing  enoounter  is  shown  in  Fig.  7.  Air  temper- 
ature was  about  -6*C  at  an  altitude  of  5,200  feet.  Airoraft  speed  was  about  l60  knots,  lbs  liquid  water 
content  as  given  by  the  Ice  Dateotor/loing  Severity  Meter  system  is  seen  to  average  about  0.2  or  0.25 
g m'3  for  the  enoounter  with  peak  values  between  0.U  and  0.5  g m“3.  A olimb  to  about  6000  feet  to  look 
for  more  severe  conditions  resulted  instead  in  leas  severe  conditions,  so  the  airoraft  returned  to  its 
original  altitude.  The  less  severe  conditions  at  increased  altitude  are  also  assn  during  the  climb  out 
at  the  end  of  the  enoounter. 

The  next  and  final  example  is  that  of  the  loe  deteotor  installation  on  a UH-1H  helicopter  equipped 
with  an  experimental  rotor  blade  de-icing  system.  In  this  installation  the  loe  deteotor  output  ramp 
voltage  is  fad  into  an  Intagrating  Rate  Unit  (IRU).  This  unit  computes  the  liquid  water  content  from  the 
ramp  rise  time  and  displays  it  digitally  as  in  the  Iolng  Severity  Meter.  In  addition  it  oomputea  from 
this  and  the  detector  oyole  time  the  equivalent  ioe  accretion  lnorement  scaled  to  represent  that  at  some 
location  on  the  rotor  blades  (usually  50%  radius).  These  increments  are  acoumulatad  until  a pre-set 
thickness  is  reached  or  exoeeded,  at  which  point  the  blade  de-lolng  oyole  is  Initiated  and  the  accumula- 
tor ia  reset  to  zero. 

Figs.  8a  and  8b,  traced  from  the  flight  recorder  film,  show  two  portions  of  the  same  blade  ioing 
cycle.  Fig.  8a  shows  the  5th,  6th  and  7th  loe  deteotor  oyoles  after  lnltlationi  the  liquid  water  content 
value  displayed  has  been  entered  beside  each.  The  trace  below  the  detector  output  trace  represents  the 
incremented  blade  ice  thiokneas  as  oomputed  by  the  Integrating  Rets  Unit)  the  thickness  value  stored  in 
the  IRU  has  been  added  beside  each  lnoraxent.  Fig.  8b  shows  the  U3rd,  U*th  and  U5th  ice  detector  cycles. 
After  the  UUth  cyole  the  blade  loe  Integrator  has  incremented  to  0.1*03  incbi  the  pre-set  ice  thickness 
for  blade  de-lolng  was  0.U  15  inch,  so  on  receipt  of  the  next  lnorement  (0.017  inch)  the  integrator  incre- 
ments momentarily  to  0.U20  inch  and  then  resets  to  zero,  at  the  same  time  providing  an  initiating  pulse 
to  the  rotor  blade  de-icing  controller. 

It  may  be  wondered  why  one  does  not  simply  sum  the  number  of  loe  deteotor  oyoles  and  initiate  de-lolng 
after,  say,  US  oyolss.  Figs.  8a  and  8b  demonstrate  why  this  procedure  is  not  always  dseirable.  It  may 
be  noted  that  loe  detector  oyole  number  6 (Fig.  8aJ  gave  an  lnorement  of  0.012  lnah  for  a liquid  water 
content  of  0.3U  g m3,  while  oyole  number  UU  (Fig.  ob)  gave  an  increment  of  0.015  inch  for  a LWC  of  0.33 
g m3,  the  reason  being  that  although  the  LWC  was  almost  the  same,  the  oyols  tins  had  increased  by  $ sec- 
onds (a  28%  increase)  beoauae  of  the  longer  probe  oooling  time  at  -2.5*0  than  at  -U*C. 

CONCLUDING  REMARKS 

nils  paper  has  attempted  to  demonstrate  the  latest  advanoes  in  helicopter  loe  detection  and  liquid 
water  content  measurements.  It  may  be  expected  that,  with  the  advent  of  micro-processors,  grsatsr 
advanoes  will  be  seen  in  the  not-too-dl stant  future,  particularly  in  the  area  of  processing  the  informa- 
tion provided  by  the  primary  sensors  to  correot  the  final  outputs  for  suoh  variables  as  amblsnt  tempera- 
ture, airspeed,  bleed  air  pressure  end  so  forth. 

Where  mixed  ioing  conditions  an  concerned,  a vast  amount  of  development  work  will  be  neosasary  to 
avoirs  a relatively  simple  system  to  measure  suoh  conditions  on  a routine  basis. 
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FIG.  1:  SIMPLIFIED  OUTLINE  OF  DYNAMIC  ICE  DETECTOR 


AIRSPEED 

FIG.  2:  INDUCED  VELOCITY  THROUGH  THE  ICE  DETECTOR 
AS  A FUNCTION  OF  AIRSPEED 
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FIG.  6:  LIQUID  WATER  CONTENT  INDICATION  AS  A FUNCTION  OF  AIRSPEED 
SHOWING  EFFECT  OF  INDUCED  VELOCITY 
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FIG.  6:  PORTION  OF  RECORDED  DATA  FROM  DHC-7  ICING  QUALIFICATION  FLIGHT 
ON  18  MARCH  1976,  SHOWING  VOLTAGE  OUTPUT  FROM  ICE  DETECTOR  AND 
DERIVED  LIQUID  WATER  CONTENTS 


FIG.  8:  TWO  PORTIONS  OF  RECORDED  DATA  FROM  UH-1H  NATURAL  ICING  TEST  FLIGHT  ON  31  MARCH  1977, 
SHOWING  ICE  DETECTOR  OUTPUT  RAMP  AND  ACCUMULATED  VALUE  IN  INTEGRATING  RATE  UNIT 
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SUMMARY 

One  of  the  large  test  chambers  of  the  Altitude  Test  Facility  at  the  NGTE  has  been  successfully  used 
for  full-scale  trials  on  the  front  fuselage  and  engine  intakes  of  the  Westland  Sea  King  and  the  Westland/ 
Aerospatiale  Lynx  helicopters  under  simulated  wet  icing  conditions. 

The  paper  describes  tests  made  to  examine  the  effects  of  air  temperature,  air  speed,  water  concen- 
tration and  aircraft  attitude  on  the  extent  of  ice  accretion.  The  possibility  of  ice  ingestion  by  the 
engines  during  melt-off  has  also  been  studied  using  high-speed  videotape  recordings. 

Data  are  also  presented  from  non-rotating  tests  on  full-scale  sections  of  main  rotor  blade  made  to 
examine  the  effects  on  accretion  of  blade  speed,  pitch  angle,  air  temperature  and  water  concentration. 

The  effects  of  cyclic  pitch  change  have  also  been  examined  and  typical  results  are  presented. 

The  paper  includes  general  comments  on  the  use  of  an  altitude  cell  for  helicopter  icing  investiga- 
tions and  the  accumulation  of  data  leading  to  clearance  for  flight  in  icing  conditions. 


1.  INTRODUCTION 

Clearing  a helicopter  for  flight  in  icing  conditions  requires  the  assessment  of  data  from  a variety 
• of  sources  ranging  from  small  test  rigs  to  full-scale  flight  trials.  The  final  judgement  on  the  limits  of 

safe  operation  naturally  leans  heavily  on  the  flight  trial  results,  but  a considerable  degree  of  confidence 
can  be  gained  and  much  time  and  expense  saved  by  carefully  designed  experiments  using  ground-based  test 
rigs  and  facilities.  These  enable  icing  to  be  studied  under  closely  controlled  conditions  and  the  influ- 
ence of  each  major  variable  examined  separately.  Ground  facilities  can  also  be  used  for  development  test- 
ing, for  example  the  optimisation  of  the  power  input  distribution  to  thermal  anti-ice  or  de-ice  systems, 
and  results  obtained  more  quickly  and  precisely  than  is  possible  in  flight. 

This  p. per  describes  the  work  done  by  the  National  Gas  Turbine  Establishment  (NGTE)  at  Pyestock  to 
study  at  full-  scale  ice  build  up  on  the  forward  fuselage  and  engine  air  intakes  of  Sea  King  and  Lynx 
helicopters  ani  the  possibility  that  ice  shed  during  melt-off  may  be  ingested  by  the  engines.  Ways  in 
which  this  problem  was  overcome  are  described.  Tests  to  examine  ice  accretion  on  rotor  blade  specimens 
are  also  covered. 

2.  TEST  FACILITIES 

NGTE  has  undertaken  icing  trials  on  aero-engines  for  the  past  15  years.  With  the  remarkable  growth 
of  engine  power  that  has  taken  place  during  this  time  and  the  corresponding  increase  in  the  size  of  the 
supporting  test  facilities,  in  particular  the  increase  of  compressor/exhauster  plant  capacity  and  cold  air 
supply,  the  Establishment  has  acquired  a unique  capability  for  investigating  helicopter  icing  problems. 

Figure  1 shows  the  largest  of  the  altitude  test  chambers  at  NGTE,  Cell  3 West.  The  test  chamber, 
which  has  a diameter  of  25  ft,  can  be  seen  on  the  left  of  the  picture;  the  air  cooler,  refrigeration  plant 
and  cold  store  are  on  the  right.  A detailed  description  of  this  Facility  is  given  in  Reference  1. 

The  Cell  can  easily  accosraodate  the  fuselage  of  a Sea  King  or  Lynx  helicopter  and  by  correctly 
positioning  it  in  a free-jet  air  stream  a large  part  of  the  forward-facing  fuselage  surfaces  including  the 
windscreen,  cabin  roof  and  engine  intakes  can  be  subjected  to  icing  conditions.  Figure  2 shows  the  test 
arrangement.  Air  drawn  from  atmosphere  passes  through  the  cooler,  the  inlet  pressure  control  valve  and 
into  the  test  chamber  through  the  blowing  nozzle.  This  produces  a jet  of  air  some  8 ft  in  diameter  which 
simulates  the  flight  velocity.  The  icing  cloud  is  generated  by  injecting  water  from  an  array  of  242  air- 
assisted  atomisers  which  can  be  seen  in  Figure  3. 

No  attempt  is  made  to  simulate  the  effects  of  rotor  downwash.  Tests  made  using  a one-twentieth 
scale  model  of  the  cell  installation  showed  that  tha  pressure  distribution  over  the  fuselage  upper  surface 
was  virtually  identical  to  that  measured  in  a conventional  wind-tunnel.  It  was  therefore  concluded  that 
the  free-jet  teat  arrangement  gives  a close  representation  of  conditions  encountered  in  forward  flight. 

Tha  test  conditions  most  commonly  used  represent  flight  at  speeds  from  90  to  160  kt  with  sir  temp- 
eratures in  the  range  -2  to  -24°C.  The  liquid  water  content  is  set  to  meet  the  requirements  of  AvP  970 
(Reference  2)  for  either  maximum  continuous  or  periodic  maximum  conditions. 

To  avoid  possible  damage  by  ingestad  ica  most  teats  are  run  without  engines  fitted,  air  being  drawn 
into  the  intakes  by  an  ajactor  system  operated  by  tha  plant  pressure  air  supply.  Correct  thermal 
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I: 

1 

• imulat ion  of  the  heat  inputs  from  the  engine  front  frame  and  the  gearboxes  is  obtained  by  electric  1 

heating.  Some  tests  are  run  with  an  engine  and  for  these  the  power  turbine  is  locked  to  obviate  the  need 

for  a dynamometer.  t 

: 

After  completing  an  icing  test  the  water  sprays  are  turned  off  and  the  air  temperature  slowly  < 

increased  to  above  freezing  point  to  establish  whether  ice  ingestion  occurs  during  melt-off.  To  do  this  • 

it  is  necessary  to  examine  high  speed  videotape  records  frame  by  frame  since  the  lumps  of  ice  move  too 
fast  to  be  followed  by  direct  viewing  of  the  TV  monitors. 

3.  INTAKE  ICING 

Figure  4 shows  typical  results  of  icing  tests  after  16  and  30  minutes  at  temperatures  of  -3  and 

-10  C on  a Sea  King  fitted  with  pitot  intakes.  The  forward  speed  was  120  knots  and  the  liquid  water 

content  was  as  specified  in  AvP  970,  see  Reference  2 for  continuous  maximum  conditions.  The  engine  and 
intake  flare  anti-icing  systems  prevented  dangerous  ice  accretion  on  the  intakes  but  with  this  type  of 
intake  the  engines  are  vulnerable  to  the  ingestion  of  ice  shed  from  the  front  fuselage. 

In  an  attempt  to  overcome  this  problem  the  so  called  "mushroom"  intakes  were  fitted  and  Figure  3 

shows  these  after  icing  tests  at  four  temperatures.  The  front  face  of  the  port  "mushroom"  was  unheated 

while  that  of  the  starboard  had  1.2  kU  of  electrical  heating.  It  can  be  seen  that  ice  formed  on  the 
unheated  "mushroom"  at  -A°C  and  on  the  heated  "mushroom"  at  -15°C.  Although  in  forward  flight  these 
intakes  prevented  ingestion  of  ice  from  the  front  fuselage,  sometimes  ice  from  the  "mushroom"  front  face 
was  ingested.  Large  ice  deposits  formed  on  the  unheated  area  behind  the  intake  cowls  caused  by  water  run- 

back  from  the  heater  mats  and  some  of  this  ice  was  also  ingested.  j 

Tests  with  a Sikorsky  Foreign  Object  Deflector  (FOD)  fitted  to  the  cabin  roof  in  front  of  the  pitot 
intakes.  Figure  6,  indicated  that  with  light  ice  deposits  no  ice  ingestion  occurred.  However  considerable 
ice  accretion  particularly  around  the  edges  occurred  after  30  minutes  icing,  Figure  7A,  and  during  melt- 
off  some  of  this  ice  was  ingested  through  the  cut  aways  at  the  top  rear  edges  of  the  deflector.  With  a 

de-icing  fluid  exuded  from  porous  strips  inserted  flush  in  the  front  face  of  the  FOD  a marked  reduction 

in  ice  accretion  was  achieved,  Figure  7C.  Some  of  the  fluid  was  carried  over  in  the  air  stream  into  and  ' 

around  the  engine  intakes  reducing  ice  accretion  in  these  areas  and  obviating  the  need  for  the  auxiliary  " 

heater  mats. 

Sideways  facing  intakes.  Figure  8,  with  electrically  heated  mats  for  anti-icing  have  given 
encouraging  results  at  temperatures  below  -15°C.  The  mats  were  divided  into  a number  of  separate  sections 

each  with  its  own  heating  intensity  control  thus  enabling  the  optimum  heat  distribution  to  be  determined.  \ 

Figure  9 shows  the  port  intake  after  30  minutes  icing  at  -16°C.  Run  back  ice  built  forward  from  the 
unheated  area  behind  the  intake  but  this  shed  clear  during  melt-off  and  did  not  constitute  an  ingestion 

problem.  Ice  shedding  from  the  front  of  the  fuselage  also  passed  clear  of  the  intakes. 

Icing  tests  were  also  made  on  the  smaller  Lynx  helicopter  fuselage  which  has  the  engines  and 
intakes  mounted  on  top  of  the  fuselage  behind  the  rotor  gearbox,  Figure  10.  Ice  accretion  along  the  top 
of  the  rotor  gearbox  cowling  on  the  original  design,  Figure  11A,  resulted  in  a number  of  instances  of 
engine  ice  ingestion  and  one  example  can  be  seen  in  the  high  speed  videotape  recording  reproduced  in 
Figure  12.  Streamlining  the  gearbox  cowling  and  smoothing  the  windscreen  "eyebrows"  reduced  very  consider- 
ably the  chance  of  ice  ingestion  from  these  areas. 

Auxiliary  air  intakes,  such  as  those  required  for  oil  cooling,  are  at  risk  in  icing  conditions  as 
well  as  the  main  engine  intakes.  Figure  13  shows  the  extent  to  which  an  auxiliary  intake  can  become 
blocked  with  ice. 

4.  ROTOR  BLADE  ICING  ( > 

Cell  3 West  is  not  large  enough  to  accommodate  a complete  main  rotor,  but  full-scale  tests  on 
complete  tail  rotors  are  feasible.  No  such  tests  have  yet  been  made  but  some  are  planned  to  examine  the 
effects  of  inertia  loading  on  ice  accretion.  It  is  thought  that  the  rig  could  be  used  to  clear  specific 
designs  for  operation  in  icing  conditions. 

NGTE  has  made  non-rotating  tests  on  sections  of  full-scale  main  rotor  blade  to  determine  the  effects 
on  accretion  of  water  concentration,  droplat  siza,  air  temperature  and  velocity;  the  influence  of  pitch  t 

angle  and  aerofoil  shape  have  also  bean  investigated.  The  test  rig  is  shown  in  Figure  14.  The  blade 
specimen  is  mounted  in  a gymbil-type  suspension  so  that  the  pitch  and  yaw  angles  can  be  varied  cyclically. 

The  airspeed  can  be  adjusted  to  represent  conditions  at  selected  spanwise  stations. 

The  photograph  shown  at  the  top  of  Figure  IS  shows  the  extensive  ice  accretion  observed  after  a 
7}  minute  test  at  -8°C  with  a water  content  of  0.65  g/m* . The  picture  is  of  a chordwise  section  through 
the  ice. 

Natural  icing  conditions  are  more  complex  than  those  represented  by  an  air  stream  containing  only 
supercooled  water  droplets  and  NGTE  has  therefore  examined  ice  deposition  in  a mixed  environment,  that  is 
one  containing  solid  ice  particles  and  suparcoolad  water  dropleta.  Figure  IS  compares  the  build  up  after 
two  10  minute  testa,  the  first  at  the  mixed  condition  and  the  second  with  the  same  water  content  but  with 
no  ice  particles.  A comparison  of  results  from  tha  two  tests  indicated  that  the  presence  of  the  solid  ice 
reduces  the  thickness  of  the  accretion  and  produces  a deposit  having  a much  smoother  contour. 

5.  CONCLUDING  REMARKS 

This  paper  has  presented  a brief  outline  of  work  done  at  NGTE  to  investigate  several  aspects  of 
the  helicopter  icing  problem.  Tha  test  facilities  available  at  the  Establishment  allow  experiments  to  be 
made  at  full  scale.  This  has  proved  to  be  far  more  affective  for  development  testing  than  either  flight 
tests  or  open  air  facilities  which  rely  on  prevailing  weather  conditions  to  provide  the  icing  environment. 


However  whilst  a test  facility  can  provide  a considerable  degree  of  confidence  in  the  ability  of  a 
helicopter  to  operate  in  icing,  the  final  judgement  on  the  limits  of  safe  operation  must  be  based  on 
flight  trials. 
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Fig.  1 Cell  3 west 
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Fig.2  Arrangement  of  Sea  King  helicopter  fuselage  in  cell  3 west 
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Fig.4  Sea  King  with  pitot  intakes 
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30  mins  AT  0.7g/m3  25  mins  AT  0 45g/m3 

-5  AND  120kt  -15°C  AND  120kt 


Fig.5  Sea  King  with  mushroom  intakes 


3 mins  AT  1.5g/m3  3 mins  AT  1 2g/m3 

-4°C  AND  90 kt  -10°C  AND  90kt 

UNTREATED  FOD  FRONT  FACE 


27  mins  AT  0 6g/m3 
3 mins  AT  1.2g/m3 
-10°C  AND  90kt 


FOD  FRONT  FACE  TREATED  WITH  DE  ICING  FLUID 

Fjg.7  Sea  King  with  Sikorsky  foreign 
object  deflector  (FOD) 


27  mins  AT  0 4g/m3 
3 mins  AT  0 8g/m3 
-16°C  AND  90kt 


Fig.9  Sea  King  with  sideways  facing  intakes 
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ORIGINAL  COWLING  30  mins  AT  0 7g/m3 

-5°C  AND  120kt 


MODIFIED  COWLING  30  mins  AT  0 7g/m3 

-5°C  AND  120kt 


Fig.  1 1 Lynx  - pre  and  post  modification 
to  the  rotor  gearbox  cowling 
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Fig.  12  Lynx  - high  speed  television  sequences  ice  ingestion 


30  mins  AT  0.7g/m3 
-5°C  AND  120kt 


Fig.  1 3 Lynx  - ice  accretion  on  oil  cooler  air  intake 


SECTION  OF  ICE  ON  JETRANGER  ROTOR  BLADE 
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Fig.  1 5 Ice  deposition  on  rotor  blade  test  sections 
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AFTER  Th  mins  AT  0.65g/m3  WATER 


8°C  AND  339  ft/s 


ICE  ACCHETION  AND  ITS  EFFECTS  ON  AERODYNAMICS  OF  UNPROTECTED 
AIRCRAFT  COMPONENTS 


4-1 


I 

\ 

| 

by 

Boris  Laschka  and  Rudolf  E.  Jesse  t 

Messerschmi tt -Bolkow-Blohm  GmbH 
Ottobrunn  / Munich,  Germany 

SUMMARY 

This  paper  gives  a survey  on  theoretical  and  experimental  methods  applied  to  determine 
the  ice  accretion  and  its  effects  on  the  respective  aerodynamics  on  unprotected  compo- 
nents of  a modern  jet  transport  aircraft. 

Together  with  an  introductory  outline  of  the  mechanism  and  the  basic  mathematical  equa- 
tions and  correlations  ruling  ice  generation,  the  experience  obtained  for  the  European 
Airbus  A3 00  is  used  to  describe  these  methods • 

In  addition, an  example  for  the  prediction  of  trajectories  of  shedding  ice  particles  as 
applied  for  the  German  jet  transport  aircraft  VFW  6l**  is  given  herein. 


INTRODUCTION 

Early  in  the  development  phase  of  the  A300,  it  was  decided  not  to  protect  the  tail 
surfaces  and  the  air  intake  system  of  the  Auxiliary  Power  Unit  (APU)  by  a de-icing  or 
anti -icing  system  in  order  to  save  weight  at  the  rear  of  the  aircraft,  electric  and 
hydraulic  power,  and  to  simplify  system  design  (see  fig.  l).  * In  order  to  ensure  air- 

worthiness under  all  possible  extreme  icing  conditions,  thorough  investigation,  to 
guarantee  full  flight  safety,  sufficient  flying  qualities  and  controllability  and  APU 
system  function  have  been  performed.  These  studies  covered  the  theoretical  and  experi- 
mental det erminat ion  of  the  amount  of  ice,  ice  build  up,  and  shape  and  the  resulting 
aerodynamic  implications,  completed  by  flight  tests  under  natural  icing  conditions 

Ct].[a]. 


1.  UNPROTECTED  TAIL  OF  THE  A300 

Icing  mechanism:  looking  into  the  flow  field  around  an  airfoil,  as  shown  in  fig.  2, 
one  may  find,  that  the  supercooled  water  droplets  due  to  their  inertial  forces  deviate 
from  the  curvature  of  the  streamlines  in  the  vicinity  of  the  leading  edge.  The  larger 
droplets  tend  to  deviate  more  from  the  streamlines  having  less  curved  trajectories  and 
then  impinging  on  the  airfoil.  So,  the  water  caption  problem  consists  in  determining 
the  droplet  trajectories  and,  more  particulary,  the  tangent,  or  limit  trajectories  as  a 
function  of  droplet  size  (d). 


Basic  equations  ruling  droplet  trajectories:  the  governing  differential  equation 
may  easily  be  derived  by  the  Newton  law 

£ = (t) 

d t 

with  m as  mass  of  the  droplet 


O 4 11  (rf)5 

*»water  3 2 


(2) 


and  the  only  force  acting  on  the-  droplet 
given  by  its  drag  ft 


CD  71  ' > 
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£ + 5 = 0 


(3) 

(4) 


equation  (l)  can  be  transferred  to 


dV  3 Sair  (0-C)  lO-frl 

dt  “ * 2 water  °D'  2 


Vector  diagramm  for  droplet  and 
(5)  air  velocity 


Since  the  droplet  Reynolds  Numbers  are  very  small, the  STOKE's  law  ( Cp  = 
can  be  used  with  good  approximation.  The  droplet  Reynolds  Number  itself 
from  the  relative  velocity  iC-Ol  of  the  droplet  with  respect  to  the  air 


Re 


droplet 


-cr1-  IV  - 01 

air 


24/Re)  for  drag 
is  determined 

(6) 


Generally,  the  differential  equation  is  solved  by  classical  difference  methods  in  combi- 
nation with  a panel  flow  field  calculation. 


* Note : 

The  decision  not  to  protect  the  wing  slat  area  between  the  fuselage  and  the  engine 
pylon  was  not  made  at  this  time,  this  was  decided  later  according  to  experimental 
results. 


i 
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Might  teg i me  of  maximum  icing:  for  most  transport  aircraft,  the  severest  icing  within 
the  flight  envelope  will  be  encountered  during  holding  where  usually  the  airplane  by 
traffic  control  reasons  is  not  allowed  to  change  altitude  or  location  arbitrarily  to 
avoid  icing  regimes  and  where  appreciable  flight  times  have  to  be  expected.  This  applies 
also  to  the  A")00  for  which  the  severe  special  requirement  of  45  minutes  flight  under 
maximum  continuous  icing  conditions  had  to  be  fulfilled.  Using  the  statistical  metereo- 
logical  data  as  given  in  FAR  25  and  taking  into  account,  that  the  average  water  captions 
pe«  time  unit  is  proportional  to  liquid  water  content  LWC,  the  true  flight  speed  V and 
a caption  factor  K (depending  on  droplet  diameter  and  tail  flow) 


water 

dt 


= K • V • LWC 


and  that  the  ambient  temperature  in  the  stagnation  region  of  the  tail  surfaces  has  to 
Se  inn  eased  by  AT  = 0,2  T„,  • M2  • cos2  9 ( 'P  being  tail  l.e.  sweep),  one  may  conclude 
that,  for  the  ATOO  tail, biggest  ice  accretion  occurs  for  ambient  temperature  from 
5t  to  - 106C  and  altitudes  around  14  000  ft,  as  shown  in  fig.  3 for  a droplet  size  of 
d 20  p (demonstration  case!).  Though  cruise  velocities  are  higher  than  at  holding,  no 
severe  icing  is  to  be  expected,  as  because  of  the  dynamic  temperature  increase  only  at 
very  low  temperatures  at  which  LWC  is  small  icing  occurs. 


Influence  of  droplet  size  on  water  caption:  for  different  airfoil  geometries, one 
will  obtain  different  droplet  sizes  for  which  maximum  ice  accretion  occurs.  Conse- 
quently, the  Influence  of  droplet  size  on  total  water  caption  of  the  Airbus  tail  sur- 
faces was  established  by  impingement  analysis.  Analysis  [3]  ,[4]  shows,  see  fig.  4,  that 
the  biggest  water  caption  will  be  with  droplets  diameter  between  20  p and  28  p both  for 
the  vertical  as  well  as  for  the  horizontal  tail.  Furthermore,  the  calculations  give  evi- 
dence that  ice  accretion  will  be  more  severe  for  the  horizontal  tail.  This  is  due  to 
the  lact  that  the  horizontal  tail  has  an  absolute  smaller  leading  edge  radius  than  the 
vertical  tail.  It  is  Interesting  to  note  that  for  the  Boeing  747  also  the  horizontal 
tail  showed  to  be  the  surface  with  the  biggest  water  caption,  but  that  the  critical 
droplet  size  around  28  p was  somewhat  bigger  than  at  the  Airbus. 


Types  of  ice  shapes:  of  utmost  importance  for  the  alternation  of  the  aerodynamic 
qualities  is  not  only  the  amount  of  ice  accretion  but  rather  its  shape,  see  fig.  5- 

Chordwise  ice  shapes 

lienerally,  two  types  of  ice  shapes  may  be  observed  on  airfoils. 

• Rime  ice  results  from  low  surface  temperatures  and  low  LWC.  The  water  droplets  freeze 
immediately  after  impinging  on  the  leading  edge  of  the  airfoil  and  form  a soft  mil- 
kywhite  smooth  and  streamline  ice  shage  around  the  leading  edge.  It  does  not  occur 
when  surface  temperatures  surpass  -12DC.  From  the  aerodynamic  point  of  view, these 
forms  are  not  too  bad,  as  they  do  not  adversely  effect  lift-  and  drag  characteris- 
tics of  the  airfoil  and,  sometimes,  even  result  in  an  increase  of  lift  (enlargement 
of  a ir foi 1 -chord ) . 

. Glaze  ice  (mushroom  ice,  double  horn)  forms  at  temperatures  just  below  freezing 
point  and  high  LWC.  Glaze  ice  is  normally  formed  at  temperatures  between  0 C and 
-8°C.  At  these  temperatures,  water  droplets  do  not  freeze  immediately  after  impinge- 
ment but  run  first  on  the  profile  upward  or  downward  from  the  stagnation  point  q 
before  freezing.  In  this  case, the  temperature  at  the  stagnation  point  is  about  0 C 
or  above,  i.e.  the  so-called  Ludlam  limit  is  exceeded  (freezing  factor  < l).  The  re- 
sult is  a deposit  of  ice  with  two  horns  like  protrusions  and  normally  is  transparent 
ind  harder  than  rime  ice.  Such  an  ice  shape  is  aerodynamically  very  unfavourable 
since  the  horn  on  the  suction  side  spoils  the  flow  and.  reduces  the  lift  whereas  the 
drag  is  increased.  Naturally  that  mainly  glaze  ice  need  to  be  considered  if  aerody- 
namically critical  icing  is  to  be  investigated. 

Spanwise  ice  shapes  for  swept  wings 

On  swept  wings, in  addition  to  the  Just  discussed  chordwise  ice  shapes  accretions  with 
considerable  segmented  "notchings"  similar  to  a "lobster  tail"  occur.  The  segments  are 
curved  in  such  a way,  that  their  ends  are  bowed  towards  the  wing  tip. 


Kmperical  method  for  determination  of  ice  height:  computations  assuming  that 
• all  droplets  impinging  the  contour  freeze  immediately 
. breaking-off  of  ice  deposits  is  neglected 

do  not  show  the  mushroom  shape,  but  a shape  similar  to  rime  ice.  Such  a calculation 
always  predict  ice  heights  considerable  in  excess  of  the  real  ones.  This  is  a real  gap 
In  knowle  diej  Nevertheless,  empirical  methods  are  available,  from  which  mushroom  ice 
shapes  may  be  derived.  This  method,  see  fig.  6,  makes  extensively  use  of  measured  upper 
horn  lengths  h,,s , lower  horn  length  hLg  and  thickness  at  the  stagnation  point  hgT  for 
different  yet  similar  profiles,  which  were  plotted  versus  a quantity  I with  EM  Is  the 
total  droplet  collection  efficiency  and  h*/c  is  the  ratio  of  profile  thickness  [3],  [5]- 
A thus  determined  ice  height  at  the  horizontal  tail  of  the  A300  is  plotted  versus  span 
in  the  figure  for  a selected  case  and  then  compared  with  B 747*  In  principle,  ice 
height  decreases  with  increasing  profile  thickness  or  - what  is  tantamount  - for  con- 
stant thickness  ratio  with  increasing  profile  chord.  Almost  the  same  ice  height  is 
found  at  the  tip  for  both  horizontal  tails  yet  a different  decrease  of  height  towards  the 
fuselage  due  to  differences  in  airfoil  thickness. 


I 
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Experimental  determination  of  ice  shapes:  it  is  obvious  that  the  highly  three-dimen- 
sional structure  of  the  ice  deposltes  as  given  by  the  mushroom  and  lobster-tail  shape 
is  not  theoretically  predictable  with  sufficient  reliability  whereas  prediction  of  the 
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total  captured  ice  can  well  be  done,  as  shown  before.  Therefore, for  the  Airbus, the  cal- 
culation procedure  was  complemented  by  experiments  in  the  icing  tunnel  of  NASA,  Lewis 
Center,  Cleveland,  Ohio.  Two  horizontal  tail-models  of  the  A300  (model  1 with  constant 
profile  thickness,  model  2 was  a portion  of  the  horizontal  tail  in  the  original  design 
configuration)  were  investigated  [6].  Ice  accretions  have  been  determined  by  making 

plaster  casts  and  photographs.  In  order  to  gain  an  insight  into  the  significance  of  the  > 

icing  parameters  ambient  temperature  Too,  flow  speed  V,  droplet  diameter  d,  liquid  water 
content  LWC , icing  duration  t and  airfoil  incidence  , have  been  varied  during  the  tests  i 

[6],  [7]«  Results  show  (see  fig.  7)s  j 

(l) When  varying  LWC  and  time,  while  droplet  diameter,  velocity  and  temperature  are 
kept  constant,  one  obtains  the  same  ice  shapes  for  all  LWC • t = const. 

© If  time  varied  with  the  same  conditions,  ice  height  will  be  approximately  propor-  < 

tinal  to  the  time  period  of  icing,  while  impingement  limits  are  independent  of  time. 

A slight  nonlinear  increase  with  time  could  be  observed  for  already  big  ice  accre- 
tions, which  may  be  contributed  to  the  increase  of  capture  area  with  time. 

(j)  1 f velocity  and  time  are  simultaneously  varied  at  constant  LWC, almost  identical 
ice  heights,  yet  different  impingement  limits  and  different  horn  configurations 
will  be  obtained  on  the  condition  V't  = const. 

©If  angle  of  attack  is  varied,  different  ice  shapes  will  be  obtained;  this  condition 
implies  migration  of  the  stagnation  point. 

From  steps  1 through  3 it  follows: 

ice  height  - LWC ■ t • V = const, 
ice  shape  - LWC • t = const. 

Conversion  of  icing  test  results  to  real  flight  conditions:  conformity  with  actual 
flight  conditions  could  not  be  achieved  completely  during  icing  tests  for  some  para- 
meters like  LWC  and  velocity.  Moreover,  since  three-dimensional  flow  effects  could  not 
be  exactly  simulated,  results  had  to  be  converted  to  actual  flight  conditions.  There- 
fore, a special  procedure  [7], [8]  has  been  developed  for  the  Airbus  to  deduce  the  ice 
accretion  making  use  of  both  calculation  and  test  results. 

• First  maximum  ice  accretion  has  been  calculated  and  substantiated  by  experiment.  The  / 

result  was  that  the  highest  ice  deposit  occur  near  the  tip  of  horizontal  and  vertical 
tail,  the  higher  one  on  the  horizontal  tail.  Its  overall  height  was  slightly  below 
3 inches.  Therefore,  this  value  has  been  adapted  as  a maximum.  This  is  in  accordance 
with  observation  of  Boeing  Company  [5]  that  within  26  million  flight  hours  on  diffe- 
rent airplanes  no  bigger  accretions  have  occured.  • 

. Spanwise  variation  of  the  amount  of  ice  has  been  assumed  to  be  according  to  impinge- 
ment analysis  (see  fig.  8).  This  has  been  substantiated  for  several  tail  surface 
sections  by  experimental  results. 

. Chordwise  and  spanwise  ice  shapes  have  been  fully  taken  according  to  the  experimen-  t 

tal  shapes.  Slight  adjustments  from  test  to  flight  conditions  have  been  made  using 
similarity  laws  [9]« 

With  this  procedure,  the  most  objectionable  ice  shapes  for  the  A300  tail  have  been  1 

evaluated  and  selected  in  order  to  be  further  investigated  in  wind  tunnel  and  flight 
test.  Fig.  9 depicts  these  ice  shapes  over  the  span  for  the  horizontal  tail.  The  cor- 
responding notches  are  the  same  as  shown  before. 

Aerodynamic  investigations:  extensive  wind  tunnel  testing  with  artificial  ice  accre- 
tions of  the  just  described  shape  was  performed  in  order  to  establish  its  influence  on 
the  most  important  aerodynamic  data.  These  tests  were  done  with  complete  models  as  well 
as  pure  tail  models.  The  investigations  Included  the  determination  of  the  influence  of  < 

ice  shapes,  ice  size,  and  notchings. 

Results^  In  fig.  10a,  the  horizontal  tail  lift  (c^)  versus  angle  of  attack  as  deter-  , * 

mined  with  a tail  model  in  the  low  speed  wind  tunnel  of  Braunschweig  [10J  shows  that 

the  deterioration  effect  of  ice  accretion  sensitively  depends  on  elevator  angle.  For 

high  nose  down  moment  elevator  angles,  the  loss  in  tail  power  at  high  negative  angles 

of  attack  - what  is  equivalent  to  high  negative  settings  - has  shown  to  be  not  too 

serious.  Maneouvers  of  interest  in  that  range,  such  as  stall  recovery  or  go  around 

where  the  pilot  wants  or  has  to  Increase  the  nose  down  moment  remain  nearly  uneffected. 

On  the  other  hand, for  high  nose  up  moment  elevator  deflections,  a loss  in  maximum  nega- 
tive tall  lift  can  be  observed  which  reduces  the  trimming  and  pulling  capability  at  > 

low  speeds.  Lift  slope  was  associated  only  by  minor  changes.  This  was  different  to  re- 
sults on  the  vertical  tail  where  an  increase  in  lift  slope  according  to  the  increase  of 
fin  area  by  the  accretion  was  found.  Ice  accretions  with  notchings  showed  to  have  a 
favourable  effect  on  stall  onset  as  compared  to  results  received  without  notchings. 

This  is  due  to  the  effect  that  through  the  slots, high  energy  flow  is  transported  into  1 

the  boundary  layer  or  into  the  reattached  bubbles,  respectively,  which  exist  immediate- 
ly behind  the  ice  protrusions.  This  may  also  be  concluded  from  fig.  10b,  where  lift  ver- 
sus drag  for  the  tail  model  has  been  plotted.  Generally,  the  Cp^-value  is  not  effected 
very  much  by  the  ice  but  rather  the  drag  due  to  high  tail  lift  as  occurs  at  operation 

with  flaps  and  slats  out.  1 

Results  of  further  aerodynamic  invertigations : 

1 no  influence  of  ice  accretions  on  flutter  characteristics  of  the  tall  surfaces 
a neglectlble  C.G. -shift  («  0,5  percent  1 p)  due  to  the  additional  weight  of  ice  amount 

(G  . 1.50  kN)  at  the  tail, 
ice 

Flight  test  with  simulated  ice:  simulator  studies  with  aerodynamic  data  taking 
account  for  ice  accretion  showed  that  no  degradatlng  effects  of  any  Importance  should 
be  expected  in  flight.  Nevertheless,  before  flight  testing  under  natural  icing  condi- 
tions was  released.lt  was  decided  to  have  a flight  program  with  simulated  ice.  The  shapes 
were  chosen  according  to  to  fig.  9,  But  because  of  the  aerodynamic  results  of  the  wind 


4-4 


tunnel  tests  it  was  decided  not  to  incorporate  the  lobster  tail  cutting  (see  fig.  11 ). 
Consequently,  these  flight  tests  have  been  performed  under  more  severe  conditions  than 
would  be  expected  under  natural  icing. 

The  scope  of  investigation  comprised 

• forward  and  aft  C.G.  locations 

• low  speed  tests  in  the  clean,  take-off  and  landing  configuration,  the  latter  inclu- 
ding one  engine  failure,  down  to  V = V + 5-6  kts 

• high  speed  tests  up  to  V „ and  M = 0,855 

At  all  these  tests, no  instability  or  loss  in  control  has  occured,  no  sensible  degrada- 
tion of  flying  qualities  was  observed. 

Flight  test  under  natural  icing  conditions:  between  March  26  and  April  9»  1973  alto- 

gether^ flights  ( 36  flight  hours)  under  natural  icing  conditions  (5  flight  hours)  were 
performed  with  the  prototype  D02  of  the  A300.  One  part  of  the  tests  was  carried  out  over 
the  Mediterranean  Sea  between  the  Northcoast  of  Tunesia,  Sardinia  and  Southern  Italy,  the 
other  part  over  the  North  Sea  and  the  North-West  of  Germany  111] • The  mean  flight  alti- 
tude was  8.000  s H ss  15.000  ft  at  ambient  temperatures  of  -3  C T ® s:  -l6°C  and  flight 

speeds  of  l40*  kts  £ V a 270  kts.  The  liquid  water  content  as  determined  during  the 
flights,  fluctuated  between  0.4  g/m^  £ LWC  ± 1.0  g/m^,  the  ti  me  periods  of  icing  were 
10  min.,  droplet  diameter  has  not  been  recorded.  Fig.  12  shows  a photo  of  a specific 
deposit  on  the  vertical  tail  found  during  flight  tests.  In  all  cases, the  maximum 
assumed  ice  accretion  was  below  three  inches.  Flight  qualities  of  the  aircraft  during 
ice  accretion  were  rated  as  very  satisfactory  by  the  pilots. 


2.  UNPROTECTED  APU- AIR -INTAKE  SYSTEM  OF  THE  A300 

The  demonstration  of  the  operational  ability  of  the  APU-air-intake  system  under  icing 
conditions,  postulating  a save  APU  function  and  no  hazardous  performance  degradation, 
also  required  thorough  theoretical  and  experimental  investigations  [2] , [12] , since  the 
APU  was  installed  to  cover  various  supply  cases  not  only  on  ground  but  also  in  flight. 

Critical  intake  components:  as  shown  in  fig.  13.  the  APU-air-intake  system  consists  of 
two  separate  intakes,  the  main  intake  supplying  the  APU  and  the  second  intake  supplying 
the  pneumatic  APU-cooling  system  equipped  with  a fan  cooling  the  APU-oil  cooler,  -gene- 
rator, compartments  and  -engine  shroud.  The  components  for  which  icing  was  expected  are 
the  intake  lip,  guide  vane,  plenum  splitters  and  spherical  protection  grid  for  the  main 
intake  and  the  flat  protection  grid  for  the  cooling  fan  intake. 


Theoretical  prediction  of  ice  accretion  during  flight:  the  ice  built-up  during  flight 
operation  of  the  APU  was  obtained  on  the  basis  of  impingement  analysis  calculating  the 
trajectories  of  the  droplets  reaching  the  APU-intakes.  This  was  examined  for  the  same 
flight  cases  and  icing  conditions  used  for  the  tail  surfaces.  Fig.  l4  shows,  how  a 
"LANGMUIR" -D-droplet  distribution  assumed  in  the  flow  infront  of  the  A300  is  changed  by 


the  local  curvature  of  the  air  flow  on. the  fuselage,  i.e.  on  the  forebody  and  the  tail 


sweep.  The  average  water  contend  in  the  air  ingested  by  the  APU-intakes  is  three  times 


lower  as  for  the  free  stream  condition. 


The  resulting  ice  accretions,  here  shown  in  fig.  15  for  the  critical  components  of  the 
APU  main  Intake,  are  negligible  in  amount  and  height,  causing  no  aerodynamic  implica- 


tions. This  has  been  substantiated  by  the  flight  tests  under  natural  icing  conditions 
[11]  mentioned  before,  where  the  essential  intake  component  (spheric  protection  grid) 
was  controlled  by  a movie  camera.  No  ice  accretion  was  found  at  all. 


<perimental  determination  of  ice  accretion  during  ground  operation:  the  demonstra- 
• of  the  save  function  of  the  APU-air-intakes  under  natural  icing  conditions  on 

g' ouud  (icing  conditions:  LWC  = 0,66  g/m^,  d = 30  pm  and  Ta  = -3°C)  was  performed  by 
expriaent,  determing  the  ice  accretion  and  its  influence  on  intake  characteristics. 

For  'his, two  icing  tests  series  have  been  carried  out  at  the  icing  facilities  of 
Bundesforschungsanstalt  Arsenal  at  Vienna  using  fu^l  scale  air  intakes  in  the  original 
design  configuration  [13]  , [l4]  . Correct  intake  mass  flows  simulated  by  fans  and  exact 
simulation  of  the  intake  flow  field  were  guaranteed.  The  ice  depositions  for  the  main 
intake  after  kO  minutes  are  shown  in  fig.  16.  The  most  interesting  ice  accretion 
was  generated  at  the  spheric  grid,  with  a compact  lump  of  glaze  ice  on  the  side  facing 
the  intake  hole,  and  with  ice  depositions  on  the  rest  of  the  sphere  growing  on  the  grid 
wires  up  to  maximum  height  of  20  ram  without  blocking  the  grid  meshes  and  giving  an  in- 
sight in  the  flow  field  infront  of  the  grid. 

For  the  cooling  fan  intake  ice  accretion  was  only  found  on  the  flat  protection  'rid,  i.e. 
nearly  all  the  water  was  captured  by  the  grid.  The  photo  series  in  fig.  17  shows  the 
temporal  ice  built-up  on  the  grid.  The  grid  center  was  completely  closed  within  a time 
of  15  minutes  leaving  an  open  ring  of  30  percent  of  the  grid  area  which  then  remained 
constant  until  the  end  of  icing  period.  The  final  glaze  ice  block  reached  a height  of 
about  70  mm  with  a deep  excavation  at  the  stagnation  point. 


Effects  of  ice  accretion  on  Intake-  and  system-performance:  the  Influence  of  the  Just 
described  ice  accretions  on  the  two  most  important  intake  parameters  as  ram  recovery  and 
flow  distortion  has  been  measured  during  the  icing  tests  after  each  icing  phase  using 
invertable  pitot-rakes.  Fig.  18  shows  a comparison  of  ram  recovery  and  distortion 
index  DI  without  and  with  ice  after  60  minutes  ground  icing  for  the  main  intake  for 
three  APU  supply  cases.  It  is  demonstrated  that  the  loss  in  ram  recovery  and  the  increase 
of  distortion  index  due  to  icing  is  not  critical,  since  for  all  cases  the  limits  qpMN 
and  DI  s 5 S as  required  by  the  APU  manufacturer  are  not  exceeded,  i.e.  the  minimum 
allowable  Intake  performance  is  guaranteed  [2] . 
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Analysis  for  the  cooling  fan  Intake  [2]  shows  (see  fig.  19) • that, after  15  minutes 
icing, the  intake  pressure  loss  due  to  ice  accretion  is  stabilizing  around  10  percent 
without  reaching  the  maximum  allowable  value  of  15  percent.  This  additional  pressure 
losses  cause  only  neglectible  mass  flow  reductions  without  any  influence  on  the 
cooling  capacity  of  the  system.  The  increase  of  distortion  index  is  unimportant. 

Ground  icing  tests  under  natural  icing  conditions:  the  demonstrations  for  the  APU 
intakes  have  been  completed  by  extensive  ground  icing  tests  with  the  second  prototype 
of  the  A300  at  Helsinki  in  Finland  Cl5].  The  ambient  temperature  found  was  between 
-1°C  a T®  a -15  C and  the  time  period  of  icing  was  several  hours.  Liquid  water  content 
LWC  and  droplet  diameter  have  not  been  measured.  The  fig.  20  shows  two  photos  of  ice 
accretions  found  on  the  main  intake  and  the  cooling  fan  intake.  No  Influence  on  APU- 
and  cooling  system  function  have  been  ascertained  during  the  tests.  According  to  the 
demonstration  procedure,  the  APU  system  of  the  A300  was  certificated  for  one  hour 
ground  operation  under  natural  icing  conditions  and  for  a following  in-flight  opera- 
tion with  iced  intake  components. 

3.  VFW  6 1 4 ICE  SHEDDING  INVESTIGATIONS 

Investigations  of  effects  of  ice  accretion  on  aircraft  components  are  not  complete 
without  looking  into  the  "post  icing  phase".  This  is  the  phase  where  the  ice-material 
bond  is  melting  due  to  heating  effects  (for  example  dynamic  temperature  rise)  and  ice 
shedding  starts.  The  possibility  that  larger  ice  particles  which  shed  from  an  unpro- 
tected aircraft  component  could  be  a hazard  for  the  aircraft  or  any  other  components 
has  to  be  fully  analysed,  i.e.  determination  of  shape  and  size  of  the  particles  and 
its  corresponding  trajectories.  Such  investigations  have  been  performed  for  the  Ger- 
man transport  aircraft  VFW  6l 4 for  which, due  to  the  exceptional  cruise  engine  installa- 
tion, as  shown  in  the  fig.  21,  ice  particles  breaking-off  from  the  unprotected  radome 
have  been  expected  to  be  hazardous  for  the  engine.  The  trajectories  of  the  ice  partic- 
les * have  been  determined  in  a low  speed  wind  tunnel  at  Bremen  with  a l:l6  scale  model 
ejecting  the  particles  from  the  radome  (dynamic  similarity  was  present  - " Froude-No . " - 
Scaling)  [l6J . The  results,  matched  by  theory,  show  that, for  all  flight  attitudes, the 
ice  particles  do  not  meet  the  engine  intakes;  an  example  is  given  on  the  photo  in 


21 . 
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ICING  CONDITION;  LWC  = 0.66gr/m3.  d=30pm.  Ta>  = -3°C 
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A THEORETICAL  AND  EXPERIMENTAL 
MEANS  TO  PREDICT  ICE  ACCRETION  SHAPES 
FOR  EVALUATING  AIRCRAFT  HANDLING  AND 
PERFORMANCE  CHARACTERISTICS 


By 

Ramon  W.  Wilder 

The  Boeing  Company 
Seattle,  Washington 


SUMMARY 

To  determine  the  need  lor  airfoil  ice  protection  systems  and  the  effects  of  large  ice  shapes  on  airplane  performance,  Boeing 
initiated  .1  research  program  lo  obtain  basic  ice  accretion  and  ice  shedding  data  on  typical  jet  transport  swept  airfoils.  This  program, 
sponsored  by  the  Federal  Aviation  Agency,  was  conducted  in  the  NASA  6-foot  by  9-foot  icing  tunnel  at  Cleveland,  Ohio. 

An  empirical  relationship  derived  from  basic  airfoil  water  catch  parameters  was  used  to  correlate  measured  ice  accretion  rates 
with  theoretical  water  impingement  parameters.  This  correlation  used  glaze  ice  heights  and  angles  measured  from  the  plaster  casts 
of  the  actual  ice  shapes  as  taken  from  the  two  swept  airfoil  test  models.  Results  of.this  correlation  show  that  ice  accretion 
characteristics  are  dependent  on  airfoil  shape,  particularly  leading  edge  radius,  camber  and  angle  of  attack. 

Complex  trends  of  the  data  obtained  in  this  test  program  precluded  a general  ice  accretion  relationship  with  other  airfoils. 
Additional  testing  of  other  airfoil  shapes  and  angle  of  attack  would  provide  for  broader  application  of  the  ice  cap  calculation  procedures 
developed  herein. 

Airplane  performance  penalties  associated  with  icing  in  terms  of  landing  weight  penalties,  and  when  these  penalties  are  assessed, 
are  also  discussed.  Destination  airport  temperatures  and  ice  shedding  characteristics  are  shown  to  be  significant  in  determining  the 
frequency  of  aerodynamic  penalties  due  to  ice. 
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NOMENCLATURE 

SYMBOL  UNITS  DESCRIPTION 


A 

Ft2 

Ice  Cap  Cross  Sectional  Area 

C 

Ft 

Airfoil  Chord  Length 

dm 

Microns 

Drop  Diameter 

''m 

Dimensionless 

Total  Water  Catch  Efficiency 

Dimensionless 

Cloud  Extent  Factor 

II 

Feet 

Airfoil  Projected  Height 

"o 

Feet 

Airfoil  Projected  Height  at  0°  Angle  of  Attack 

Ha 

Feet 

Airfoil  Projected  Height  at  Angle  of  Attack 

1 

— 

Icing  Parameter  Defined  by  Equation  7 

K0 

Dimensionless 

Modified  Inertia  Parameter  Defined  in  References  2 and  5 

LS 

Airfoil 

Airfoil  Lower  Surface 

ST 

— 

Airfoil  Stagnation  Point 

°F 

Free  Stream  Air  Temperature 

US 

— 

Airfoil  Upper  Surface 

V 

Knots 

Free  Stream  Velocity 

v„ 

— 

Air  Velocity  Perpendicular  to  the  Airfoil  Leading  Edge 

Wp 

Lb/Hr-Ft^ 

Local  Water  Catch  Rate 

K 

Lb/Hr-Ft  Span 

Total  Water  Catch  Rate 

71 

Lb/Ft-3 

Ice  Density 

*i 

Feet 

Ice  Thickness 

n 

— 

Freezing  Fraction 

e 

Hours  or  Minutes 

Icing  Time 

u 

Grams/Cubic  Meter 

Cloud  Liquid  Water  Content 

A 

Degrees 

Airfoil  Sweet  Angle 

AS 

Feet 

Surface  Distance  Between  Upper  and  Lower  Water  Drop 
Impingement  Limits 

0 

— 

Local  Water  Collection  Efficiency 

INTRODUCTION 

Jel  i-npitit'  aircraft  are  less  susceptible  to  icing  for  several  reason*.  The  kinetic  temperature  rite,  due  to  higher  flight  speeds,  reduces 
the  icing  temperature  envelope  and  prevents  ice  from  forming  in  the  higher,  liquid  water-content  clouds  occurring  just  below  32°  P. 

The  aerodynamic  characteristics  of  swept  airfoils,  use  of  movable  horizontal  stabilizers  from  trim,  rapid  rate  of  climb  and  descent  of 

I 


jet  aircraft,  .ill  contribute  to  minimizing  exposure  to  icing  conditions  and  to  reducing  the  effects  of  ice  on  airplane  performance.  In 
addition,  the  size  of  current  jet  transports  and  the  large  power  reserve  of  jet  engines  reduce  the  drag  effects  of  airplane  icing. 

As  a ivsult  of  the  greater  tolerance  to  icing  of  jet  aircraft,  there  arises  the  question  of  the  need  for  ice  protection  on  certain  areas 
of  the  airplane.  Inboard  wings  have  not  been  protected  on  some  piston  engine  aircraft.  Military  aircraft  sucb  as  the  B-52  and  KC-135 
have  had  tail  and  wing  ice  protection  systems  removed  as  a result  of  adverse  weather  flight  test  programs. 

It  is  logical,  then,  to  assess  the  need  for  airfoil  ice  protection  in  terms  of  airplane  handling  characteristics  and  performance. 

To  do  this  on  a consistent  basis  requires  a wind  tunnel  and  flight  test  program  whereby  the  effects  of  ice  can  be  simulated  since 
natural  icing  is  not  consistent,  is  difficult  and  time  consuming  to  obtain,  and  is  subject  to  shedding  before  the  testing  can  be 
completed  I bis  simulation  of  icing  effects  has  been  accomplished  on  jet  transports  through  the  use  of  artificial  ice  shapes  attached 
to  unheated  portions  of  the  wing  and  tail  surfaces.  This  paper  summarizes  some  of  the  methods  and  procedures  used  to  determine 
these  ice  shapes  and  presents  a new  method  of  obtaining  ice  shapes  based  on  recent  ice  accretion  test  data  correlated  with  theoretical 
impingement  data. 


ANAI.YI  If  Al.  I’REDJCTION  OF  ICF.  SHAPES 

Ice  shapes  on  airfoils  are  classified  as  "glaze”  or  “rime”  although  many  ice  shapes  will  be  a mixture  of  the  two.  Rime  ice  is 
tumid  at  combinations  of  low  temperatures  and  low  liquid  water  contents.  The  water  droplets  freeze  on  impact  resulting  in  a 
milky-white  ice  shape,  as  shown  in  figure  I , dictated  by  the  airfoil  impingement  characteristics  in  terms  of  water  catch,  distribution, 
and  impingement  limits.  The  drag  and  other  aerodynamic  effects  of  rime  ice  are  much  less  severe  than  glaze  ice  because  of  the 
streamline  form,  therefore,  this  type  of  icing  has  not  been  considered  in  determining  the  airplane  tolerance  to  icing. 


Figure  1 — Typical  Rime  ice  Formation  on  Swept  Airfoil 


Glaze  ice  forms  at  combinations  of  high  liquid  water  contents  and  surface  temperatures  near  freezing.  Not  all  of  the  water 
droplets  freeze  on  impact,  but  some  run  a short  distance  before  freezing,  resulting  in  a blunt  or  double  hom  ice  shape  as  shown  in 
liptirc  2.  This  ice  is  usually  clear  since  there  is  little  or  no  air  entrained  in  the  ice.  The  upper  surface  hom  acts  as  a spoiler,  increasing 
drag  ami  reducing  lift.  This  icc  shape,  because  of  its  "spoiler"  action,  can  affect  aerodynamic  characteristics  in  terms  of  increased 
drag  and  decreased  maximum  lift  coefficient.  The  lower  surface  hom  usually  has  little  effect  because  of  low  local  air  velocities.  On 
highly  swept  airfoils,  glaze  ice  tends  to  form  us  a scries  of  discontinuous  cup  shapes  due  to  the  spanwise  velocity  component. 

The  shapes  and  type  of  ice  is  dependent  on  the  surface  temperature,  which  is  a complex  function  of  airspeed,  liquid  water  content, 
ambient  air  temperature,  local  flow  field  around  the  object,  and  altitude.  A correlation  of  ice  shapes  on  a small  2-inch  diameter 
unswepl  cylinder,  with  these  parameters,  is  shown  in  figure  3 where  q is  the  freezing  fraction  defined  in  reference  3 as  that  part  of 
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■ In:  water  catch  Hut  freezes  where  it  impinges.  In  general  terms,  glaze  ice  is  likely  to  occur  at  total  air  temperatures  between  25°  and 
32°  I between  10°  and  25°  F,  a mixture  of  rime  and  glaze  ice  will  occur  with  glaze  predominant  in  the  stagnation  area  and  for  high 


lupml  w.ilei  content  conditions. 


Figure  2-Typical  Glaze  Ice  Shape  on  Swept  Airfoil 


THUt  Ain  SPEED  (KNOTS) 


Figure  3-Ice  Shape  Types  es  a Function  of  Speed  and  Ambient  Temperature 
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CALCULATION  PROCEDURES 

Procedures  previously  developed  for  determining  the  ice  cap  size  are  based  on  calculating  the  airfoil  water  catch  and  adjusting 
(lu>  water  i .itch  or  ice  thickness  into  a typical  mushroom  icc  shape.  For  rime  ice,  the  icc  thickness  and  shape  can  be  calculated 
directly  from  the  local  water  catch  and  ice  density  since  the  water  freezes  immediately  on  contact  with  the  airfoil.  The  local  ice 
tliichness,  5, , is  related  to  the  local  water  catch  rate,  Wp,  by  the  following  equation  where  9 is  the  time  in  icing  and  7,  is  the  ice  density. 

(1) 

Ti 

Thus,  a theoretical  icc  shape  for  rime  ice  can  be  estimated  by  calculating  the  local  water  catch  rate  for  various  positions  in  the 
impingement  area.  The  equation  for  the  local  water  catch  rate  is  given  below  where  0 is  the  local  water  collection  efficiency,  V the 
free  stream  velocity,  A the  airfoil  sweep  angle  and  co  the  cloud  liquid  water  content. 

W0  = 0.380  (VXCOS  A)  (w)  (2) 

This  procedure,  with  some  modifications,  was  used  for  calculating  the  ice  shapes  for  the  Boeing  707  and  727  series  airplanes.  To 
account  for  glaze  icing  and  the  effects  of  ice  shape  on  water  impingement,  the  maximum  local  collection  efficiency,  0 max,  was 
assumed  to  be  1 .0.  This  means  all  of  the  water  in  the  cloud  in  the  path  of  the  airfoil  would  impinge  on  the  airfoil.  This  assumption  was 
made  to  account  for  the  effects  of  ice  shape  on  water  collection  efficiencies.  This  assumption  is  approximately  correct  for  small 
airfoils  where  the  maximum  local  collection  efficiency  is  near  0.8,  but  becomes  overly  conservative  for  larger  inboard  wings  and 
stabilizers  where  the  collection  efficiencies  are  less  than  0.5.  The  ice  thickness  calculated  using  0 max  equal  to  1.0  was  then  adjusted  to 
give  a typical  mushroom  or  glaze  ice  shape  using  the  known  impingement  limits. 


In  other  procedures,  references  2 and  6,  the  ice  shape  is  determined  from  the  calculated  total  water  catch  rate  (Wc)  from  bare  airfoil 
impingement  data  using  the  following  equation: 


Wc  = 0.38  V (COS  A)  Em  (H/C)  (C)  u 


<3) 


In  this  equation,  Em  is  the  total  water  collection  efficiency  based  on  the  airfoil  projected  height  (H),  and  (C)  is  the  airfoil  chord 
length.  This  water  catch  rate  is  in  terms  of  pounds  of  ice  per  hour  per  foot  of  span.  The  cross  sectional  area  of  the  ice  cap  is  then 
determined  from: 


144  Wc  9 
7i 


(4) 


A double-peaked  glaze  icc  shape,  as  shown  in  figure  4 which  contains  the  calculated  theoretical  cross  sectional  area,  is  then  drawn  by 
trial  and  error. 


46  MIN  HOLD  AT  6000  ft  ALTITUDE 


Figure  4-Typicel  Ice  Shepet  for  Unprotected  Tell  Surf ecet  from  Ref  2 
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The  above  calculation  procedures  only  approximate  the  actual  ice  shape  since  accurate  ice  accretion  data  for  most  airfoils  is  not 
available.  Also,  water  drop  sizes  and  resultant  collection  efficiencies  must  be  checked  for  the  maximum  water  catch.  Figure  5 shows 
a typical  plot  ol'  the  effect  of  drop  size  on  collection  efficiency  and  cloud  liquid  water  content.  This  calculation  will  generally  show 
that  the  maximum  catch  rate  will  occur  with  20  to  2S  micron  droplets. 


G>clouo  liquid  water  content  based  on  icing 
limit  TEMPERATURE  OF  2t°F  CORRESPONDING  TO 
HOLDING  SPEEDS 


MEDIAN  DROPLET  OIA  IMICRONSI 


Figure  5-Effect  of  Drop  Size  on  Water  Catch  Parameters  Holding  Flight  Speeds 


EXPERIMENTAL  DETERMINATION  OF  ICE  SHAPES-NACA  ICE  SHAPE  DATA 

A correlation  of  measured  ice  shapes  from  icing  tunnel  tests  with  theoretical  impingement  parameters  nas  been  reported  in 
reference  4 (NACA  TN  41  51 ) for  an  unswept  symmetrical  4 percent  thick  airfoil.  'Ilte  data  in  this  report  lor  mushroom  ice, 
plotted  m a slightly  different  manner,  is  shown  in  figure  6.  It  indicates  the  measured  maximum  ice  height  is  leu  than  the  theoretical 
height  based  on  the  maximum  local  water  catch.  This  is  not  surprising  since  all  of  the  impinging  water  does  not  freeze  at  the  point  of 
maximum  collection  efficiency  under  these  icing  conditions.  Unfortunately,  there  arc  no  similar  published  data  for  highly  cambered 
swept  airfoils.  Therefore  Boeing  initiated  a research  program  to  obtain  basic  ice  accretion  data  on  two  airfoils  representative  of  the 
inboard  wmg  and  horizontal  stabilizer  configuration  of  typical  current  jet  transport  aircraft.  This  program  was  conducted  in  the  NASA 
6-foot  by  9-foot  icing  tunnel  at  Cleveland,  Ohio. 
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Figure  6-Corralatlon  of  fee  Haight  with  Theoretical  Impingement  Unswept  4%  Thick  Airfoil 


BOEING  ICE  ACCRETION  TEST  SETUP  AND  PROCEDURE 


GEOMETRIC  CHORD  LINE 
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Figure  7- Airfoil  Cross  Section  Streamwise  Direction 


Tlic  models  used  in  this  test  program  were  chosen  to  represent  the  type  of  airfoils  used  on  present-day  Boeing  commercial  jet 
uircrjll.  Cross  sections  of  these  airfoils  are  shown  in  figure?.  One  model,  designated  BAC  450,  was  typical  of  an  inboard  wing  section 
with  a 1 2.5  percent  thickness  ratio  and  27.5  degree  sweep  angle.  The  other  model,  designated  BAC  470,  represented  outboard  wing 
and  horizontal  stabilizer  sections  with  an  8 percent  thickness  ratio  and  40  degree  sweep  angle.  These  models,  shown  in  figure  8, 
have  a 6-foot  streamwise  span  and  chord  length  and  were  mounted  vertically  in  the  Lewis  icing  tunnel  at  Cleveland,  Ohio.  Each  model 
was  fabricated  from  aluminum  with  six  thermocouples  and  three  electrical  heater  pads  marked  A,  B,  and  C in  figure  8,  bonded  to  the 
inner  surface  of  the  leading  edge.  The  thermocouples  were  used  to  monitor  the  leading  edge  surface  temperature  during  the  icing  runs 
and  the  heater  pads  were  used  to  assist  in  removing  ice  samples. 


r-BAC  150  AIRFOIL 
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Figure  8-Ice  Accretion  Teat  Models  BAC  460  end  BAC  470 

Icing  conditions  in  the  tunnel  were  obtained  and  measured  according  to  previously  established  techniques  and  calibrations.  Icing 
cloud  factors  considered  in  this  study  were  cloud  liquid  water  content,  icing  time,  airspeed,  and  temperature.  Cloud  droplet  »l/o  could 
not  lie  varied  independently  due  to  the  limitations  in  the  design  of  the  tunnel  spray  system.  As  a result,  the  cloud  droplet  sizes 
increased  with  increased  liquid  water  content  and  decreased  with  an  Increase  In  alra|teed. 


Aii  Kin;.!  run  consisted  of  a 7.5  or  IS  minute  exposure  of  the  model  to  a particular  set  of  icing  conditions  and  airfoil  angle  of  attack. 
On  completion  of  the  prescribed  icing  time,  the  tunnel  was  stopped  and  measurements  taken  of  the  final  ice  shape.  One  or  more  sections 
ol  the  i»e  cap.  approximately  o inches  in  length,  were  removed  using  the  heater  pads  and  the  ice  scraper  used  in  reference  4 tests.  This 
scraper  is  steam  healed  and  has  an  internal  vacuum  chamber  to  assist  in  drawing  off  the  melted  ice  so  as  not  to  affect  the  sample  weight 
or  shape.  I lie  removed  sections  were  then  dipped  repeatedly  in  liquid  paraffin  to  create  a mold  of  the  ice  shape.  When  the  ice  in  the 
mold  melted,  it  was  poured  into  a beaker  and  weighed.  The  resultant  cavity  was  then  refilled  with  tap  water  which  was  then  poured  into 
the  beaker  and  weighed  to  determine  the  specific  gravity  of  the  ice.  Representative  values  of  the  specific  gravity  of  the  ice  cap  measured  in 
this  manner  were  from  0.75  to  0.91. 

A replica  of  the  ice  sample  was  then  made  by  pouring  plaster  of  parts  into  the  mold  and  later  melting  the  wax  in  an  oven.  Photographs 
of  the  equipment  used  to  obtain  plaster  casts  of  the  ice  cap  are  shown  in  figures  9 and  10.  Sample  plaster  ice  cap  photographs  are  shown 
in  figures  1 1 and  12. 


Figure  9 -Test  Equipment  for  Weighing  end  Obteining  Ice  Stupe  Molds 


Figure  10-Meltlng  Fereffln  From  Fleeter  Ice  Cest 
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Figure  1 1 -Coating  Ice  Cap  With  Melted  Paraffin 


Alter  I lie  ice  sample  was  removed  from  the  mold,  the  cron  section  of  the  ice  cap  was  photographed  against  a one-quarter  inch  wire 
inesli  grid  held  normal  to  the  airfoil  leading  edge.  A typical  ice  cap  cross  section  is  shown  in  figure  13.  Due  to  the  difficulty  in  discerning 
the  cross  sectional  outline  from  the  photographs,  all  ice  cap  measurements  were  made  from  the  plaster  casts  as  defined  in  figure  14. 
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Figure  13- Installation  of  Wire  Grid  Over  the  Ice  Cap 
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Figure  14-Ice  Shape  Dimension 


The  majority  of  the  ice  samples  were  of  the  glaze  type  with  typical  double  hom  protuberances.  A few  rime  ice  samples  were  obtained 
and  the  remainder  were  mixtures  of  the  two.  Figure  1 5 shows  representative  ice  shapes  obtained  during  this  test  program.  The  effect  of 
the  winit  swccpback  was  most  noticeable  in  the  discontinuous  cup  shapes  formed  in  the  spanwise  direction  as  shown  in  figures  16  and  17. 
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Figure  1 7-Glue  Ice  Shape  BAC  470  Airfoil  Frontal  View 


The  separations  were  most  prominent  in  the  glaze  ice  formations  and  became  larger  as  the  ice  thickness  increased.  In  those  samples 
which  tended  toward  rime  ice,  the  distance  between  individual  cups  decreased  and  in  some  cases  formed  a continuous  solid  ice  cap  along 
the  span  of  the  airfoil. 

Tlie  leading  edge  temperatures  for  the  glaze  icing  runs  varied  between  25°  and  29°  F.  Due  to  the  relatively  high  thermal  conductivity 
of  the  aluminum  leading  edge  (typical  of  current  aircraft  configurations)  the  chordwise  temperature  gradient  under  the  ice  cap  was  negligible. 
Had  the  leading  edge  been  made  of  an  insulating  material,  the  ice  cap  could  have  had  a significantly  different  shape  since  the  impinging  water 
would  see  a warmer  stagnation  point  surface  temperature.  This  would  allow  the  water  to  flow  further  along  the  surface  before  freezing.  No 
attempt  was  made  in  this  study  to  determine  the  magnitude  of  this  effect  of  leading  edge  material  properties  on  ice  shape. 

TEST  AIRFOIL  IMPINGEMENT  CHARACTERISTICS 

The  theoretical  impingement  characteristics  in  terms  of  water  drop  collection  efficiency  and  impingement  limits  are  shown  in  figure  18 
lor  the  two  angles  of  attack  eonsidered  in  this  study.  These  data  were  obtained  from  a water  droplet  trajectory  computer  program  and  are 
plotted  as  a function  of  a modified  inertia  parameter  K0.  This  parameter,  described  and  defined  in  references  2 and  5,  allows  the 
impingement  data  taken  for  one  flight  condition  to  be  used  for  other  flight  conditions. 


TOTAL  COLLECTION  EFFICIENCY 


IMPINGEMENT  LIMITS 
BAC  450 


*0 

IMPINGEMENT  LIMITS 


BAC  450 


Figure  18-Teet  Airfoil  Impingement  Cheracterittict 
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TlicVrojected  height  of  the  airfoil  as  shown  in  figure  19  was  used  in  this  study  as  the  characteristic  dimension  for  overall  water 
drop  collection  efficiency  and  in  the  ice  shape  correlation  parameters.  This  is  in  line  with  the  impingement  data  as  presented  in  reference  2. 


Figure  19— Airfoil  Projected  Height  BAC  450  and  470 


ICE  SHAPE  CORRELATION 


Correlation  of  the  measured  ice  cap  dimensions'  was  based  on  the  water  catch  parameters  which  are  indicative  of  the  ice  thickness. 
Under  conditions  of  glaze  icing,  the  ice  thickness  is  more  a function  of  total  water  catch  or 


where  the  total  water  catch  is  calculated  from 


(5) 


Wc  = 0.38  V (COS  A)  Em  (C)  w (6) 

And  AS  is  the  chordwise  surface  distance  measured  from  the  lower  water  drop  impingement  limit  to  th*  upper  water  drop 
impingement  limit.  An  icing  parameter  indicative  of  the  ice  thickness  in  terms  of  the  water  catch  can  be  defined  from  the  relationship  in 
equation  (5)  as: 


Wc  6 
TjAS 


« Hi  V (COS  A)  — ) (u>)8 

Tj  AS  1C  \cj 


(7) 


P 

The  ice  density  under  glaze  icing  can  be  assumed  relatively  constant,  and  further,  the  term  m/AS.  >*  * single  valued  function  for 
a given  K0  value.  Since  AS  is  difficult  to  measure  or  obtain  accurately  and  is  a particular  value  for  C En,.  only  the  total  collection 
efficiency  was  used  in  this  correlation.  Thus,  the  icing  parameter  indicative  of  a theroetical  average  ice  thickness  reduces  to: 

' ‘^S  ’ (8) 

Where  X is  the  exponential  factor  used  to  correlate  angle  of  attack  effects. 

Measured  ice  thicknesses  in  terms  of  the  stagnation,  and  upper  and  lower  ice  cap  protuberances  are  plotted  against  the  icing  parameter 
defined  by  equation  (8)  in  figures  20  and  21  for  both  the  BAC  470  and  450  airfoils.  The  effect  of  angle  of  attack  was  resolved  for  the  BAC 
450  airfoils  by  raising  the  projected  height  to  the  fifth  power.  No  factor  was  found  necessary  for  the  BAC  470  airfoil  since  the  projected 
height  to  chord  ratio  (H/C)  provided  sufficient  angle  of  attack  correction.  These  factors  are  emperical  for  the  range  of  the  test  data  or  from 
0 to  4.5  degrees.  However,  some  extrapolation  beyond  these  limits  is  felt  to  be  valid  and  within  the  accuracy  of  the  test  data. 
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These  curves  define  the  ice  thickness,  but  not  their  relative  positions.  The  angles,  which  the  upper  and  lower  pinnacle  dimensions 
made  with  the  geometric  chord  line,  were  measured  and  correlated  in  a similar  manner.  This  correlation  is  shown  in  figure  20  for  the  470 
airfoil.  These  data  were  not  as  orderly  as  the  ice  thickness  data,  but  they  provide  a definite  trend  which  is  usefi'1  in  defining  the  ice  shapes. 

It  is  mil  iccoiiimciidcd  that  this  ice  angle  data  be  extrapolated  much  beyond  the  airfoil  angles  of  attack  from  which  the  data  was  obtained 
without  correlating  impingement  data. 

Tlie  frontal  heights  obtained  through  these  ice  angles  are  maximum  dimensions  since  the  test  data  was  obtained  at  or  near  maximum 
icing  temperatures  where  the  impinging  water  can  run  farther  aft  before  freezing.  The  range  of  flight  and  icing  conditions  during  holding  or 
approach  are  within  the  test  icing  conditions,  so  these  dimensions  will  also  be  valid  for  flight  conditions. 

Some  significant  trends  are  readily  apparent  from  the  curves  obtained  by  this  correlation  (figures  20  and  21 ).  The  most  significant  of 
these  is  the  distinct  reduction  in  slope  which  occurs  once  the  ice  attains  a thickness  between  1-1/2  to  2 inches.  This  change  in  the  slope  is 
indicative  of  a reduced  rate  of  growth  for  the  ice  shapes.  This  can  be  explained  by  the  physics  of  the  ice  accretion.  As  the  ice  pinnacles 
associated  with  the  glaze  ice  shape  builds  to  these  heights,  the  droplets  impinging  near  the  stagnation  area  become  trapped  and,  assisted  by 
the  spanwise  velocity  component  of  the  swept  wing,  freeze  within  the  pinnacles  instead  of  contributing  to  additional  growth  pinnacles.  This 
is  apparent  in  the  more  blunt  ice  shapes.  For  the  larger  ice  shapes,  the  ratio  of  stagnation  ice  thickness  to  pinnacle  heights  approaches  unity 
as  shown  by  the  following  tabulation: 


Icing  Parameter,  I 

(6  ST)„ 

(6US)„ 

(6ST)^/(6US)„ 

3 

.54 

.76 

.710 

8 

1.5 

1.75 

.857 

20 

2.33 

2.61 

.893 

A comparison  of  the  curves  for  the  two  airfoils  also  illustrates  the  influence  of  camber  and  leading  edge  radius  on  the  maximum  ice 
thickness.  The  BAC  450  airfoil,  being  blunter,  has  a larger  actual  ice  accretion  area.  This  resulted  in  a larger  ice  cap  both  in  thickness  and 
area  over  the  less  blunt,  highly  cambered  BAC  470  airfoil. 


APPLICATION  OF  ICE  SHAPE  CORRELATION  CURVES 

Knowing  the  ice  thicknesses,  angles  and  impingement  limits,  an  ice  shape  can  be  constructed  which  more  closely  duplicates  icing  of 
swept  wing  aircraft.  Knowledge  of  the  water  drop  trajectories  impingement  limits  and  stagnation  point  will  be  illustrated  for  the  747 
horizontal  stabilizer  (BAC  470  airfoil)  shown  in  figure  22.  Also  shown  on  this  figure  is  the  707  horizontal  stabilizer.  The  following 
flight  design  condition  will  be  used,  which  is  representative  of  that  used  on  early  707  ice  protection  system  deletion  studies: 


747 


n Compenaon  of  747  and  707  Horizontal  Stabilizer  Sizes 
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FLIGHT 

0 = 30  minute  hold 

ALT  = 15,000  feet  altitude 

V = 377  inph  True  Airspeed  (256  knot' indicated  airspeed) 

ICING  CONDITIONS  (FAR  PART  25,  APPENDIX  C) 

tQ  = 2 1°F  Icing  Limit  temperature  (ambient  air  temperature) 

dm  = 1 5 Micron  drop  diameter  (“A”  drop  distribution) 

F = 150  Mile  cloud 

to  = 0.266  g/m^  cloud  water  content  ( 1 50  mile  cloud) 

AIRFOIL  AND  IMPINGEMENT  DATA 


Stab. 

Angle  of 

Sweep 

Chord 

Ko 

(H/C)» 

(H/C) 

^m 

Sta. 

Attack 

Angle 

Length 

(BAC  470) 

(Fig.  19) 

(Fig.  19) 

(Fig.  19) 

1*17.5 

-2.38‘> 

43° 

280  in. 

.005 

.087 

.095 

.020 

410 

-2.38° 

43° 

280  in. 

.0109 

.087 

.095 

.051 

'Projected  height  for  0 degree  angle  of  attack. 

Knowing  the  impingement  and  flight  conditions,  the  icing  parameter,  I,  can  now  be  calculated  for  the  horizontal  stabilizer  station  SBL 
410  as: 

I V(COSA)(IM(h>)(Em)ij- 

- 377  (cos  43°)  (30)  (0.266)  (.05 1 ) (.095) 

1 = 10.66 


From  figure  20,  the  ice  pinnacle  heights  are: 


Stab. 

Icing 

Ice  Thickness 

(Normal  to  L.E.) 

Ice  Thickness  Streamwise 

Sta. 

Par 

Upper 

Stag. 

Lower 

Upper 

Stag. 

Lower 

(SUS)„ 

(6ST)„ 

(6LS)„ 

(«US)s 

(6ST)s 

(8LS)j 

197.5 

4.1 

1.0 

.75 

.99 

1.37 

1.025 

1.36 

410. 

10.66 

2.03 

1.84 

2.10 

2.78 

2.52 

2.78 

Tlte  ice  angles  can  now  be  calculated  from  the  ratio  of  ice  thicknesses  and  figure  20. 

Upper 

Upper  (1) 

Lower 

Lower  (2) 

Stab. 

Projected 

Ice  Thickness 

Ice  Angle 

Ice  Angle 

Ice  Angle 

Ice  Angle 

Sta. 

Ht.  Ratio 

Ratio 

Par. 

(Fig.  20) 

Par. 

(Fig.  20) 

/£st\  /6st 

\ /S 

1?I\  (MV 

/6ST\  /Ha> 

2 

VHo 

VUS/t,  V1!? 

'v  ' 

US/,  Uo / 

V^LS/j)  \ Hc / 

1 

197.5 

1.09 

.75  .76 

.89 

46° 

.63 

-44° 

410 

1.09 

.908  .878 

1.08 

26° 

.74 

-38.5° 

( 1 ) Measured  clockwise  from  geometric  chord 

(2)  Measured  counterclockwise  from  geometric  chord. 

From  the  above  ice  angles  and  pinnacle  heights,  the  resultant  ice  shape  is  drawn  as  shown  in  figure  23.  The  corresponding  nominal 
3 inch  707  ice  shape  is  also  shown  for  reference  baseu  on  the  previous  method  whereby  all  of  the  impinging  water  was  assumed  to  collect 
on  the  airfoil.  Reasonably  good  agreement  between  the  two  methods  is  shown  at  the  tip. 

However,  due  to  the  lower  collection  efficiency  at  the  inboard  end  of  the  stabilizer,  the  latest  calculation  procedure  gives  a much 
smaller  ice  thickness.  This  is  in  line  with  observations  on  natural  icing  flight  tests.  The  difference  in  the  location  of  the  ice  caps  is  due 
to  the  higher  angle  of  attack  of  the  less  cambered  707  airfoil.  This  angle  of  attack  is  determined  from  airplane  trim  requirements. 

The  mclliod  of  determining  ice  shapes  developed  from  the  Cleveland  test  data  is  applicable  to  any  airfoil.  Extrapolation  from  this 
data  will  generally  be  satisfactory  if  the  airfoils  are  similar.  Consideration  should  be  given  to  obtaining  ice  accretion  data  if  the  airfoil  is 
considerably  different  than  that  tested  here,  or  the  surface  is  considered  particularly  critical  from  a flight  safety  or  performance  standpoint. 

No  attempt  was  made  in  this  test  program  to  obtain  ice  shapes  for  extended  leading  edge  devices  such  as  flaps  or  slats  since  these  are 
normally  retracted  during  icing  conditions.  Ice  shapes  on  extended  leading  edge  devices,  if  required,  can  be  and  have  been  determined  from 
estimated  water  collection  rates. 
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Figure  23-747/707  Horizontal  Stabilizer  ice  Shapes 


SOME  C OMMENTS  ON  FLIGHT  TESTING  WITH  SIMULATED  ICE  SHAPES 
WIND  TUNNEL  TESTING 

A preliminary  evaluation  of  the  effects  of  ice  on  airplane  performance  and  stability  is  usually  obtained  from  wind  tunnel  testing  with 
ice  shapes  attached  to  the  indicated  airfoil  sections.  Wind  tunnel  test  results  showing  the  effect  on  the  lift  coefficient  and  pitching  moment 
curves  of  a typical  Boeing  jet  transport  are  shown  in  figure  24.  At  angles  of  attack  associated  with  level  flight,  the  ice  has  little  effect  other 
than  to  increase  drag.  This  effect  on  airplane  drag,  with  flaps  extended,  is  shown  in  figure  25.  As  indicated  in  this  figure,  the  drag  increase 
due  to  ice,  increases  sharply  with  landing  flap  settings.  It  is  under  these  conditions  where  ice  has  the  most  effect  on  airplane  performance. 


•-ICE 
O-NO  ICE 


Figure  24- Effect  of  lea  on  Lift  Drag  and  Pitching  Moment  fora  Typical 
Jat  Transport  From  Wind  Tunnel  Data  t Flaps  Up) 
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Figure  25-Effect  of  Ice  on  Dreg  for  e Typical  Jet  Transport  I Flaps  Extended ) 


FLIGHT  TESTING 

The  final  check  of  the  airplane  tolerance  to  icing  is  usually  obtained  through  flight  testing  with  simulated  ice  shapes  unless  previous 
data  is  available  or  sufficient  to  determine  aerodynamic  characteristics  or  penalties.  These  simulated  ice  shapes  are  made  up  from  a styrofoam 
core  attached  to  a fiberglass  glove  which  fastens  over  the  airfoil  leading  edge.  Epabond  or  body  putty  is  applied  over  the  styrofoam  ice 
shape  to  give  a rough  outside  texture.  A typical  ice  shape  used  on  the  727  horizontal  stabilizer  is  shown  in  figure  26.  Exact  simulation  of  the 
cup-like  protrusions,  shown  in  figure  1 7,  is  not  considered  necessary  since  the  hom  or  ice  pinnacle  acts  as  the  spoiler  which  results  in  airflow 
separation.  In  some  conditions,  the  cup-like  protrusions  could  act  like  vortex  generators  and  promote  reattachment  of  the  airflow  to  the 
airfoil.  Therefore,  only  the  external  ice  shape  is  duplicated,  including  the  sharpness  of  the  ice  pinnacles. 

The  flight  test  evaluation  by  FAA  pilots  includes  basic  airplane  performance  data  as  well  as  flight  characteristics  for  the  following  at  the 
most  critical  center  of  gravity: 

1.  Demonstration  of  stalls,  sideslip,  heading  changes  under  engine-out  conditions  and  windup  turns. 

2.  Demonstration  of  elevator  effectiveness  through  pushovers  starting  at  1.4  Vs  to  flap  placard  speeds  and  down  to  1.1  V,. 

3.  Demonstration  of  high  speed  handling  characteristics  including  lateral  control  evaluation,  pushovers,  and  pullups. 

4.  Flap  retraction  and  extension  during  simulated  go-around. 

5.  Demonstration  of  approach  and  cruise  stability  through  trimming  the  airplane  and  stabilizing  on  airspeed  through  elevator  only. 


Rgure  26-727  Artificial  toe  Shape  For  Horizontal  Stabilizer 
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This  biter  condition  usually  is  the  most  difficult  in  that  the  aircraft  must  have  sufficient  elevator  control  available  to  attain  speeds 
between  I . I Vs  and  1 .7  Vs  for  approach  without  retrimming  the  aircraft.  Using  a movable  horizontal  stabilizer  for  trim  usually  means  | 

more  elevator  control  available,  especially  for  fully  powered  or  directly  controlled  elevators.  Aircraft  with  fixed  stabilizers  have  exper- 
ienced severe  nose  down  pitching  due  to  loss  of  elevator  control  from  the  effects  of  ice  on  the  stabilizer  leuding  edges.  This  is  usually 
associated  with  a change  in  flap  positions  or  speed,  which  requires  the  tail  surfaces  to  work  harder.  These  surfaces,  already  being  at  a 
high  lift  coefficient  condition,  "stall  out”  resulting  in  loss  of  elevator  effectiveness. 

f 

PERFORMANCE  PENALTIES 

In  addition  to  evaluating  flight  characteristics,  airplane  performance  must  also  be  considered  due  to  the  ice  effects  on  drag.  Figure 
27  shows  a typical  drag  polar  obtained  from  flight  tests  on  a Boeing  jet  transport.  These  data,  along  with  wind  tunner  test  data,  are  used 
to  evaluate  the  flaps-down  approach  and  landing  climb  performance  as  required  for  use  in  the  flight  manual.  Figure  28  shows  the  effect  of 
ice  in  terms  of  climb  gradient  and  aircraft  gross  weight.  The  increment  between  the  ice-on  and  ice-off  curves  must  be  subtracted  from  the 
aircraft  maximum  allowable  landing  weight.  On  short-haul  aircraft,  where  the  maximum  landing  weight  is  near  the  takeoff  weight,  this 
weight  decrease  can  come  out  of  the  payload.  This  performance  penalty  has  been  assessed  any  time  icing  conditions  en  route  are 
anticipated,  and  the  destination  airport  temperature  is  below  the  temperature  at  which  the  ice  will  shed. 


Figure  27-Effect  of  Ice  on  Low  Speed  Dreg  Polars-25°  Flaps— All  Engines  Operating 


NOTES 

1.  APPROACH  CLIMS  CONFIGURATION 

A 26*  T.  E.  FLAPS 

b.  ONE  ENGINE  OUT 

c.  GEAR  UP 

tf.  CL  • (I/O)  MAX 

2.  LANDING  CUM!  CONFIGURATION: 

A <0*  T.  E FLAPS 
tt  GEAR  OOWN 


t iguro  28  Effoct  of  Ice  on  Approach  and  Landing  Climb  Gradient-7000  Ft  Altitude' 
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ICt  SHtUUINC 

Ice  shedding  characteristics  can  also  be  of  interest  in  determining  the  need  for  airfoil  ice  protection  since  they  determine  the  frequency 
ol  application  of  I lie  airplane  performance  penalty.  This  penalty,  in  terms  of  the  airplane  maximum  landing  weight,  is  presented  in  the 
performance  section  of  the  I AA  flight  manual  in  terms  of  destination  ambient  air  temperature.  To  determine  this  maximum  allowable  air 
temperature  requires  a knowledge  of  the  airplane  descent  flight  conditions,  the  resultant  airfoil  surface  temperature  under  the  ice  cap  and  the 
surface  tumperature  at  which  the  ice  sheds. 


Typical  jet  transport  time-temperature  surface  temperature  relationships  for  a descent  flight  profile  are  illustrated  in  figure  29.  As  the 
airplane  stalls  the  descent  phase  from  cruise  conditions,  the  ambient  air  temperature  starts  to  increase  due  to  the  temperature-altitude  lapse 
rate  as  shown  by  the  lower  line  in  figure  29.  The  total  temperature,  upper  line,  follows  directly  as  a function  of  velocity.  The  surface 
temperature  starts  to  follow  the  air  temperatures  but  lags  due  to  the  thermal  inertia  of  the  ice  and  airplane  leading  edge.  At  some  point  in  the 
descent,  the  frictional  heating  (ram  temperature  rise)  is  sufficient  to  melt  the  ice-metal  bond  and  the  ice  sheds.  During  this  time,  the  surface 
temperature  either  remains  constant  or  climbs  slowly,  due  to  the  latent  heat  of  fusion  being  supplied  to  the  ice-metal  interface.  After 
shedding,  there  is  an  abrupt  rise  in  surface  temperature  as  insulation  due  to  the  ice  is  lost  and  the  airfoil  is  heated  by  aerodynamic  heating 
effects. 


TEMP 


Figure  29— Typical  Time-  Temperature  Profile  During  Descent 


Tlie  airfoil  surface  temperature  under  the  ice  cap  must  be  known  in  order  to  determine  the  time  at  which  the  ice  will  shed  during  a 
descent.  This  temperature  differs  from  the  total  of  recovery  temperature  during  descent  because  of  the  ice  cap  thermal  capacitance  (thermal 
lag).  A method  of  predicting  the  time  the  ice  will  shed  has  been  developed  by  Boeing,  using  laboratory  test  data  in  conjunction  with  a 
heal  Iraiixfci  analysis.  Tills  method,  reference  2,  has  been  verified  using  datu  from  a 737  natural  icing  flight  test. 

CONCLUSIONS 

The  determination  of  the  need  for  airfoil  ice  protection  requires  the  establishment  of  the  effects  of  large  ice  accretions  on  airplane 
IH'rfnrmautv.  The  Ice  accretion  test  program  conducted  by  Boeing  provides  a basis  for  accurately  determining  the  sire  and  shape  of 
these  ice  accretions.  This  has  been  accomplished  by  correlating  the  measured  ice  cap  dimensions  with  theoretical  impingement  parameters. 
Use  of  this  correlation  will  enable  the  calculation  of  artificial  ice  shapes  representative  of  those  which  would  accrete  on  the  airplane  during 
a severe  icing  encounter. 
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1.  EUROPEAN  INQUIRY  ON  ICING  FACILITIES  AND  TECHNIQUES 

The  European  inquiry  was  conducted,  as  suggested  by  the  R.T.D.  chairman,  by  Marcel  Pierre,  previously  Director 
of  the  Modane  Wind  Tunnels  Department  and  presently  Senior  Technical  Advisor  at  ONERA  for  new  facility  projects. 
The  inquiry  results  are  presented  in  tabular  form. 

1.1.  Organizations  consulted 

Table  I gives  the  organizations  consulted  with  their  locations  plotted  on  the  map,  Table  II.  Among  them,  some 
flight  test  centers  and  the  aircraft  and  helicopters  constructors  have  not  their  own  icing  ground  facilities.  They  use 
the  existing  ones,  carry  out  flight  tests  for  the  certification  of  aircraft  or  helicopters  and  sometimes  develop  icing  tech- 
niques. 

It  is  to  be  noticed  that  the  atmospheric  icing  conditions  are  studied  by  LDMA  at  Clermont-Ferrand  in  France 
with  cloud  characteristics  measurements  in  flight. 


Table  I — European  icing  facilities  and  techniques 


Country 

Town 

Organization,  acronym,  sponaor 

Austria 

Wian 

Fahrzaugvarsuchsanlaga  - Arsanal  (ground  vahiclas  tast  station) 

Franca 

Saint-Cloud 

Avions  Marcal  Dassault  - Brtguet  Aviation  AMD-BA 

Toulousa 

Cantra  d'Essais  Atronautiques  da  Toulousa  - CEAT 

Saclay 

Centra  d'Essais  da  Propulseurs-CEPr  (angina  units  tast  cantar) 

BrStlgny 

Cantra  d'Essais  an  Vol  - CEV  (Flight  Tast  Cantar) 

Bourgas 

Etablissamant  Tachniqua  das  Armamants  Tarrastras  - ETBS 

Clarmont-  Farrand 

Laboratoira  da  Dynamiqua  at  Microphysique  da  1' Atmosphere  * 

LDMA  (dynamic  A microphysics  of  atm.  lab.) 

Modana 

Off  lea  National  d'Etudas  at  da  Racharchas  Atrospatiales  (ONERA) 

Toulousa 

Aerospatiale  - SNIAS  - Aircraft  Division 

Marignanna 

" " * Haiicoptar  Division 

Oar  many 

KAIn 

Dautscna  Forschungs-und  Varsuchsanstalt  fur  Luft-und  Raumfahrt 

DFVLR 

Braman 

Varainigta  Flugtachnischa  Werke  Fokkar  - VFW 

Wolfsburg 

Kllma  Windkanal  Volkswagen  (climatic  W.T.) 

Italy 

Orbaasano 

Centro  Rleerche  Flet  (Flet  test  center) 

Norway 

Kjallar 

Defence  Research  Establishment 

S wad  an 

Bromma 

Flygtakniska  FArsAksanstalten-FFA  (aeronautical  research  institute) 

Unitad 

Artington 

> Lucas  Aerospace  Ltd 

Kingdom 

Burn  lay 

1 

Boscomba  Down 

Aeroplane  A Armament  Experimental  Establlshment-A  A AEE 

Farnborough 

National  Oeaturblne  Establishment  • NOTE 

Farnborough 

Royal  Aircraft  Establishment  - RAE 

Yeovil 

Westland  Helicopters  Ltd 

Facil.  Tw;hn. 

n test 
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Table  II  — Locations  of  European  facilities. 


1.2.  Icing  facilities  in  Europa 

Table  III  lists  the  European  tunnels  where 
icing  testing  is  possible,  with  some  details  about 
these  facilities.  In  fact,  these  facilities  comprise  : 

— conventional  wind  tunnels  with  special 
equipment  for  icing  testing, 

— altitude  cells  or  test  chambers,  where  the 
flow  is  induced  by  the  operation  of  the  tested 
engine  itself,  or  with  a free  jet,  which  are  in  fact 
open-jet  wind  tunnels.  For  example,  the  cell  n° 3 
west  of  NGTE  has  a chamber  7.6  m in  diameter 
equipped  with  a free  jet  2.4  m in  diameter. 
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Table  III  — Characteristics  of  European  icing  testing  facilities. 


Na 

tion 

Facility 

Dim  (m) 

Max.  spaad 
(m/s) 

Min. 

tamp. 

<°C) 

Icing  MVD 
(/Ml) 

Ramarks 

AU 

Ground  vahiclss  tact  station  - Vianna 

4.9  x 4.9 

32 

D 

Snow  a king  tasts 

CEAT  • Toulousa 

0 0.25 

0<  1 

250 

80 

20  -►  30 

• 

Two  tast  factions 
* Max.  flow  rata  : 10  g/s 

CEPr  • Saclay  - R2  Call 

R3  Call 

< o 1.1 

02 

150 

! s s 

i i i 
i 

j 

15  -►  30 

16  -*■  30 

0 -►  6 

0 10 

Etablistamant  Tachn.  Bourgas 

m 

22  -►  44 

(acc.  to  araa) 

-40 

- 

- 

Staam  injaction 

ONER  A Modana  tan  carrtar 

0 8** 

150 

—6  to 

-16* 

10  ->  20 

0.4  10 

* According  to  ambiar-t  cond. 

• * Icing  cloud  araa  : 4 m2 

sat  u rat  ad  tast  sac  t Ion 

— 

— . 

1 

Flat  rasa  arch  can  tar 

4x4 

45 

(Qj 

- . 

- 

No  Icing  tast  until  now 

— 

FRO 

Volkswagan  climatic  W.T. 

8-  37  m3 

66 

-36 

- 

- 

No  icing  tast  until  now 

U.K. 

A & AEE  - Bo  sco  mb  a Down 
Blowartunnal  a icing  facility 

0 1.2 

0 1.8 

140 

110 

-30 

-30 

30 

3 

LN2  injaction 

Artlngton  icing  tunnal 

0.30  x 0.30 

100 

-40 

7 

7 

Burn  lay  altftud.  tMt  facility 
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The  FFA  wind  tunnel  at  Bromma  is  not  really  an  icing  test  tunnel  : it  is  used  to  study  the  aerodynamics  of  ‘j  1 

profiles  including  ice  accretion  simulated  by  moulding.  In  fact,  all  the  conventional  wind  tunnels  can  be  used  simi-  , 

larly. 

j I 

i 

The  smallest  tunnels,  i.e.  the  tunnels  of  CEAT  at  Toulouse  or  Artington  icing  tunnel,  whose  testing  area  is  less  , 

than  1 m*.  are  used  to  test  small  aeronautical  equipments.  Beside  those  listed,  there  are  other  very  small  section  tunnels 
which  have  limited  application,  mainly  to  instrumentation  development,  such  as  that  of  the  Engineering  Physics 
Department  at  RAE  Farnborough  or  at  DFVLR. 

For  information  are  mentioned  the  tunnels  specialized  for  tests  on  automotive,  cars,  ground  vehicles  and  so 
forth.  Among  them,  the  Fiat  and  Volkswagen  climatic  tunnels  are  not  actually  used  for  other  tests.  The  tunnels  of 
Vienna  and  Bourges  offer  the  possibility  to  test  aircraft  parts  but  it  is  to  be  noticed  that  the  speed  is  very  low,  i.e. 
less  than  45  m/s. 

Some  wind  tunnels,  or  rather  some  cells  or  altitude  chambers  are  used  for  tests  on  engines  and  thus  are  in  the 
scope  of  the  PEP.  Nevertheless,  these  facilities  are  also  used  for  icing  tests  on  air  intakes  and  small  parts  of  aircraft 
or  helicopters  ; that  is  the  case  for  the  two  cells  of  the  CEPr  at  Saclay. 

The  A & AEE  tunnel  at  Boscombe  Down  is  an  open  jet  blower  tunnel.  The  two  smaller  nozzles  that  have  been 
recently  modified  for  icing  work  are  4 ft  and  6 ft  in  diameter.  Liquid  nitrogen  is  used  to  lower  the  temperature  of  the 
airstream.  A grid  of  water  spray  nozzles  provides  the  icing  cloud  ; the  area  in  front  of  the  tunnel  is  large  enough  to  * 

take  any  size  of  aircraft. 

The  Burnley  altitude  test  facility  has  been  used  for  icing  tests  since  the  late  1950  s.  Initial  anti-icing  tests  have 
been  carried  out  on  engine  intakes  for  single  engines,  but  for  more  advanced  tests  a larger  facility  is  needed. 

t: 

The  NGTE  Engine  Test  Facility  has  two  altitude  chambers  in  which  icing  tests  can  be  made.  The  majority  of 
tests  simulate  wet  icing  conditions,  but  some  work  has  been  done  to  investigate  icing  in  mixed  conditions,  that  is  an 
environment  containing  solid  ice  particles  and  supercooled  water  droplets.  Equipment  exists  in  cell  3 for  the  production 
of  solid  ice  particles  and  a similar  equipment  is  under  construction  for  use  in  cell  3 west  and  will  become  available 
during  1978.  Both  cells  have  been  used  to  clear  engines  for  operation  in  icing  conditions.  The  two  cells  can  also  be 
run  in  a free  jet  mode  so  that  icing  tests  can  be  made  on  complete  propulsion  units  (intake  and  engine)  or  on  full 
scale  air  frames, i.e.  a complete  helicopter  fuselage  (without  rotor)  with  engine  air  intakes.  To  support  the  icing  test 
programme,  NGTE  is  undertaking  work  on  instrumentation  so  that  the  cell  test  conditions  can  be  more  closely  defined 
and  controlled.  The  aim  is  to  provide  on-line  output  of  droplet  size  and  size  distribution,  liquid  water  content  and  ice/ 
water  ratio. 

The  largest  European  tunnel  for  icing  tests  is  the  SI  Modane  tunnel.  This  wind  tunnel  offers  a test  section  8 m 
in  diameter  in  which  the  whole  range  of  subsonic  velocities  can  be  obtained.  Thanks  to  climatic  conditions  and  cooling 
by  exchange  of  air  with  the  outside  atmosphere,  the  temperature  goes  down  often  in  winter  below  —5,  even  —10 
and  sometimes  — 15° C.  It  is  then  possible  to  create  in  front  of  the  models  an  artificial  icing  cloud  corresponding  to  the 
i standard  requirements.  To  this  end,  a 4 m1  pulverization  grid,  equipped  with  444  compressed  air  injectors,  distributes 

an  homogeneous  cloud.  The  air  and  water  supplies  of  that  grid  can  be  automatically  controlled  in  order  to  obtain  alter- 
nate icings,  while  keeping  the  droplet  diameter  constant  with  periods  from  3 to  several  tens  of  seconds.  The  gradual 
implementation  of  the  facility  allowed  a number  of  various  icing  tests.  The  first  ones  concerned  actual  aircraft  elements 
equipped  with  de-icers  : elements  of  wings  or  tails  (even  complete  tails),  fuselage  radomes.  In  some  cases,  and  in  parti- 
cular for  Concorde,  Airbus  or  for  helicopter  rotors,  tests  on  reduced  scale  models  were  mandatory.  By  associating  the 
correctly  adapted  spray  grid  and  the  system  for  testing  rotors,  icing  on  a model  rotor  up  to  4 m in  diameter  is  possible. 

Up  to  the  present,  the  SI  Modane  tunnel  can  be  used  for  icing  tests  only  during  the  winter,  when  adequate  temperature 
levels  can  be  reached.  To  palliate  these  drawbacks,  ONERA  is  studying  means  that  would  allow  artificial  cooling  of  the 
test  section.  The  realization  of  such  a project  would  make  it  possible  to  reduce  the  air  temperature  in  the  test  section 
by  15°C.  This  would  allow  one  to  do  away  with  local  atmospheric  conditions  for  a large  number  of  tests  and  to  carry 
out  the  very  low  temperature  tests  necessary  to  qualify  some  materials,  such  as  de-icing  and  anti-icing  devices,  in  the 
whole  required  range. 


1.3.  Bibliography 


Table  IV  gives  the  tested  models  in  e3ch  facility,  the  particular  techniques  developed  by  each  organization  and 
the  published  bibliography. 


Table  IV  — Information  gathered  on  icing  teat  techniques  in  Europe. 


Facility 

Tested  models/items 

GVTS  - Vienna 
(AU) 

Full  scale  : automotive-rallw. 
cars  ground  vehicles  • Aircraft 
parts 

A & AEE  - 
Boscombe  (UK) 
Blower  tunnel  & 
icing  facility 

Aircraft  wings  - Engines  & Air 
intakes  - Weapon  systems  • 
Helicopters. 

Artington  icing 
tunnel  (UK) 

Small  aeronautical  equipment 

Burnley  altitude 
test  facility  (UK) 

Air  intakes  - Small  helicopter 
engines  ■ De-icers 

NGTE  • Cell  3 
Farnbo rough  (UK) 

Air  intakes  - Aircraft  & Helicop- 
ter engines 

FFA  - 3.6  m W.T. 
Bromma  (S) 

2 D aerofoils,  c«  0.65  m 

CEAT  Toulouse 
(F) 

Small  aeronautical  equipment 

CEPr  Saclay  (FI 

Full  scale  : engines  - Air  intakes 
Nacelles 

ETBS  - Bourges 
(F) 

Large  ground  veh.  & equipment 
helicopter  engines  & sir  intakes 

ONE  RA  Modane 
(FI  SI  MA  W.T. 

Full  scale  aircraft  parts  : flap- 
ped wings  - Tail  surfaces  - Redo 
mas  - Reduced  scale  models  : 
half-models  (1/6  & 1/12).  Heli- 
copter rotors  (1/3) 

LDMA  - Aubltre 
(F) 

Flight  cloud  calibration  - Icing 
spray  grids  ground  callb. 

(ONERA,  CEAT.  CEPr) 

Remarks 

Bibliography 

Very  long  test  section 

Booklet  : Problems- Versuche-Forschung 

Free  air  tunnel  - LN2  injection 

Leaflets  (A  & AEE)  : 

- Canberra  WV  787  icing  tanker 

- Environmental  test  centre 

- Helicopter  iring  section 

- Blower  tunnel  & icing  facility 

(t  Booklet  : ) 

Extension  of  engine  test  range 
under  study 

Short  description  (Aug.  77)  Aircraft  ice 
protection  conference  1962  - D.  Napier 
engineering  report  DEV/R/ 1596/942. 
April  1973  - Lucas. 

Short  description  (Aug.  77) 

Icing  conditions  : wet  or  mixed 
(ice*  water) 

- Short  description  (Aug.  77) 

- Rotary  wing  icing  symposium  Edwards 

AFB.  June  1974. 

Simulated  ice  accretion 
shapes  - Moulding  technique 

FFA  Tech.  Notes  : AU  816  (1970)  - 
AU  902  (1972)  - AU  995  (1974)  - AU 
1169  (1977)  - Memo  90  (1973). 

Adaptable  TS  0.05  to  5 mJ 

Note  Technique  CEAT  15/N/77 

Duplication  of  flight  conditions 

"Fiches  sommairet"  on  R2  & R6  benchei 

Icing  tasting  possibla  (not  yat 
performad) 

Deter.  Note  : "Chambre  et  toufflerie  ell- 
matlque". 

Spray  grid  in  large  test  taction 
Force  meat.  Ice  accretion 
Moulding  - Red.  scale  testing  - 
Min.  temp,  according  to  am- 
bient 

TP  ONERA  No  1054  (1972) 

Rtf.  1 to  this  report. 

Appendix  to  this  report. 

Study  of  atmospheric  icing 
conditions  : droplet  dlam. 
spectra  & LWC  IKnollenberg 
& Rutkln  tensors) 

Scientific  reports  LDMA  : No  8 (1974) 

No  10  (1975),  No  11  (1976). 

Thetis  (3d  cycle)  1976. 
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IN  PERFORMING  ARTIFICIAL  ICING  TESTS  IN  WIND  TUNNELS 


In  this  paragraph,  the  particulars  of  the  artificial  icing  tests  in  wind  tunnels  are  emphasized.  The  difficulties  en- 
countered, easily  or  hardly  solved,  and  the  problems  which  remain  to  master  are  pointed  out,  relying  on  examples  taken 
from  SI  Modane  wind  tunnel  [1], 

2.1.  Facilities 

The  test  section  dimensions  allow  one  to  carry  out  tests  on  real  elements,  and  especially  to  determine  the  zones 
which  are  to  be  protected  against  the  ice,  and  to  study  the  efficiency  of  de-icing  or  anti-icing  devices.  Nevertheless, 
even  with  large  dimensions,  like  SI  Modane,  one  must  resort  to  scaled  down  models  for  the  helicopter  rotors  or  trans- 
port aircraft.  For  such  tests,  some  particular  similarity  rules  are  to  be  applied  (see  Section  2.3). 

The  interest  of  favourable  climatic  conditions  is  obvious  to  reach  sufficiently  low  temperatures  : so  SI  Modane 
wind  tunnel  benefits  of  alpine  conditions  at  medium  altitude.  Nevertheless,  to  extend  the  season  and  hours  suitable  for 
icing  tests  and  to  reach  lower  temperatures,  a supplement  of  cooling  is  necessary.  Two  techniques  can  be  considered  : 
cold  fluid  mixing  or  heat  exchanger. 

The  cold  fluid  mixing,  with  liquid  nitrogen  or  airstream  cooled  by  expansion,  does  not  bring  any  pressure  losses 
but  needs  a storage.  The  disadvantage  of  this  method  resides  in  the  ice  formation  at  the  cool  points,  which  destroys 
the  droplet  supercooling  and  thus  alters  the  icing  phenomena. 

This  disadvantage  disappears  with  cooling  by  exchanger  with  liquid  ammonia  or  freon.  Exchanger  icing  must  be 
avoided  by  dry  airstream  supply  and  an  open  circuit  operating  layout.  The  exchangers  generate  a pressure  loss  and 
consequently  require  extra  power.  This  method  is  rather  expensive  in  investment  cost  and  needs  a cooling  fluid  storage 
larger  than  the  liquid  nitrogen  arrangement. 

If  the  wind  tunnel  permits  conventional  tests,  the  icing  tests  will  take  advantage  of  the  existing  instrumentation. 
That  is  the  case  of  the  SI  Modane  tunnel,  used  for  icing  tests  but  not  primarily  intended  for  them.  All  the  test  para- 
meters and  their  variations  can  be  recorded.  The  photographic  motion  picture,  TV  and  magnetoscope  apparatus  are 
conventional. 

Trajectography  of  shedded  ice  parts  and  stroboscopy  for  rotor  tests  may  be  used  without  difficulties.  The  parti- 
cular icing  tests  method  consisting  in  the  duplication  of  the  ice  accretion  shape  by  casting  (Fig.  1)  is  easy  : this  method 
requires  a quick  access  to  the  model  after  icing,  to  avoid  alterations  of  the  ice  deposits. 
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Fig.  1 - Plaster  ice  cap  replicas  (other  presentation  of  fig.  9 of  the  appendix). 


2.2.  Artificial  icing  cloud  generation 

The  generation  of  the  artificial  icing  cloud 
must  satisfy  all  the  specifications  required  by  the 
standard  regulations  for  aircraft  certification  in 
icing  flight  conditions.  These  specifications,  esta- 
blished by  the  official  services,  define  the  icing 
cloud  standard  characteristics  from  statistical  me- 
teorological data  (Fig.  2).  The  respect  of  these 
specifications  implies  to  master  separately  three 
parameters  : 


continuous  max 


intermittent  max 


— the  liquid  water  content  (LWC) 

— the  mean  volumetric  diameter  (MVD) 

— the  spectral  distribution  by  diameter 
around  the  MVD. 


L WC  ( g/m3)t 
2 


Fig.  2 - Icing  standard  conditions. 


Compressed  air  injectors,  with  fixed  geometry,  can  be  controlled  with  only  two  parameters  : the  water 
and  air  flow  rates  fix  the  LWC  and  MVD.  To  obtain  a spectral  distribution  as  close  as  possible  to  the  standard 
specification  over  the  whole  range,  from  full  scale  to  1/12  scale,  the  injector  geometry  must  be  carefully  studied. 
For  this  purpose,  ONERA  developed  a new  type  of  spray  grid  and  atomizers.  Nevertheless,  some  problems  remain 
critical  : 


— small  MVD  (~  10  pm),  associated  with  large  LWC,  need  high  air  pressures,  then  high  expansion  rates, 
from  which  follows  a risk  to  freeze  prematurely  the  droplets.  This  inconvenience  can  be  overcome  by  heating 
the  air,  sometimes  up  to  70°C  ; 

— small  LWC,  combined  with  low  wind  speed,  lead  to  small  water  flow  rates  and  eventually  a poor  pulve- 
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Fig.  1 - Plaster  ice  cap  replicas  (other  presentation  of  fig.  9 of  the  appendix). 


2.2.  Artificial  icing  cloud  generation 

The  generation  of  the  artificial  icing  cloud 
must  satisfy  all  the  specifications  required  by  the 
standard  regulations  for  aircraft  certification  in 
icing  flight  conditions.  These  specifications,  esta- 
blished by  the  official  services,  define  the  icing 
cloud  standard  characteristics  from  statistical  me- 
teorological data  (Fig.  2).  The  respect  of  these 
specifications  implies  to  master  separately  three 
parameters  : 

- the  liquid  water  content  (LWC) 

- the  mean  volumetric  diameter  (MVD) 

- the  spectral  distribution  by  diameter 

around  the  MVD.  Fig.  2 — Icing  standard  conditions. 

Compressed  air  injectors,  with  fixed  geometry,  can  be  controlled  with  only  two  parameters  : the  water 
and  air  flow  rates  fix  the  LWC  and  MVD.  To  obtain  a spectral  distribution  as  close  as  possible  to  the  standard 
specification  over  the  whole  range,  from  full  scale  to  1/12  scale,  the  injector  geometry  must  be  carefully  studied. 
For  this  purpose,  ONERA  developed  a new  type  of  spray  grid  and  atomizers.  Nevertheless,  some  problems  remain 
critical  : 

- small  MVD  (~  10  pm),  associated  with  large  LWC,  need  high  air  pressures,  then  high  expansion  rates, 
from  which  follows  a risk  to  freeze  prematurely  the  droplets.  This  inconvenience  can  be  overcome  by  heating 
the  air,  sometimes  up  to  70°C  ; 

- small  LWC,  combined  with  low  wind  speed,  lead  to  small  water  flow  rates  and  eventually  a poor  pulve- 
rization. 
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A very  difficult  problem  is  to  get  a good  uniformity  of  the  distribution  in  the  useful  icing  cloud  section. 
That  involves  : 

- a good  atomization  for  each  injector,  avoiding  the  local  segregation  of  the  droplets  in  the  injector  field  ; 
this  entails  a thorough  maintenance  and  repeated  injector  calibrations ; 

- an  accurate  identity  of  the  injectors,  which  requires  a close  fabrication  tolerance  (~  10  pm)  and  an  accu- 
rate mounting  of  the  injector  pieces  (at  SI  Modane  the  spray  grid  contains  444  injectors) ; 

— the  control  of  water  supply,  separately  for  each  horizontal  grid  bar  ; 

— a good  design  of  the  grid  bar  which  holds  and  supplies  the  injectors  from  both  aerodynamic  and  mecha- 
nical points  of  view,  in  order  to  reduce  the  grid  drag  and  avoid  the  bar  bending  vibrations. 

To  avoid  the  droplet  evaporation  between  grid  and  model  or  along  the  model,  which  is  mainly  due  to  fine 
pulverization,  an  auxiliary  installation  makes  it  possible  to  reach  saturation  in  the  test  section  by  increasing  the 
hygrometric  factor  from  60  to  nearly  100  %.  At  SI  Modane,  a warm  water  injection,  downstream  of  the  model, 
ensures  the  constancy  of  icing  cloud  characteristics  along  the  test  section. 

To  cope  with  the  various  shapes  and  locations  of  the  models  in  the  test  section,  the  spray  grid,  when  it 
does  not  cover  the  whole  test  section  area,  has  to  be  adapted  in  location  and  configuration.  At  Modane,  the 
spray  grid  can  cover  a 4 m2  area  with  various  arrangements  of  the  grid  bars  : figure  4 of  appendix  shows  an 

example  for  a half  model  of  Concorde  mounted  on  the  test  section  floor. 

Another  example  (fig.  6 of  appendix)  shows  the  spray  grid  array  for  rotor  models.  The  aerodynamic 
field  is  surveyed  by  smoke  visualizations.  The  references  must  be  taken  behind  the  grid  to  take  into  account 

the  injector  air  mass  flow.  Comparative  tests,  with  and  without  grid,  guarantee  the  validity  of  the  results. 

To  obtain  the  alternance  of  the  continuous  and  intermittent  icing  conditions,  between  them  or  with  cleat 
sky,  specified  by  certification  icing  standards,  the  liquid  water  content  must  vary  by  sharp  steps  (less  than  1 sec), 
while  keeping  constant  the  mean  volumetric  diameter.  This  problem  has  been  solved  by  doubling  the  whole 
array  of  water  and  compressed  air  supply,  presetting  the  desired  rate  values  for  alternate  conditions  and  switching 
from  one  to  the  other  condition  and  reversely  by  fully  automatic  control  of  quick  acting  valves,  taking  into 
account  the  response  time  of  the  piping. 

The  calibration  of  the  icing  cloud  is  an  uneasy  problem  for  all  the  icing  tests,  because  it  should  be  made 
"in  situ"  near  the  model. 

The  total  water  flow  rate,  divided  by  the  airspeed  and  the  icing  cloud  area,  gives  the  liquid  water  content. 
This  assumes  that  no  evaporation,  no  coalescence,  small  marginal  effects  and  a good  homogeneity  are  ensured. 

Local  measurements  are  often  less  reliable. 

The  mean  volumetric  diameter  and  the  droplet  spectra  are  determined  by  photometric  analysis  or  mecha- 
nical catching  of  samples. 

The  catching  on  oiled  plate,  developed  by  Napier,  is  a slow,  cumbersome  and  questionable  process  because 
of  the  low  sampling  rate  and  the  variation  of  the  specimens  between  sampling  and  analysis.  Moreover,  this 
process  underrates  the  smaller  droplets,  if  no  catching  factor  correction  is  made. 

The  catching  on  the  softened  surface  of  a film,  along  the  Bertin  method,  gives  a stable  stamp  ; the 
difficult  calibration  of  the  stamp/droplet  diameter  ratio  was  overcome  by  an  elegant  method,  using  salty  solution 
of  known  concentration. 

During  the  analysis  in  situ,  it  is  difficult  to  avoid  alterations  of  the  cloud,  due  to  the  probe  itself  and  parasitic 
condensations.  So  far,  the  agreement  between  simultaneous  calibration,  by  the  several  known  methods,  is  not  quite 
satisfactory. 
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2.3.  Scaled  down  model  tests 
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As,  in  some  cases,  tests  on  reduced  scale  models  are  compulsory,  similarity  rules  are  then  to  be  applied  (see 
appendix).  In  fact,  one  has  to  call  up  on  a restricted  similitude  which  had  been  established  for  those  cases,  by  resear- 
chers from  Douglas,  Lockheed,  NAE, and  later  of  BAC  and  ONERA.  This  similitude  respects  the  catching  rate  and  the 
surface  heat  exchanges  through  the  parameters  : LWC,  MVD,  air  speed  and  icing  duration.  Moreover,  for  helicopter 
rotors,  after  a certain  icing  time,  the  importance  of  ice  deposit  stabilizes,  due  to  the  ice  shedding  under  the  effect  of 
centrifugal  force.  For  tests  on  reduced  scale  rotors,  it  is  necessary  to  simulate  correctly  this  phenomenon,  so  that  the 
stresses  on  the  ice  accretion  under  the  centrifugal  force  be  the  same  as  in  flight  : this  entails  the  identity  of  tip  speeds, 
thus  of  forward  speeds  [2].  The  respect  of  similitude  rules  is  then  obtained  by  acting  on  the  droplet  size,  the  liquid 
water  content  and  the  rotor  r.p.m.  During  all  the  time  when  the  rotor  is  maintained  in  icing  condition,  the  rotor  r.p.m. 
is  kept  constant  by  acting  on  the  driving  motors  output. 


Such  a restricted  similitude  allows  the  overall  study  on  models  for  the  captation,  accretion  and  ice  shedding,  and 
for  the  determination  of  the  zones  which  are  to  be  protected  against  ice.  The  Reynolds  number  and  sometimes  the 
Mach  number  are  not  respected.  Nevertheless  the  similitude  ensures  reasonably  well  the  thermodynamic  exchanges  go 
verning  the  formation  of  ice  deposits,  because  the  flight  speeds  are  rather  low  and  the  main  phenomena  occur  around 
the  leading  edge  where  the  flow  pattern  is  less  affected  by  viscosity  effects.  The  altitude  and  the  thermodynamics  are 
only  approximately  respected.  For  the  rotor  blades,  the  aeroelastic  characteristics  were  not  simulated  until  now. 


fig.  3 — Accretion  on  cylinders 
at  full  and  1/6  scale. 
(Combination  of  fig.  3 and  12 
of  the  appendix). 


According  to  these  considerations,  ONERA  has  assessed  the  validity  of  such  icing  tests  on  scaled  down  models 
with  an  experimental  control  of  the  similarity  rules  by  applying  it  in  the  wind  tunnel  to  test  pieces  respectively  at 
full  and  reduced  scales  : as  an  example,  figure  3 shows  the  shapes  of  accretion  obtained  on  two  cylinders  at  scales  1 
and  1/6.  The  accretions  are  geometrically  quite  similar  with  respect  to  the  similarity  rules.  Naturally,  if  the  two  cylin 
ders  are  subjected  simultaneously  to  the  same  icing  conditions,  for  instance  full  scale  conditions,  the  results  obtained 
on  the  1/6  model  are  quite  insignifiant. 

The  assessment  with  delta  wings  with  sharp  leading  edges  at  incidence,  and  on  conventional  profiles,  gives  good 
agreement  on  both  the  extent  and  the  shape  of  the  accretion. 


-i 


6-10 


I 

The  flight/wind  tunnel  comparisons  offer  an  essential  complement  to  validate  this  method.  Unfortunately,  few 
flight  results  can  be  compared  to  wind  tunnel  tests  carried  out  previously  (naturally,  it  would  be  easier  to  reproduce  ( 

in  tunnel  the  flight  tests  conditions). 

I 

> 

According  to  wind  tunnel  tests  on  a 1/6  half  model  of  Concorde  at  SI  Modane,  the  deposits  on  the  upper 
surface  were  almost  non  existent,  whereas  on  the  leading  edge  the  ice  deposits  were  located  at  the  lower  surface  with 
a stratified  shape  generating  natural  tear-offs  of  ice  blocks.  Such  concretions  were  found  again  in  flight,  and  the  tra- 
jectories of  the  naturally  shedded  ice  parts  were  similar  to  those  seen  in  wind  tunnel. 

I 

Icing  tests  on  a 1/12  model  of  Airbus  led  to  the  determination  of  the  captation  zones  on  the  wing  and  tail 
leading  edges  and  on  the  nacelles  lips.  In  both  tunnel  and  flight,  the  stratified  ice  deposit  shape,  so-called  "lobster 
tail"  [3j,appears  figure  4 : it  is  characteristic  of  the  leading  edge  icing  on  a swept  back  wing. 


Fig.  4 — Lobster  tail  ice  shape  of  A300  in  flight  and  in  tunnel  tests. 

(other  presentation  of  fig.  20  of  the  appendix). 

In  the  particular  case  of  rotor  blades  during  an  initial  period  lasting  5 to  6 minutes  on  the  model  of  a 3-blade, 
1/3  scale  rotor,  the  rotor  coefficients  are  rapidly  and  significantly  altered  by  ice  accretion  (fig.  23  of  appendix).  After 
that  time,  the  ice  starts  breaking  up  locally  under  the  action  of  the  centrifugal  force  and  as  an  average  there  is  equality 
between  the  ice  volume  deposition  by  accretion  and  that  evacuated  by  centrifugation.  Rotor  coefficients  then  vary 
slightly  around  mean  values.  The  torque  increase  necessary  to  maintain  the  rotation  rate  constant  is  of  the  same  order 
of  magnitude  as  that  mentioned  by  Lake  and  Bradey  [4]  in  flight. 

As  the  tunnel  test  temperature  was  not  very  low  (— 4°C),  a fraction  of  water  not  congealed  at  impact  runs  along 
the  blade  surface  under  the  effect  of  centrifugal  forces  (fig.  22  of  appendix).  Beyond  a relative  radius  of  0.5,  the  sur- 
faces are  clear  of  ice,  probably  due  to  kinetic  heating.  These  observations  agree  with  those  made  in  flight  by  Stallabrats 
[5]  for  temperatures  between  -1  and  -4°C. 

This  set  of  comparisons,  though  made  on  non-similar  rotors,  shows  agreements  that  augur  well  for  the  represen- 
tativity  of  phenomena  studied  in  the  wind  tunnel  on  rotary  wings  over  the  whole  flight  domain. 
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SUMMARY 

Since  1962,  the  large  SI  wind  tunnel  at  Modane  has  been  used  to  perform  icing  tests  on  aircraft  parts  (wing,  tail 
surfaces,  radome  ...)  at  full  scale. 

During  the  last  years,  studies  and  experimental  verifications  have  been  done  to  promote  a technique  of  icing  tests 
on  down-scaled  fixed  wings.  The  equations  defining  the  non  dimensional  ratios,  which  have  to  be  satisfied  to  obtain  an 
acceptable  similitude  of  water  catching  and  development  of  ice  accretion,  are  given. 

Lately  this  down-scaled  icing  test  technique  has  been  extended  to  rotary  wings.  The  evolutions  of  the  rotor  aero- 
dynamic coefficients  in  icing  conditions  and  the  ice  distribution  on  the  blades  are  very  similar  to  those  observed  in 
actual  flight. 


♦ 

TECHNIQUES  ET  MOYENS  D'ESSAIS  DE  GIVRAGE  UTILISES  A L'ONERA  (CENTRE  DE  MODANE  I 
RESUME 

Depuis  1962,  I'ONERA  utilise  la  grande  soufflerie  SI  du  Centre  de  Modane  pour  realiser  des  essais  de  givrage  sur 
des  elements  d'avions  k echelles  grandeur  (ailes,  empennages,  raddmes). 

Au  court  des  derniires  annles,  des  etudes  et  des  verifications  experimental  ont  6t6  r£alis6es  dans  le  cadre  de  la 
mise  au  point  d'une  technique  de  givrage  & tchelle  rtduite  sur  des  voilures  fixes.  Les  equations  definissant  les  rapports  k 
respecter  pour  obtenir  une  similitude  acceptable  de  la  captation  et  du  developpement  des  concretions  de  glace  sont  rap 
peMes. 

Dernierement,  cette  technique  de  givrage  k echelle  reduite  a ete  etendue  aux  voilures  tournantes.  Les  resultats 
obtenus  dans  ce  domaine,  tant  du  point  de  vue  des  evolutions  des  coefficients  aerodynamiques  des  rotors  en  condition 
givrante  que  de  la  repartion  de  la  glace  sur  les  pales,  sont  tret  voitint  des  observations  faites  en  vol  reel. 

♦ 

VE RE ISUNGSTECHNI K UND  - VERSUCHSEINRICHTUNG  DER  ONERA 
ZUSAMMENFASSUNG 

Seit  1962  benutzt  die  ONERA  den  grossen  Windkanal  SI  des  Modaner  Versuchszentrums  zur  Durchfiihrung  von 
Vereitungtversuchen  an  Flugzeugbauteilen  wahrer  Gross*  (Flugel,  Leitwerk,  Radom). 

Im  Rahmen  der  Entwicklung  einer  Vereisungstechnik  im  verkleinerten  Mass-stab  wurden  im  Verlauf  der  vergangenen 
Jahre  experimented*  Versuche  und  Nachprufungen  am  TragflOgel  durchgefiihrt. 

Es  warden  die  Gleichungen  angegeben,  welch*  die  zu  beachtenden  Grotaenverhaltnisse  definieren,  urn  eine  annehmbare 
Ahnlichkeit  bei  der  Entstehung  und  Ausbildung  der  Eitanlagarung  zu  gawihrlaitten. 

In  ietztar  Zeit  wurde  diet*  Vereisungstechnik  auf  Drehflugel  erweitert.  Die  hierbei  erzielten  Ergebnitse  sowohl 
hinsichtlich  des  Verlauft  der  aerodynamischen  Rotorbeiwerte  alt  auch  bezuglich  der  Eisentwicklung  entlang  der  Blatter 
stfmmen  mit  den  Im  Flug  gemachten  Beobachtungen  gut  uberein. 
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NOTATIONS  ( 


a Median  volume  diameter  of  droplets  (pm) 

b Heat  relative  factor  [3] 

c Blade  chord 

CD  Actual  drag  coefficient  of  a droplet 

Cf  Rotor  torque  coefficient,  also  rotor  power  coefficient.  If  Cj  is  the  torque  and  l/V  the  shaft  power  : 

C1  = 100  Cj/lpNcR'U2  = 100  W'-LpNcR'uU2 

k Model  scale 

LWC  Liquid  water  content  in  clouds  (g/m3) 

Mg^  Mach  number  at  advancing  blade  tip  : Mg^  = {V*  U)  / sound  velocity  in  the  test  section 
n Icing  fraction  [3] 

N Number  of  blades.  Here  /V-  3 

Pg  Ambient  atmospheric  pressure  (Pa) 

PQ  Static  pressure  in  the  test  section  (Pa) 

PPge  Partial  water  vapour  in  the  atmosphere  (Pa) 

R Rotor  radius  (m) 

R Relative  radius  (distance  to  axis,  relative  to  rotor  radius) 

Rv  Reynolds  number  relative  to  droplet  diameter  and  droplet  velocity 
f<»,  T.  Air  static  temperature  at  infinity  upstream  (°C,  K) 

er  a 

U Rotor  tip  velocity.  U-  u>R  (m/s) 

V Rotor  forward  velocity  (m/s) 

V x Velocity  at  infinity  upstream  (m/s) 

Rotor  propulsive  force  coefficient  (parallel  to  wind).  If  F is  the  propulsive  force  : X = 100  x/^pNcRU1 

Rotor  lift  coefficient  (vertical).  If  Z is  the  lift  : Z = 100  z/—pNcRU 2 

2 ... 
a q Pitch  angle,  or  inclination  of  the  plane  described  by  the  blade  tips  on  the  horizontal  plane  ; positive  for  pitch 

up  (degrees). 

Uq  Inclination  of  the  rotor  shaft  on  the  vertical  ; positive  for  pitch  up  (degrees) 
d c Cyclic  pitch  (degrees) 

0o  Collective  pitch  (degrees) 

A Advance  ratio.  A = V/U 
p Air  density  (kg/m3 ) 

o Shear  stress  in  ice  (Pa) 

r Icing  time  (s) 

« Rotor  angular  velocity  (rad/s) 

Subscripts 

f Actual  aircraft  or  full  scale  model 

k Ratio  uf  the  values  of  a parameter  lor  scale  k and  full  scale  (ex 
m Model 
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INTRODUCTION 


There  exist  different  means  for  defining  the  shapes  and  extents  of  ice  deposits  on  aircraft  frames  and  on  heli- 
copter blades,  and  for  qualifying  de-icing  and  anti-icing  devices  : 

- flight  tests  in  natural  icing  conditions  ; 

— flight  tests  in  artificial  icing  conditions  : 

, of  actual  parts  mounted  on  a carrier  aircraft,  behind  a water  pulverizing  grid, 

. of  the  aircraft  itself,  flying  in  an  artificial  cloud  behind  a tanker  aircraft  ; 

- ground  tests  in  stationary  operation  (helicopters)  ; 

— wind  tunnel  tests  : 

. of  actual  parts  or  at  full  scale, 

. of  full  or  half  models  at  reduced  scale. 


The  study  in  artificial  cloud  implying  the  creation  of  clouds  of  the  same  nature  as  those  encountered  in  flight, 
it  is  mandatory  to  know  the  character, sties  of  the  latter. 

Statistical  meteorological  data  enabled  the  Official  Services  to  define  the  dimensional  characteristics  of  icing 
clouds  liable  to  be  encountered  by  aircraft  in  flight,  as  well  as  the  technical  conditions  for  granting  the  icing  item  its 
airworthiness  certification. 

In  artificial  icing  test  procedures,  among  them  wind  tunnel  tests,  two  cases  should  be  considered  : maximum 
continuous  icing,  and  maximum  intermittent  icing. 


The  atmospheric  conditions  corresponding  to  maximum  icing  are  defined  for  water  droplets  of  20  pm  median 
volume  diameter  in  the  following  conditions  : 


Atmospheric  temperature  (°C) 

Liquid  water  content  (g/m3) 

Continuous 

Intermittent 

0 

0.8 

2.5 

- 10 

0.6 

2.2 

- 20 

0.3 

1.7 

- 30 

0.2 

1.0 

- 40 

- 

0.2 

The  temperature  and  pressure-altitude  ranges  where  these  types  of  cloud  may  be  encountered  are  defined  on 
figure  1 . 


Fig.  1 - Maximum  icing  condition!,  a)  Continuous  b)  intermittent. 
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The  extent  of  these  types  of  cloud  is  considered  as  unlimited  ; however,  for  demonstration  tests  it  is  taken 
equal  to  30  minutes  for  cruise  or  stand-by  flight  configurations,  and  to  the  actual  flight  time  for  climb  or  descent 
flight  configurations. 

For  intermittent  icing,  the  extent  of  the  cloud  is  supposed  to  be  5 km  ; however,  for  demonstration  tests  the 
corresponding  cloud  formations  are  made  of  5-km -wide  clouds,  separated  either  by  clear  sky  intervals  or  by  clouds  of 
maximum  continuous  type,  defined  as  above. 


1 - STUDY  OF  SIMILITUDE  LAWS 

The  theoretical  study  of  the  laws  governing  icing  in  similitude  conditions  requires  the  analysis  of  various  physical 
phenomena.  The  first  one  is  that  of  the  mechanics  of  water  droplet  pick-up,  which  involves  the  investigation  of  the 
trajectories  of  these  droplets  in  the  aerodynamic  field  of  the  flow  around  the  icing  obstacle.  A study  of  this  problem 
has  been  made  by  Langmuir  and  Blodgett,  of  General  Electric  Company  [3]. 

The  second  of  these  phenomena  is  that  of  the  distribution  of  water,  then  of  ice,  on  the  profiles  ; it  involves  the 
water  content  of  the  cloud,  the  icing  time,  and  also  the  convection  heat  transfer  coefficient,  defined  by  Gelder  and 
Lewis  in  ref.  [ 2). 

The  third  phenomenon,  perhaps  the  most  difficult  to  treat,  is  that  of  the  thermal  balance  of  a wall  in  icing  condi- 
tions. A very  detailed  study  of  this  problem  has  been  performed  by  Messinger,  of  Lockheed  Aircraft  Corp.  [3],  in 
which  he  brings  to  light  the  concepts  of  "heat  relative  factor,  b"  and  of  "icing  fraction,  n" . 

This  various  studies  were  put  together,  with  a view  to  define  the  relations  between  the  parameters  governing  wind 
tunnel  similitude  tests,  by  Mauger,  Englar  and  Reaser,  of  Cbuglas  Aircraft  Comp.  [4]  and  by  Brun , of  Paris  University 

15]- 

For  industrial  tests  in  similitude  conditions,  and  for  similitude  icing  tests  on  simple  models  intended  to  check 
the  theory,  the  parameters  defining  icing  clouds  have  been  calculated  from  the  works  by  Googan,  Jackson  and 
Hubbold,  of  British  Aircraft  Corp.  [6,  7].  These  works  made  it  possible  to  carry  out  the  transposition  from  the  air- 
craft flight  range  to  that  simulated  in  the  wind  tunnel,  within  the  framework  of  icing  restricted  similitude.  This  simi- 
litude is  restricted,  as  it  is  not  possible  to  respect  together  the  Mach  and  Reynolds  number  conditions  for  the  aerody- 
namic flow  around  the  obstacle  and  the  Reynolds  number  condition  for  the  droplets.  But  it  happens  that,  in  the 
flight  conditions  where  icing  might  the  dangerous,  the  aerodynamic  field  in  the  ice  pick-up  zones  is  little  sensitive  to 
these  two  parameters,  as  the  compressibility  and  viscosity  phenomena  in  the  boundary  layer  remain  inimportant 
around  the  leading  edges  at  low  or  moderate  speed. 

The  analysis  of  the  aerodynamic  and  inertial  forces  acting  on  the  water  droplets  close  to  the  obstacle  leads  to 
eq.  (1)  relating  the  various  parameters  of  the  droplet  trajectory  similitude  : 


In  eq.  (1),  x is  chosen  so  as  the  approximate  formula  : 

C0.Rv/24  = K(RV)X 

represents  at  best  the  evolution  of  the  actual  droplet  drag  as  a function  of  the  Reynolds  number  Rv  related  to  its  velo- 
city relative  to  the  flow,  within  the  range  of  values  of  Ry  where  this  relative  velocity  intervenes  determinantly  on  the 


fc 
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trajectory  deflexion  from  the  streamlines,  at  full  scale  and  at  reduced  scale  in  wind  tunnel  tests.  Figure  2 shows  that 
for  a droplet  Reynolds  number  range  between  6 and  120,  covering  both  the  flight  domain  and  the  1/6  and  1/12 
scale  domain,  the  value  of  x represented  by  the  slope  of  the  straight  line  passing  at  best  through  the  points,  can  be 
taken  as  0.39. 

10 

6 

< 

2 

1 

4 6810  20  40  60  80  tIO  200  300 

Fig.  2 - Ratio  between  the  actual  drag  of  a droplet  and  that  given  by  the  Stokes  law 

at  a function  of  Reynolds  number. 
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1.1.  Icing  similitude  on  aircraft  parts  t 

Icing  similitude  between  aircraft  and  model  is  ensured  if  a number  of  conditions  is  respected,  concerning  the  aero- 
dynamic field,  the  thermodynamic  exchanges  between  air,  water  and  ice,  the  droplet  trajectories,  the  water  content  in 
icing  clouds  and  the  crossing  duration  of  these  clouds.  » 

1 

For  an  aircraft  flying  in  given  conditions,  the  flight  is  characterized  by  the  velocity  (V,^),  the  angle  of  attack 
taf)  and  the  flight  altitude  (Hf).  At  this  altitude  corresponds  a static  pressure  ( Pgf)  and  an  ambient  temperature  (tfl  ). 

The  icing  cloud  is  characterized  by  its  w'ter  content  (LWC)f,  the  droplet  median  volume  diameter  (a^l  and  the  water 
vapor  partial  pressure  (PP^)-  The  ice  deposit  on  the  wing  depends  on  the  above  parameters  and  the  cloud  crossing 
duration  (if). 

For  a model  of  this  aircraft  at  scale  k,  tested  in  the  wind  tunnel,  the  equivalent  data  (subscript  m)  have  to  be 
determined  in  relation  to  those  if  the  aircraft  in  flight  or  at  full  scale  (subscript  f). 

1. 1. 1.  Aerodynamic  similitude 

The  aerodynamic  similitude  is  first  ensured  by  the  similitude  of  shapes  and  the  equality  of  angles  of  attack  : 

a = otf. 
m t 

The  flight  velocities  in  icing  conditions  being  usually  low,  the  aerodynamic  field  does  not  imply  transonic  phe 
nomena,  and  is  not  modified  by  a velocity  variation.  The  Reynolds  number  effect  being  usually  negligible,  there  is 
not  necessarily  an  equality  between  flight  and  wind  tunnel  velocities.  On  the  other  hand,  the  aerodynamic  coeffi- 
cients are  the  same  in  both  cases. 

1. 1.2.  Thermodynamic  similitude 

The  parameters  characterizing  thermodynamic  exchanges  between  air,  water  and  ice  satisfy  the  relation  (5).  The 
four  parameters  P ^ f PPX  and  Vx  determine  the  conditions  in  which  thermodynamic  exchanges  take  place,  so  that 
eq.  (5)  is  reduced  to  a homographic  relation  between  n and  b : 


(6) 


Aj * B f b 
b 


This  relation  fixes  the  nature  of  the  thermodynamic  exchanges.  For  the  similarity  to  be  ensured,  the  set  of 
values  for  parameters  P tg,  PP^  and  should  lead  to  the  same  relation  (6),  for  the  flight  and  the  wind  tunnel,  so 
that  : 


(7) 


Af*  Bf  b A * B b 
T t _ m m 

b b 


(8) 


or 


A — A t 
m t 


B = Be 
m f 


The  values  of  the  four  parameters  are  known  in  flight,  so  that  Af  and  Bf  are  determined.  In  the  wind  tunnel, 
the  stagnation  pressure  and  temperature  are  known,  and  there  exists  a relation  between  V',,  tg  and  Pg  . Lastly, 
there  is  a relation  between  the  water  vapour  partial  pressure  (PPxm)  and  the  static  temperature  itgm).  So  there  exist 
four  relations  between  Pa^y.  tam-  ^<xm  an<^  these  are  all  determined  if  the  thermodynar  x simifitude  has  to  be 
achieved. 
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1. 1.3.  Similitude  of  droplet  trajectories 

Aircraft  and  model  being  geometrically  similar,  at  scale  k,  we  must  make  sure  that  the  droplet  trajectories  are 
also  similar  in  flight  and  in  the  wind  tunnel  '(he  aerodynamic  field  being  the  same,  the  trajectories  depend  only  on 
the  droplet  size.  The  trajectory  similitude  is  ensured  if  the  median  volume  diameters  in  flight  (a^  and  in  the  wind 
tunnel  (am)  satisfy  eq.  ( 2). 


(2') 


k P: 


0.39 


1.61  _ 


<>k 


V , 0-6 1 
°°k 


In  this  relation,  all  the  terms  on  the  right  hand  side  are  fixed  by  the  previous  conditions  and,  af  being  known, 
am  is  determined.  In  fact,  the  experimenter  will  choose  the  parameters  so  as  the  similitudes  be  respected  and  am 
might  be  actually  obtained  with  the  existing  pulverization  device. 

Equation  (2)  must  be  respected  for  the  droplet  trajectories  to  be  correct  and  the  ice  pick-up  on  the  wing  be 
realistic. 

1. 1.4.  Similitude  of  ice  deposits 

The  importance  of  the  ice  deposits  depends  on  the  quantity  of  liquid  water  in  the  icing  clouds  (LWC)  and  the 
cloud  crossing  duration  <T). The  analysis  of  the  various  parameters  acting  on  the  volume  and  distribution  of  ice  on  the 
profiles  leads  to  the  relations  of  similitude  defined  in  eq.  (3)  and  (4). 


The  set  of  relations  (5),  (2),  (3)  and  (4)  allows  the  solution  of  icing  problems  at  reduced  scale  on  models  of 
fixed  wing  vehicles. 

1.2.  Icing  similitude  for  rotary  wing  aircraft 

Their  laws  are  defined  from  these  valid  for  fixed  wing  aircraft 

1.2. 1.  Aerodynamic  similitude 

It  is  ensured  of  the  aerodynamic  field  for  each  blade  section  is  the  same  during  a revolution  in  flight  and  in 
wind  tunnel,  on  both  the  full  scale  and  the  scaled-down  model.  To  this  end,  the  shaft  inclination  tag),  the  collective 
pitch  (6q)  and  the  cyclic  pitch  ( 6C ) are  adjusted  so  that  the  disc  inclination  («g)  be  that  of  the  flight  and  that  the 
geometric  incidences  be  realistic.  The  rotation  speed  (to),  tip  speed  (U)  and  advancing  speed  (V)  are  adjusted  so  that 
the  velocity  composition  be  the  same  in  flight  and  wind  tunnel  on  any  point  of  the  disc  ; for  this,  the  advance  ratio 
(A  = V/U)  must  be  the  same.  The  Mach  number  of  the  advancing  blade  tip  (Mgb)has  to  be  maintained  equal  only  if 
compressibility  phenomena  are  significant.  This  was  not  the  case  for  the  tests  at  SIMA,  but  as  equality  of  speeds  was 
required  for  respecting  equality  of  stresses  in  ice  (see  section  1.2.4),  equality  was  respected  for  V,  U,  and  A. 
Consequently,  the  angular  velocities  were  in  the  inverse  ratio  of  k : 


Um  = um  Rm  = Uf  * uf  Rf 
(9) 

uk  = = 1/Rk  = 1/k' 

The  correct  adjustement  of  these  various  parameters  is  verified  globally  if  the  lift  and  torque  rotor  coefficients 
are  the  same  in  the  wind  tunnel  as  in  flight. 
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1.2.2.  Thermodynamic  similitude 

Equation  (5),  valid  for  fixed  wings,  applies  here  in  the  same  way,  the  parameters,  Pg,  t ^ PP0 „ n and  b having 
the  same  definition,  and  being  the  instantaneous  velocity  on  the  profile  section  considered,  resulting  from  the 
composition  of  the  rotating  and  advancing  linear  velocities. 

Thus,  the  thermodynamic  similitude  consists  in  making  to  every  flight  velocity  correspond  a wind  tunnel 
model  velocity  Vaom  which  is  a priori  different.  There  is  there  a contradiction  with  the  aerodynamic  similitude,  which 
had  been  chosen  so  as  to  impose  equal  local  velocities  in  wind  tunnel  and  in  flight.  Practically,  this  difficulty  was 
solved  by  keeping  to  a velocity  calculated  to  satisfy  the  aerodynamic  similitude.  The  error  resulting  in  the  thermody- 
namic similitude  in  the  tests  at  S1-MA  is  about  0.5° C,  a value  of  the  same  order  of  magnitude  as  that  of  the  measu- 
ring precision  of  temperature  (f.)  during  the  tests  ; for  that  reason,  it  has  been  considered  as  negligible. 
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1.2.3.  Water  droplet  trajectory  similitude 

The  aerodynamic  similitude  ensuring  the  equality  of  the  aerodynamic  fields  on  similar  profiles,  the  trajectory 
similitude  is  ensured  in  the  same  condition  as  for  fixed  wings,  and  the  same  eq.  (2)  applies.  The  values  of  velocities 
Vgg  \iary  with  the  profile  section  and  the  blade  position  during  a turn  ; but  the  equality  of  the  advance  ratio  implies 
that  the  velocity  ratio  is  the  same  at  every  point.  In  the  particular  case  of  the  S1-MA  tests,  the  velocities  being 
equal,  the  ratio  is  1. 


1.2.4.  Ice  deposit  similitude 

In  the  same  way,  the  same  laws  apply  as  regards  liquid  water  con'ents  and  cloud  crossing  duration.  But  moreover 
in  the  case  of  rotary  wings  account  must  be  taken  of  the  centrifugal  forces  due  to  blade  revolution,  that  provokes  shedd 
ing  of  ice  deposits.  For  these  sheddings  to  be  the  same  in  flight  and  in  the  wind  tunnel,  it  is  necessary  first  that  shear 
stresses  in  ice  be  the  same  in  both  cases.  But,  all  the  parameters  being  fixed,  the  ratio  existing  between  these 
stresses  is  an  established  fact.  Assuming  the  ice  density  to  be  the  same  in  flight  and  in  the  wind  tunnel  : 


(10) 
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So  the  shear  stresses  in  ice,  in  flight  and  in  the  wind  tunnel  and  the  same  in  similar  sections. 

This  allows  one  to  hope  that  natural  ice  sheddings  due  to  centrifugal  forces,  as  obtained  in  the  wind 
tunnel,  are  well  representative  of  those  in  flight.  In  fact,  the  natural  sheddings  of  ice  are  certainly  under  the 
strong  influence  of  the  blade  flexibility.  In  the  case  of  the  helicopter  model  used  in  the  tests  presented  below, 
the  model  blades  were  stiffer  than  those  of  the  real  rotor.  The  wind  tunnel  test  results  are  thus  pessimistic,  as 
ice  will  tend  to  shed  more  easily  in  flight. 


2 - TESTING  MEANS  AND  TECHNIQUES 


The  specifications  of  the  standards  concerning  icing  conditions  in  flight  required  the  fabrication  of  a set-up 
providing  a wall  defined  icing  cloud  in  test  section  n°  1 of  the  S1-MA  wind  tunnel. 

The  pulverization  grid  can  comprise  444  air  fractioning  injectors  distributing  an  homogeneous  cloud  of  droplets 
whose  median  volume  diameter  can  be  adjusted  between  10  and  20  /jm  over  a 4-m2  surface.  The  liquid  water  content 
can  be  adjusted  between  0.4  and  10  g/m3  for  velocities  reaching  100  m/s. 

The  air  and  water  supplies  to  that  grid  can  be  automatically  controlled  so  as  to  ensure  alternate  icings  (condi- 
tions of  maximum  continuous,  maximum  intermittent  and  clear  sky,  according  to  official  specifications)  while  keeping 
the  droplet  diameter  practically  constant,  with  periods  varying  from  3 seconds  to  several  tens  of  seconds. 


Preliminary  tests  had  shown  that  the  droplets  evaporated  if  the  air  of  the  main  flow  was  not  saturated.  This 
evaporation  modified  the  droplet  diameter  and  the  liquid  water  content  of  the  artificial  cloud,  which  thus  caused  ice 
deposits  not  representative  of  the  set  conditions.  To  palliate  this  difficulty,  the  S1-MA  wind  tunnel  was  equipped  with 
a saturating  device  made  of  8 pulverization  nozzles  providing  up  to  9C0  dm3/h  each,  and  thus  ensuring  saturation  in 
the  test  section  by  raising  the  hygrometric  degree  from  60  to  nearly  100%.  The  permanent  monitoring  of  the  hygro- 
metric  degree  is  ensured  by  an  automatic  dew-point  hygrometer  operating  by  continuous  pick-up  in  the  test  section 
(fig.  3). 
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Fig.  3 - Icing  and  saturation  davicat. 
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The  pulverization  grid  can  be  placed  in  various  points  of  the  test  section  according  to  the  shape  and  position 
of  the  model  (figs.  4,  5,  6). 
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The  position  of  the  pulverizing  grid  relative  to  the  model  is  determined  after  a preliminary  test,  during  which  the 
fluid  filets  surrounding  the  model  are  visualized  by  smoke  streaks  ; this  positioning  is  later  checked  with  the  water  dro- 
plets. 

Tests  made  it  possible  to  determine  the  positive  or  negative  pressure  drop  in  the  wind  tunnel,  behind  the  grid, 
as  a function  of  air  flow,  and  to  take  it  into  account  for  the  determination  of  the  velocity  at  infinity  upstream. 

The  study  of  the  aerodynamic  flow  behind  the  grid  was  carried  out  by  combing  and  by  comparing  wall  visua- 
lizations, pressure  distributions  and  aerodynamic  weighings  on  the  profiles  with  and  without  the  grid.  The  perturbating 
effect  of  the  grid  is  very  low  ; an  example  of  this  effect  is  given,  in  the  case  of  a rotor,  on  figure  7.  The  curves  on  the 
left  present  the  results  obtained,  in  the  same  operating  conditions,  with  and  without  grid  in  the  test  section.  The  curves 
i on  the  right  present  another  experiment,  where  the  rotor  has  been  weighed  behind  the  grid  with  and  without  air  injec- 

tion. The  differences  are  small,  of  the  same  order  of  magnitude  as  the  precision  of  the  measuring  instruments  used  for 
this  experiment.  We  must  then  conclude  that  the  presence  of  the  icing  grid,  in  front  of  the  rotor,  does  not  cause  any 
appreciable  perturbation  in  the  operation  of  the  latter. 


Ftg.  7 — Efhctt  of  tht  grid  and  air  Injection  on  rotor  coafflclanti. 
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The  8-m-dia.  S1-MA  test  section  permits  the  experimentation  of  large  models  with  a low  blockage  factor  ; these 
models  can  easily  be  reached  five  minutes  after  the  test,  which  makes  possible  the  recording  of  ice  deposit  extent  and 
shape  for  photography  and  by  direct  measurements.  Interpretation  of  the  photographs  is  made  easier  by  the  use  of 
cross-lined  patterns  or  of  very  precise  marks  painted  on  the  models,  identifying  chord  percentages  or  angles  (fig.  8). 

In  case  of  important  ice  deposits,  their  castings  can  be  performed  (fig.  9).  The  tunnel  floor  comprises  a balance  allowing 
an  angular  setting  and  the  measurement  of  aerodynamic  forces  on  the  model  before  and  after  icing,  as  well  as  a com- 
pressed air  supply  to  feed  with  primary  air  the  ejectors  simulating  the  engine  air  intakes  in  some  models. 


Fig.  8 — Profile  of  an  ice  deposit  measured 
with  a 1 cm-square  pattern. 


The  test  section  sidewalls  can  be  equipped  with  floodlights  and  high-speed  cameras  (1200  frames  per  second) 
necessary  for  studying  the  trajectory  of  ice  chunks  shed,  either  naturally  or  artifically,  from  the  forward  parts  of  the 
model. 


The  particular  test  rig  for  helicopter  rotors  is  described  in  detail  in  ref.  [8,  9,  10].  This  set  is  put  into  operation 
under  the  control  of  a monitoring  system  ensuring  the  safety  of  personnel  and  installations.  This  monitoring  and  safety 
system  comprises  : 
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- the  visualization  in  two  large  multitrace  screens  of  the  signals  provided  by  strain  gages  mounted  cn  the  blades 
and  the  signals  from  the  balance  and  the  torquemeter  ; the  main  parameters  are  permanently  recorded  on  an  analogue 
magnetic  tape  ; 

- an  independent  electric  supply  ensuring  the  operation  of  some  essential  functions  (lubrication,  pitch  control) 
in  case  of  breakdown  of  the  main  electric  network  ; 

; - a closed  television  circuit,  with  stroboscopic  lighting  and  magnetic  recording  for  the  permanent  monitoring  of 

the  rotor  ; 

- a particular  interphone  between  the  rotor  control  room  and  the  monitoring  room  ; order  transmission  and 

p exchange  of  essential  safety  informations  are  thus  ensured  without  any  risk  of  interference  with  the  routine  test  com- 

munications. 

Let  us  add  that  the  test  section  and  its  windows  are  reinforced  and  that,  before  each  test,  a thorough  study  of 
the  dynamic  characteristics  of  the  test  rig  and  the  blades  is  carried  out,  this  study  being  liable  to  lead  to  modifications 
of  the  natural  frequencies  of  parls  of  the  rig. 

The  acquisition  unit  allows  on-line  calculations,  for  the  setting  of  the  test  parameters  as  well  as  of  the  results 
on  aerodynamic  weighings  ; these  can  be  presented  on  a cathode  ray  tube  as  a curve. 

The  recording  on  magnetoscope  and  by  photography  of  the  observations  made  under  stroboscopic  lighting 
allows  the  restitution  of  ice  accretions. 


3 - VALIDITY  OF  ICING  AT  REDUCED  SCALE 

Having  yet  no  element  for  comparing  icings  in  flight  and  at  reduced  scale  in  the  wind  tunnel,  it  has  been  consi- 
dered necessary  to  check  the  calculations  concerning  icing  in  similitude  by  comparative  wind  tunnel  tests  at  full  and 
reduced  scale.  So  ONERA  undertook,  in  parallel  with  industrial  tests,  similitude  icing  tests  either  on  simple  shapes 
such  as  cylinders  or  on  wing  shapes  similar  to  those  concerned  with  industrial  tests.  These  tests  were  performed  first 
at  1/6  scale,  then  extended  to  1/12  scale.  In  this  case,  the  calculations  are  the  same  as  for  flight-wind  tunnel  similitude, 
but  as  it  is  not  possible  in  the  wind  tunnel  to  adjust  the  pressure  and  temperature,  the  theoretical  thermodynamic 
icing  similitude  can  not  be  strictly  obtained. 

In  the  cases  where  the  icing  similitude  laws  were  not  respected,  the  ice  deposits  appearing  on  the  models  are  very 
dissimilar  as  regards  both  extent  and  shape.  Figure  10  presents  a set  of  photographs  in  which  the  left  hand  ones  (a) 
represent  icings  obtained  on  models  associated  with  correct  clouds,  and  the  right  hand  ones  (b)  models  associated  with 
clouds  not  respecting  the  similitude  laws.  Icings  on  cylinders  (fig.  10-1 -a  and  b)  show  very  different  glaciation  angles 
(80  and  140°).  As  for  ice  thicknesses,  brought  back  to  the  same  scale  (full  scale),  they  are  also  widely  different  : 30 
and  180  mm.  On  delta  wings  (fig.  10-2-a  and  b),  the  icing  on  the  leading  edge  that  appears  in  normal  test  is  practically 
non  existent  on  the  large  model.  In  the  case  of  icing  on  segments  of  thick  profiles  (fig.  10-3-a  and  b),  the  differences 
are  still  more  dramatic  : the  correct  icing  extends  up  to  2.5%  of  the  chord  while  on  the  small  model  it  extends  up  to 
15%,  with  very  different  shapes  of  accretions. 
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Icings  at  small  scale  are  more  regular  than  those  at  large  scale,  when  lateral  protuberances  appear  ; on  the  other 
hand,  the  zone  localised  around  the  stagnation  line  of  the  flow  is  made  of  transparent  ice  in  both  cases. 


Figures  11  and  12  show,  as  photographs  and  scaled  sketches,  results  of  tests  on  cylindrical  models  in  1/6  and  1/12 
scale  similitude.  The  glaciation  angles  are  very  similar,  and  the  thicknesses,  brought  back  to  the  some  scale,  are  not 
very  different,  as  shown  on  Table  I : 


Table  I 


— Results  of  similitude  tests 


Fig.  10 


— Scale  influence  on  icing  - Tests  with  the  same  cloud  and  models  at  scales  1 and  1/6. 
al  Large  model  - Similitude  respected 
b)  Small  model  - Similitude  not  respected. 

1.  Cylinder  2.  Delta  wing  3.  Swept  leading  edge 
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Fig.  11  - Comparison  of  ice  deposits  on  cylinders  in  similitude  conditions. 


Fig.  12  - Comparison  of  lea  deposits  on  cylinders  in  similitude  conditions. 
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Figure  13  shows  the  shapes  and  extents  of  ice  deposits  at  full  and  1/6  scale  on  two  half  delta  wings  of  75° 
sweep  and  at  15°  angle  of  attack.  On  these  wings,  the  ice  appears  only  on  the  leading  edge  lower  surface,  as  a corru- 
gated deposit  whose  section  perpendicular  to  the  leading  edge  has  the  shape  of  an  isosceles  triangle  attached  by  one 
of  its  apexes.  This  shape  is  observed  at  both  scales. 

Other  comparisons  have  been  made  on  thick  profiles  at  incidence  [11 1. 

The  whole  set  of  experiments  made  in  conditions  respecting  the  similitude  laws  shows  a good  agreement  as 
regards  both  deposit  shapes  and  extents  ; however,  a slightly  smaller  deposit  is  generally  observed  at  small  scale,  which 
may  be  explained  by  the  formation  of  a small  percentage  of  ice  crystals  at  the  grid  outlet.  This  percentage  could  be 
higher  in  tests  at  small  scale,  which  requires  smaller  droplets  and  a higher  injection  air  pressure. 

More  recently,  comparisons  between  similitude  tests  in  the  wind  tunnel  with  full  scale  tests  of  an  SST  air 
intake,  and  with  actual  icing  in  flight  (Concorde,  Airbus),  provide  an  essential  complement  to  the  validity  of  the 
method.  Moreover,  it  must  be  noted  that  the  natural  shedding  of  ice  chunks,  due  to  the  stratified  nature  of  the 
deposits,  has  been  also  observed  on  Concorde  in  flight. 


at  Scale  1 b)  Scale  1/6 


Fig.  13  - Companion  of  ica  depotitt  on  delta  wings  in  timilitude  conditions. 
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4 - TESTS  PERFORMED 


Industrial  tests  performed  to-date  in  the  S1-MA  wind  tunnel  concern  aircraft  parts  at  full  scale  and  aircraft 
complete  models  at  reduced  scale. 

4.1.  Tests  at  full  scale 

The  case  of  icing  at  full  scale  can  be  considered  as  similitude  scale,  with  /r  • 1.  This  it  is  only  a particular  case 
of  similitude  tests. 

The  tests  are  usually  performed  on  actual  aircraft  parts,  either  to  test  de-icing  devices  or  to  define  the  tones 
and  shapes  of  accretions  for  non  de-iced  parts. 
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Figure  14  shows,  at  the  end  of  a test,  the  leading  edge  of  an  aircraft  fin  equipped  with  a pneumatic  de  icing 
device.  A non  de  iced  reference  2one  carries  the  ice  accumulated  during  the  whole  test  ; the  rest  of  the  leading  edge, 

completely  free  of  ice,  demonstrates  the  efficacy  of  the  de-icing  device.  ( 

Figure  15  shows  the  evolution  of  ice  deposition  in  the  central  part  of  non  de  iced  fins  of  an  aircraft  during  tests  ( 

at  Mach  0.3  lasting  10,  20  and  30  minutes.  During  these  tests,  the  drag  rise  has  been  measured. 

Figure  16  shows  the  effect  of  the  actuation  of  the  horizontal  control  surface  of  a fin  after  a 30  mn  icing  test. 

During  the  movement,  the  effort  required  to  break  the  ice  has  been  measured. 

Other  tests  have  been  performed  with  thermal  de-icing  devices. 


Fig.  IS  - evolution  of  ica  dapoiitt  on  a tail  tin  during  an  icing  tait. 
al  altar  10  mn  b)  aftar  20  mn  cl  aftar  30  mn 
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Fig.  16  - Actuation  of  a horizontal  fin  control  surface, 
al  before  actuation  bl  after  actuation 

4.2.  Tests  at  reduced  scale 

These  tests  concern  wings  and  fins,  as  well  as  rotors. 

4.2. 1.  Wing  models 

The  study  of  icing  on  Concorde  has  been  performed  in  cooperation  with  the  British  and  French  Official 
Services,  on  a 1/6  scale  half  model  built  by  the  British  Aircraft  Corporation.  These  tests  aimed  at  the  study  of  the 
volume  and  extent  of  ice  deposits  likely  to  occur  on  the  upper  surface  and  the  leading  edges  of  the  wing  and  air 
intakes  during  stand-by  and  approach  flights.  The  tests  revealed  that  upper  surface  deposits  were  practically  non  exis- 
tent, while  the  deposits  on  the  wing  leading  edge  were  limited  to  the  leading  edge  lower  surface,  with  a stratified 
shape  promoting  natural  sheddings  by  rupture  (fig.  17,  18).  These  very  brittle  deposits,  as  well  as  the  natural  sheddings 
whose  trajectories  had  been  filmed  in  the  wind  tunnel,  have  been  observed  during  Concorde  flight  tests.  Figure  18 
also  shows  the  influence  of  ice  deposits  on  the  longitudinal  aerodynamic  coefficients. 

Similitude  icing  tests  have  been  performed  on  a 1/12  scale  model  of  Airbus  ; these  tests  revealed  the  positions 
of  accretion  zones  on  the  wing  and  fin  leading  edges  and  on  the  engine  pod  lips  (fig.  19).  A flight-wind  tunnel  compa- 
rison on  the  shape  of  wing  leading  edge  deposits  is  given  on  figure  20  : the  photograph  of  the  flight  tests  is  taken  in 
the  Laschka  and  Jesse  paper  [12].  Both  in  flight  and  in  the  wind  tunnel  there  appears  a stratified  ice  shape,  called 
"lobster  tail",  that  characterizes  icing  on  a swept  leading  edge. 


Fig.  17  - Ice  deposit  on  a 1/6  scale 
half  model  of  Concorde. 
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4.2.2.  Rotor  models  <, 
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Durirg  the  1973-74  winter,  ONERA  performed,  in  cooperation  with  the  French  official  services  and  the 
Aerospatiale  Company,  the  first  similitude  icing  tests  on  a helicopter  rotor  1 13].  These  tests  were  performed  on  a 
4 m-dia.  3-blade  rotor,  supplied  by  Aerospatiale  ; the  blade  chord  was  210  mm  and  the  root  profile  a NACA  0012. 

The  hub  allowed  a collective  pitch  adjustable  from  0.9  to  26.5°,  as  well  as  a cyclic  pitch  setting.  The  rotor  blades 
were  made  of  a glass  fiber  reinforced  resin,  with  a metallic  longeron. 

The  test  technique  consists  in  adjusting  to  their  nominal  values  the  rotor  and  the  wind  tunnel  before  icing 
injection  of  icing  cloud.  As  soon  as  these  settings  are  obtained,  the  cloud  in  injected  and  the  rotor  rotation  speed  is 
maintained  constant  by  action  on  the  driving  motors.  During  the  whole  icing  test,  aerodynamic  weighings  and  magne 
toscope  recordings  are  performed.  At  the  end  of  the  test,  photographs  of  ice  deposits  are  taken  in  the  test  section. 


An  example  of  these  photographs  is  presented  on  figure  21.  The  four  parallel  white  lines  are  located  at  5,  10, 
20  and  30%  of  the  chord  respectively.  The  two  white  lines  numbered  3 and  5 are  respectively  at  0.3  and  0.5  R.  The 
photographs  of  figure  21  show  that  ice  deposits  appear  in  the  central  zone  of  the  rotor,  close  to  the  hub,  and  that 
the  blades  are  free  of  ice  for  any  relative  radius  above  0.5  R.  This  absence  of  ice  may  be  explained  by  the  increase 
in  kinetic  heating  along  the  span,  due  to  that  of  the  local  aerodynamic  velocity  /(/♦  U Rl,  and/or  by  the  increase  of 
centrifugal  force  along  the  span. 


Fig.  21  - Ice  deposits  on  the  hub  and  blade  upper  surface. 


The  photographs  of  figure  22  show  the  three  views  defining  the  extents  of  ice  deposits  near  the  root  of  one  of 
the  blades.  The  views  of  the  lower  and  upper  surfaces  show  that  the  ice  deposit  on  the  leading  edge  is  localized  in  the 
first  5%  of  the  chord.  In  this  zone,  the  maximum  thickness  is  about  5 mm,  corresponding  to  8.5  mm  at  the  full  scale 
of  a 6.7-m-dia.  rotor  (k  - 0.6).  There  is,  especially  on  the  upper  surface,  a rather  thin  deposit  of  transparent  ice  provo- 
ked by  the  residue  of  water  picked  up  at  the  leading  edge,  not  transformed  into  ice  at  impact  and  flowing  under  the 
effect  of  aerodynamic  and  centrifugal  forces  ("run-back"  phenomenon).  The  test  temperatuie  was  not  very  low 
(-4°C),  and  the  fraction  of  water  transformed  into  ice  at  impact  is  relatively  small,  so  there  is  an  important  percen- 
tage of  water  able  to  run  along  the  surfaces. 
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The  leading  edge  deposit  is  not  very  smooth  ; there  appear  on  the  photographs  of  figure  22  some  dents  (marked 
by  arrows)  corresponding  to  ice  sheddings  that  occurred  during  actual  icing.  Visual  stroboscopic  observations,  as  well 
as  restitutions  from  the  magnetoscope,  revealed  that  after  an  icing  duration  of  some  five  minutes  during  which  ice 
settles  regularly  on  the  blade  leading  edge,  cycles  take  place  during  which  ice  settles,  grows  and  is  torn  off.  The  period 
of  such  cycles  is  difficult  to  determine,  but  is  somewhere  between  10  and  20  seconds,  corresponding  to  20  and  40 
seconds  at  full  scale.  The  icing  test  could  last  a very  long  time  without  the  cumulated  amount  of  ice  becoming  more 
important  than  during  a single  cycle  ; as  an  average,  there  is  as  much  ice  deposited  as  evacuated  by  centrifugation. 


Lower  surface 


Leading  edge 


Upper  surface 


Fig.  22  - Ice  deposits  on  a rotor  blade. 


Fig.  23  - Variations  in  time  of  rotor  coefficients. 


The  cycle  phenomenon  appears,  as  regards  aerodynamic  weighings  of  the  rotor,  as  constant  values  after  5 to  6 
minutes  of  icing  (see  fig.  23)  : on  the  three  main  parameters,  i.e.  lift,  propulsive  thrust  and  torque,  there  is  a rather 
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rapid  evolution  from  the  icing  onset,  then  a tendency  to  stabilize,  when  the  above  cycling  phenomenon  takes  place. 

This  stabilization  is  obtained,  in  the  test  case  mentioned  here,  with  a lift  drop  of  2.5%,  a rise  in  the  torque  required  ( 

from  the  constant  speed  motor  of  25%,  and  a decrease  in  propulsive  thrust  of  65%.  During  the  whole  icing  duration, 
the  monitoring  of  the  strain  gages,  and  also  of  the  dynamic  values  of  the  signals  provided  by  the  dynamometric  unit,  f 

did  not  show  any  significant  increase  as  regards  vibrations. 

No  flight  test  resu'ts  are  yet  available  that  are  strictly  comparable  to  these  performed  in  the  wind  tunnel,  but  it  is 
possible,  as  regards  the  comparison  on  the  shape  and  extent  of  ice  deposits  and  on  the  evolutions  of  aerodynamic 
parameters  during  flight  in  icing  conditions,  to  confer  to  documents  ref  [14  through  19],  In  [14],  Stallabrass  shows 
that  there  is  a limit  to  ice  extension  on  the  blade  leading  edge,  which  depends  on  external  static  temperature  and  is 
due  to  kinetic  heating.  Still  in  [14],  the  author  shows  that  there  is  another  limit  of  ice  thickness,  due  to  centrifugal 
shedding  of  ice  chunks.  These  theoretical  considerations  on  the  spanwise  limitations  are  confirmed  by  some  flight 
tests  presented  in  [15  through  18].  Icing  behind  the  leading  edge  zone,  due  to  the  "run-back"  phenomenon  (see  fig. 

22),  is  also  found  in  flight  [14,  15].  The  results  obtained  in  flight,  as  described  in  ref.  [ 1 4,  15,  17,  19], show  increases 
in  the  torque  and  power  required  from  the  engine  during  flights  in  icing  conditions  quite  comparable  to  those  mea- 
sured in  the  wind  tunnel  (30%  increase  in  torque),  as  well  as  the  appearance  of  a stabilization  of  the  torque  when  the 
cyclic  phenomenon  of  centrifugal  ice  shedding  is  reached.  Lake  and  Bradey  [19]  mention  the  great  sensitivity  of  the 
rotors  to  small  ice  deposits,  which  may  be  explained  by  the  very  fast  increase  of  the  torque  required  from  the  driving 
motor  in  the  wind  tunnel  from  the  very  onset  of  icing  conditions. 

This  set  of  comparisons,  though  performed  on  rotors  different  from  that  tested  in  the  wind  tunnel,  show  such 
similarities  that  there  is  no  doubt  that  similitude  icing  tests  in  the  wind  tunnel  give  a very  close  representation  of  tests 
in  actual  flight. 

5 - EXTENSION  POSSIBILITIES  OF  THE  TESTS 

The  large  SI  wind  tunnel  of  Modane  used  for  icing  tests  is  cooled  by  aic  exchange  with  the  atmosphere,  and  can 
be  used  to  this  end  only  during  cold  weather  in  the  winter  months.  This  time  limitation  in  the  general  planning  of  the 
wind  tunnel  operation  is  thus  complemented  by  a random  limitation  due  to  meteorology.  To  palliate  these  difficulties, 

ONERA  is  studying  means  for  artificially  cooling  the  test  section.  A project  combining  a low  temperature  exchanger 
and  expansion  of  air  at  high  pressure  should  provide  the  frigories  necessary  to  reduce  the  air  temperature  in  the  wind 
tunnel  by  15°C.  The  realization  of  that  project  should  permit  us  to  ignore,  for  a larqe  number  of  tests,  the  local  atmos- 
pheric conditions,  and  to  perform  tests  at  very  low  temperature  necessary  to  qualify  some  devices  such  as  de-icing  and 
anti-icing  systems. 

6 - CONCLUSION 

The  development  of  artificial  icing-cloud-producing  devices  in  large  wind  tunnels  on  the  one  hand,  and  that  of 
theoretical  laws  of  similitude  applicable  to  icing  problems  on  the  other  hand  made  it  possible  to  perform  in  wind 
tunnels,  and  in  particular  in  the  ONERA  SI  tunnel  of  Modane,  icing  tests  that  provided  quite  significant  results. 

These  results  concerned  the  qualitative  pile-up  of  ice  deposits  on  aircraft  parts  at  full  scale  and  of  complete  air- 
craft models  at  reduced  scale,  down  to  1/12  ; similar  results  were  obtained  on  much  more  complicated  structures  such 
as  helicopter  rotors.  Measurements  taken  during  these  tests  provided  quantitative  results  quite  comparable  to  those 
obtained  in  flight. 

It  is  thus  possible  to  use  in  a practical  way  the  wind  tunnel  to  predict  the  performance  of  new  fixed  or  rotary 
wing  aircraft  from  wind  tunnel  tests.  Future  technological  developments  should  increase  the  yearly  amount  of  such 
tests. 
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ICHC  TEST  FACILITIES  IN  CANADA 
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STMURT 

Tha  icing  simulation  faoilltlaa  for  raaaaroh,  development  and  testing  at  tha  National  Raaaarch  Council 
of  Canada  arc  described.  The  major  faoilltlaa  inoluda  a low  speed  icing  wind  tunnel,  a high  speed  variable 
density  tunnel,  a refrigerated  Jet  angina  teat  oall  and  a helicopter  icing  spray  rig.  An  Appendix  is 
included  that  Hats  reports  and  paper*  published  on  icing  slnoa  1935. 


INTRODUCTION 

Icing  raaaaroh  and  taat  facilities  ware  in  us*  in  Canada  during  the  early  UOs,  although  at  that  time 
■oat  of  tha  work  waa  conducted  during  flight  testing  and  only  Halted  ground  taat  facilities  ware  availa- 
ble. The  construction  of  an  open  circuit,  naturally  refrigerated,  icing  wind  tunnel  to  be  powered  by  an 
aircraft  engine  and  propeller  waa  carried  out  between  191*2  and  191*5,  although  it  waa  never  placed  in  serv- 
ice. Ultimately  tha  building  containing  this  open  oirouit  tunnel  became  tha  nucleus  of  a much  larger 
building  designed  for  icing  and  cold  weather  raaaarch. 

The  NEC  Low  Taaperature  Laboratory  of  tha  Mechanical  Engineering  Division  waa  opened  in  191*8^.  The 
physical  plant  of  this  laboratory  consisted  of  a large  refrigeration  system  in  addition  to  three  cold 
chambers  and  an  icing  wind  tunnel.  The  facility  that  subsequently  has  beooae  known  aa  tha  Low  Spaed  Icing 
Wind  Tunnsl  waa  originally  placed  in  operation  in  1950.  In  the  early  60s  a high  apaad  lolag  tunnel  waa 
Installed  in  an  extension  to  tha  Low  Taaperature  Laboratory  building. 

In  the  late  l*0s  the  Engine  Laboratory  of  tha  NBC  Mechanical  Engineering  Division  operated  a field 
station  at  Churchill,  Manitoba,  in  which  Jet  anginas  could  be  tasted  under  natural  low  taaperature  condi- 
tions. Icing  of  Jet  anginas  was  first  carried  out  at  tha  Port  Churchill  location.  Ultimately  icing  teats 
were  transferred  to  Ottawa  in  one  or  aor*  of  the  angina  teat  oalls  at  tha  Engine  Laboratory. 

In  tha  early  50a  following  work  on  propeller  and  wing  protection  for  fixed  wing  aircraft,  tha  Low 
Taaperature  Laboratory  undertook  oonet ruotion  of  aa  lolag  research  facility  in  which  a haliooptar  could  bo 
hovered.  While  the  original  design  waa  not  satisfactory  tha  redesigned  haliooptar  spray  rig  has  now  bean 
operated  for  over  twenty  years  on  two  altos. 

In  addition  to  the  loing  facilities  of  tha  Low  Taaperature  Laboratory  and  tha  Engine  Laboratory,  the 
Fuels  and  Lubrloante  Laboratory  has  conducted  tests  on  aircraft  angina  carburetor  icing  with  a specialised 
faoility(2)\3).  Outside  of  tha  Rational  Raaaarch  Counoil  there  exists  one  other  icing  wind  tunnel  in  ,,  . 
Canada*  the  Uhl  varsity  of  Alberta  that  ha*  bean  used  primarily  for  non-ae  ronaut  1 cal  icing  investigations'  . 

LOW  SPEED  ICING  WIND  TUNNEL 

this  is  a dosed  return  oirouit  tunnel  vented  to  atmosphere  at  tha  working  section^  The  layout 
of  the  tunnel  with  a number  of  representative  croaa  aeotiona  is  shown  in  figure  1.  Tha  tunnel  ahell  is  of 
steel  plat*  reinforced  with  steel  riba  and  supported  on  concrete  piers.  The  fan  section  is  the  only  fixed 
portion  while  the  remainder  of  tha  tunnsl  slides  on  tha  supports  to  allow  for  expansion  and  contraction, 
figure  2 shows  the  tunnel  ahell  construction. 

A 1000  HP  synchronous  speed  eleotrlo  motor  driving  through  a variable  speed  electro -magnetic  ooupling 
provides  power  to  an  eight-blade  fan.  The  drive  assembly  is  shown  in  figure  3* 

The  air  stream  is  ooolad  by  direct  evaporation  of  f onis  in  the  bank  of  oooling  coils  loos ted  in  the 
lowest  velocity  son*  of  the  tunnel.  The  refrigeration  system  is  a three-stag*  ammonia  compression  plant 
of  1000  HP  oapaelty.  At  marlmim  air  speed  sod  under  summer  oondansing  temperatures,  a m1n1m>»  temperature 
of  -20  *C  can  be  obtained.  At  lower  speeds  temperatures  of  -1*0  *C  oan  be  obtained. 

Air-atomising  water  sprays  looated  upstream  of  the  working  seotion  produce  the  cloud  water  droplets. 

A liquid  water  oontent  up  to  3 gm/m3  oan  be  obtained.  Median  droplet  diameters  down  to  15  microns  oan  be 
obtained.  Tap  water  is  used  to  supply  the  sprays  of  this  tunnel.  The  oloud  liquid  water  oontent  is  con- 
trolled by  the  supply  of  water  to  the  nossles  which  is  monitored  by  visual  indicating  flowmeters. 

The  Low  Speed  loing  Tunnel  was  designed  with  a mm  bar  of  features  to  allow  for  loing  that  re duo*  the 
aerodynamio  effectiveness  and  result  in  a low  energy  ratio.  Between  the  diffuser  and  the  fan  a large  low 
veloolty  settling  dumber  was  provided  so  that  iolng  of  the  fan  would  not  result.  The  divergent  seotion 
between  the  fan  and  the  oooling  coll  section  has  a greater  angle  than  desirable  because  of  spao*  limitations. 

This  tunnel  has  been  used  for  a wide  variety  of  iolng  research  and  test  programs  for  the  past  quarter 
oenturyw).  in  the  early  50*  the  research  and  test  work  on  the  electro- thermal  de-icing  system  for  the 
wings  and  e^enneg*  of  the  CP  100  was  conducted.  The  iolng  of  engine  inlets  of  various  types  has  been 
investl^ted  a a is  indioated  by  figure  U.  The  iolng  of  swept  wirg’s  and  delta  wing  shapes  waa  investigated. 
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The  performance  of  hot  gas  leading  edge  protection  has  been  Investigated  In  this  tunnel  since  auxiliary 
compressed  air  is  available  at  the  test  section  from  sources  in  other  laboratories. 

Non-aeronauti cal  icing  research  on  electrical  transmission  lines > communications  equipment  and  ship 
components  has  been  carried  out.  Provision  has  been  made  to  simulate  freezing  rain  in  the  working  section 
of  the  tunnel  for  some  of  the  non-aeronautlcal  icing  work. 

HIGH  SPEED  ICING  WIND  TUNNEL 

This  wind  tunnel  is  a dosed  oircuit  variable  density  icing  tunnel.  The  vertical  cross-section  of 
the  tunnel  is  shown  in  Figure  £.  It  is  powered  by  a 600  HP  D.C.  motor  driving  a two-stage  compressor 
through  a speed  increaser.  The  original  design  provided  maximum  velocities  up  to  Mach  0.9  in  a 1-ft. 
square  by  18-in.  long  working  section.  Recently  a semi -permanent  insert  ahead  of  the  working  section  has 
been  installed  that  limits  speeds  to  approximately  300  knots  while  improving  the  low  speed  icing  simula- 
tion. 


The  tunnel  shell  is  a sealed  and  reinforced  steel  structure,  thus  permitting  a variation  in  the  static 
pressure  at  the  working  section  from  sea  level  to  the  equivalent  of  30,000  ft.  altitude  by  application  of 
a 150  cfm  vacuum  pump.  This  tunnel  is  refrigerated  by  the  same  three-stage  ammonia  plant  as  the  low  speed 
tunnel,  thus  simultaneous  operation  of  the  two  tunnels  is  not  possible.  The  high  speed  tunnel  is  cooled 
using  a secondary  refrigerant,  trichloroethylene,  pumped  through  a finned  tube  heat  exchanger.  Tempera- 
tures as  low  as  -]*0#C  can  be  attained  in  the  tunnel  with  temperature  control  in  the  range  of  +0.5*C. 

Icing  1b  produced  in  the  working  section  by  pneumatic  atomized  water  droplets  introduced  into  the 
air  stream  by  nozzles  upstream  of  the  working  aeotlon.  The  ioing  simulation  can  be  varied  over  concen- 
trations up  to  2 gm/m3  with  a droplet  size  of  20  microns. 

Due  to  the  small  size  of  the  working  section  of  this  tunnel,  it  has  been  employed  primarily  for  ioing 

instrumentation  research,  calibration  and  development  tests.  It  has  been  ub ed  to  test  high  speed  aircraft 

components  subject  to  icing  only  to  a limited  extent. 

HELICOPTER  ICING  SPRAY  RIG 

In  1953  interest  in  helicopter  icing  arose  as  a result  of  a request  from  the  Canadian  naval  services. 

A prototype  spray  rig  described  in  detail  in  NRC  publication  453°'° J was  operated  in  the  1953-54  winter. 

A second  spray  rig  was  designed  and  constructed  during  1951*  that  subsequently  has  seen  modified  and  moved 
to  a second  test  site™).  This  sprmy  rig  has  been  operated  since  the  winter  of  1954-55  with  a somewhat 
variable  utilization  eaoh  winter. 

Hie  test  site  of  this  icing  facility  is  at  the  northwest  extremity  of  Uplands  International  Airport 
in  Ottawa.  As  shown  in  Figure  6,  it  has  access  to  one  runway  and  thus  some  fixed  wing  aircraft  have  been 
able  to  taxi  onto  the  teat  site  for  ioing  tests,  although  the  road  was  provided  to  tow  out  unserviceable 
helloopters. 

The  spray  rig  provides  a cloud  75  feet  wide  by  15  feet  deep  from  a nozzle  array  of  161  steam  atomiz- 
ing water  nozzles  mounted  on  a steel  frame.  The  steel  frame  can  be  positioned  either  at  the  bottom  of 
the  supporting  mast  or  in  the  uppermost  position  so  that  the  aloud  is  approximately  5°  feet  above  ground 
level.  A wind  velocity  of  approximately  10  mph  is  required  to  move  the  cloud  away  from  the  array  and 
provide  an  adequate  streaming  cloud  for  a helicopter  to  hover  in  and  have  good  ooverage  of  the  main  and 
tail  rotor. 

The  maximum  L.W.C.  that  can  be  obtained  in  the  spray  rig  is  approximately  0.8  gm/m^.  The  droplet 
size  is  normally  30  microns  median  volvae  diameter.  The  Ambient  temperatures  available  in  Ottawa  permit 
testing  to  -25 ‘C,  although  lower  tumperatures  are  usually  experienced.  The  ioing  season  in  Ottawa  is 
restricted  to  the  tine  interval  between  15  November  and  15  March. 

(a) 

Details  on  the  helicopter  spray  rig  are  given  in  NRC  report  IA-l86A'7/.  Although  this  publication 
describes  the  spray  rig  as  installed  on  its  initial  site,  the  major  change  that  resulted  was  an  increase 
in  the  area  of  the  test  site  and  greater  safety  for  both  the  flight  crews  and  the  rig  operators. 

ENGINE  ICING  TEST  CELL 


Sinoe  the  late  40s  the  Engine  Laboratory  of  the  National  Research  Counoil  has  bean  engaged  in  Jet 
engine  ioing  research,  development  and  testing  to  assist  manufacturer*  and  operators.  Originally  an  ioing 
test  faelllty  was  built  at  Fort  Churchill,  Manitoba.  It  oonsisted  of  a single  straight-through  test  oell 
of  15  x 15-foot  oxoss-seotion.  Hie  facility  was  operated  during  the  winter  months  utilising  the  natural 
refrigeration  available  at  this  looation.  For  icing,  a mobile  water  spray  rig  was  positioned  in  front  of 
the  engine  inlet. 


In  the  early  50s  the  icing  test  work  was  transferred  to  the  main  Engine  Laboratory  in  Ottawa.  Curr- 
ently, number  4 teat  oell  is  equipped  for  ioing  tests  of  turbopropellsr  and  turbojet  ermines' 10 ).  The 
test  oell  has  a 25-foot  diameter  working  section  and  will  accept  engines  up  to  500  lta/a  airflow  and 
50,000  lbf  thrust.  This  teat  cell  has  an  air- trichloroethylene  heat  exohanger  to  oool  the  inlet  air. 

The  refrigeration  capacity  is  suoh  that  a 15*F  depression  can  be  obtained  with  a mass  flow  of  100  lbm/s. 
This  oell  is  operated  from  approximately  1 December  to  31  Maroh  when  temperatures  are  suitable. 

A 3^-nozzle  spray  rig,  looated  in  the  intake  duot, simulates  the  icing  oloud.  To  simulate  ice  crys- 
tals, two  banka  of  stacked  circular  saw  blades  out  ioe  in  blocks  into  partloles  that  are  blown  into  the 
intake  duot  by  centrifugal  fans. 


An  Appendix  is  included  listing  the  unclassified  reports  and  papera  on  ioing  prepared  by  four  labor- 


stories  of  the  NSC  Mechanical  Engineering  Division  during  the  past  four  decades 
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FIG.  4:  EXPERIMENTAL  ENGINE  INLET  WITH  ICING 
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FIG.  6.  HELICOPTER  ICING  TEST  FACILITY 
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J4jAbstract 
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safety  and  performance  and  to  decide  on  the  need  for  anti-or  de-icing  systems.  The 
Fluid  Dynamics  Panel  of  AGARD  sponsored  a round-table  discussion  on  the  subject  of 
Aircraft  Icing  , on  30  September  1977,  in  Ottawa,  Canada,  in  conjunction  with  an  FDP 
Symposium  on  Unsteady  Aerodynamics.  The  seven  papers  presented,  covering  a broad 
spectrum  of  topics,  are  presented  in  this  Advisory  Report. 
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